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WAVE FORCES ON SUBMERGED STRUCTURES 


Ernest F. Brater,! M. ASCE, John S. McNown,2 M. ASCE 
and Leslie D. Stair3 
(Proc. Paper 1833) 


SYNOPSIS 


The magnitude and characteristics of forces resulting from oscillatory | 
waves were determined for models of submerged barge-like structures. The od 2 
wave profiles and the variation with time of the resulting horizontal and - 
vertical forces were obtained for various wave heights, wave periods, and 
locations of the barge with respect to the water surface. The basic model 
was a right parallelepiped having the proportions of a typical barge. Modi- 
fications of this model were made to determine the effects of rounding the 
sides and of the inclusion of a slot in one end of the barge. Horizontal force 


measurements were also made for a flat plate having the same dimensions * | 
as the long vertical face of the barge. These studies provided design data a 
needed for the application of analytical methods in the determination of forces | 
on some types of submerged structures. 


INTRODUCTION | 
The design of off-shore structures such as oil drilling platforms and 
radar stations requires knowledge of the magnitude of the wave forces. 
Methods of calculating such forces have not been fully determined. Labora- 7 
tory investigations have been made of wave forces on vertical piles, (1,2,3,4) F 


horizontal struts(5,6) and on larger and more complicated structures.(7) 


Note: Discussion open until April 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1833 is wa 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 6, November, 1958. 
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No laboratory investigations of forces on simple barge-like shapes have been . 
reported. 
Oe The forces produced by oscillatory waves result from the orbital motion 
of the water. Acceptably accurate estimates can be made of orbital velocities E 
and accelerations, and the associated pressure variations can be determined 3 
wii analytically. Also required for the computation of the drag and inertial forces My 
: are values of the coefficient of drag (Cg) and the coefficient of inertial re- : 
: sistance (Cm). The drag coefficients for various bodies have been determined 2g 
} experimentally for steady flow, but little is known about their magnitudes for " 
z unsteady flow. Values of C,, have been determined theoretically as well as 
experimentally for certain shapes, but here again the applicability of such 
values to the conditions existing during wave motion has not been tested. 
These studies were planned to determine the necessary coefficients and a 
Suitable method for estimating wave forces on structures placed in the sea. 


Testing Facilities 


The experiments were conducted in the Lake Hydraulics Laboratory(8) of 
ie the University of Michigan. The wave tank used for these tests is 86 feet 
- long and 40 feet wide and has walls 3‘ 10°’ high. The testing frame was 
b located in a channel 13' 10°’ wide separated from the remainder of the tank 
by a temporary wall as shown in Fig. 1. Waves were created by a wave a 
% machine of the plunger type located at one end of the channel. The wave oP 
= machine could be regulated to produce waves of various heights and periods. = 
td The model was centered in the channel approximately 45 feet from the 4 
plunger. After passing the model, the waves were reflected out of the channel aa 
a into the larger portion of the tank where their energy was dissipated by means ms 
Pes of an artifial beach (Fig. 1). Wave heights were measured by means of re- 
cording electric resistance gages. 

The testing frame was mounted on two heavy I-beams supported on the 
? walls of the wave channel. A photograph showing the method of supporting 
a the model and some of the instrumentation is shown in Fig. 2. Details of the 
. dynamometer frame are shown in Fig. 3. Two inverted U-shaped frames (U) - 
served as the supporting structure for the models. A horizontal plate (P) : 
was suspended from the U-frames by four slender aluminum rods (r) 4 feet 
long. The dynamometer plate was restrained from horizontal movement by 
“ 3 dynamometers designated in Fig. 3 as hg, hy, and N. The first two pre- 
Be vented the plate from moving in the direction of wave travel, and the third 
j supplied horizontal restraint in the direction normal to wave motion. The 
: models were suspended from the aluminum plate by two struts (s) stream- 
bi. . lined in the direction of wave motion. The length of the struts was adjusta- 
ie ble, so that the vertical position of the model could be varied. 

Dynamometers hg and hj, consisted of small aluminum cantilevers fas- 

tened to one of the I-beams, and connected to the dynamometer plate by 


i means of small bars in such a manner as to be unaffected by vertical forces. : 
Dynamometer N was fastened to a cross member as shown in Fig. 3. Elec- “a 
ye trical resistance gages were cemented on opposite sides of the cantilevers is 
MY and connected to opposite arms of a Wheatstone bridge in a Brush strain am- 4 
le plifier. The dynamometers were calibrated by subjecting them to known 4 
be horizontal forces and recording the corresponding change in resistance by ; 
a means of an oscillograph. The relationship between the force on the dyna- P 
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TESTING ARRANGEMENT 
Fig. 2 


mometers and the force on the model was determined by applying loads at 
various model positions and noting the resulting forces on the dynamometers. 
Vertical forces were measured by means of two dynamometers, vg and 
Vp, Which consited of slender horizontal beams through which the four sup- 
porting rods were connected to the testing frame as shown in Figs. 2 and 3. 
Strain gages were placed on top and bottom of the center of these beams. 
Vertical loads on the model produced bending moments in the beam and 
changes in the resistances of the strain gages which were recorded by means 


of an oscillograph. These dynamometers were so arranged that they were 
not affected by horizontal forces. 


Test Procedure 


An extensive series of preliminary investigations was made before be- 
ginning the principal tests in order to determine the behavior of the waves 
in the testing area. One phase of the preliminary investigation was planned 
to determine the effect of the presence of the model on wave heights in the 
vicinity of the model. The placing of the model in the water caused vari- 
ations in the wave height for some distance upstream. However, a point 11.0 
feet upstream from the model was found to be free from this effect for all 
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conditions, and was, therefore, selected as the wave gage location. 

The variation in wave height in the direction perpendicular to the wave 
motion was also studied. For certain wave frequencies, standing transverse 
waves were developed which produced a variation in wave height from one 
side of the model to the other of as much as 40 per cent. For such conditions, 
a single measurement of wave height at the center of the model provided a , 
value which differed as much as five per cent from the average effective 
wave height. This phenomenon was caused by the establishment of a resonant 
motion in the channel and tank, and no method was found to eliminate it. 
However, these standing transverse waves did not appear until approximately i 
seventeen waves had passed the model following the beginning of a test. Also, 
the first 4 or 5 waves of a group changed in size and length as they progressed 
along the channel from the wave rod to the model. Both difficulties were 
avoided by using only those waves between the 5th and 17th in a series. 

The dynamometers were calibrated by loading the model with known 
weights. The periodic loading resulting from wave motion was simulated by 
manual manipulation of the weights. The calibration coefficient for the sum 
of the two principal horizontal dynamometers differed by 2.5 per cent from 
the beginning until the end of the testing program. For the vertical dyna- 
mometers there was no measurable difference in the calibration coefficient 
from the beginning until the end of the testing program. 

The wave rod was calibrated by raising and lowering the rod known 
amounts and noting the corresponding oscillograph reading. It was necessary 
to conduct the calibration in the wave tank itself because the calibration was 
affected by the chemical composition of the water or the grounding conditions. 
The calibration curve was nearly linear, and, although it remained virtually 
constant, a new calibration was made each day. 

The desired water depth was maintained by adding water to compensate 
for evaporation. The water depth was measured at the front of the model. 

The concrete floor of the tank was relatively smooth but there were undu- 
lations as great as 1/4 inch in the region between the model and the wave 
machine. 

Prior to starting a test run, the wave machine was adjusted to generate 
a wave of the desired height and period, and the model was placed at the 
desired elevation. During the first run simultaneous measurements were 
made of the wave height at a point 11 feet upstream from the model and of 
either the horizontal or the vertical forces. Such a test required three os- 
cillograph channels because both the horizontal and vertical forces were 
obtained by adding the forces in two dynamometers. Upon completing such a 
run, one dynamometer was switched off and a record was obtained from a 
wave rod located at the side of and opposite the center of the model. The 
wave heights measured by this gage were distorted by the presence of the 
model and were used only to obtain a phase relationship between the wave 
and forces. In the early tests, records were also taken from dynamometer 
N to determine forces normal to the direction of wave motion. This was 
Giscontinued because these forces were found to be very small. 


> 


The Experimental Program 


The testing program was planned to determine horizontal and vertical 
forces on barges of several shapes for a variety of wave characteristics and 
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Rectangular Barge Rounded Barge 


Slotted Rectangular Barge Flat Plate 


MODELS 
FIGURE 4 


for various barge locations. Representative results were sought throughout 
so as to establish trends rather than to provide complete results for all 
possible combinations of variables. The depth of water (d) was kept constant 
at one foot throughout the tests. This depth, as well as the model dimensions 
and the wave characteristics, were selected as being typical of representative 
prototype values reduced in a linear scale ration of 1 to 100. 
The basic model was a rectangular parallelepiped 2' 6'' long, 10 inches 
wide, and 1-7/8 inches high as shown in Fig. 4. Two modifications of the 
J basic model were studied. One modification consisted of the inclusion of a 
slot in one end of the barge to simulate a typical drilling structure, and the 
second was the rounding of the upstream and downstream sides as shown in 
Fig. 4. Tests were also made on a flat plate having the same area as one 
side of the barge and a thickness of 1/4 inch. All tests were conducted with 
the top edges of the models horizontal and with the long side (2' 6"' x 1-7/8"') 
normal to the direction of wave motion. 
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For each barge position and wave period, waves of three sizes were pro- 
jected against the model. The smallest wave height (H) in such a series was 
usually about 0.1 foot, the intermediate one approximately 0.2 foot and the 
largest greater than .25 foot. Three wave periods (T) were used in the tests 
on the basic model. These were approximately 1.4 seconds, 1.1 seconds, and 
0.85 second, the corresponding wave lengths (L) being 7.1 feet, 5.2 feet, and 
3.5 feet, respectively. Values of d/L thus varied from 0.14 to 0.29. Actual 
wave periods for individual tests differed slightly, the values given above 
being approximately the averages for the three groups of tests. For the 
period of 1.1 seconds, the model was tested at 10 locations varying from a 
position with the center of the barge 0.17 inches below the water surface to 
a location as near the bottom as possible without touching the bottom. For 
periods of 1.4 seconds and 0.9 second, all tests were conducted with the barge 
at mid-depth. The slotted and rounded barges were tested only for a wave 
having a period of 1.1 seconds. 


Oscillatory Wave Motion 


Oscillatory waves are so called because the water particles oscillate in 
nearly closed orbits, there being substantially no forward motion of the fluid. 
In relatively deep water (d/L >0.5), the orbits are circular at the water 
surface, becoming slightly elliptical and much smaller as the bottom is ap- 
proached. In shallow water (d/L< 0.5), the particle orbits are elliptical 
throughout che depth. The analysis of the test data requires the use of theo- 
retical relationships which define the form and orbital motion of the waves. 
Several rather complete summaries of wave theory are available.(9,10, 11,12) 

Equations for the shape of and motion in oscillatory waves were developed 
by Airy.(13) Although his derivation was based on the assumption that wave 
heights are small, his results are sufficiently accurate for many purposes 
even thorgh used to describe waves with appreciable amplitudes. Stokes(14) 
developed equations which include the effect of wave height, but they are ex- 
ceedingly cumbersome. The small refinement provided by them was not 
justified in the analysis of the test results. Consequently, the relationships 
utilized and presented in the following paragraphs are those credited to Airy. 

The wave length (L), the distance from crest to crest of adjacent waves, 
is related to the other characteristics of a wave by the following expression: 


2 
L 37, tanh 


in which T is the wave period in seconds and d is the depth of water. The 
relationship between wave period, wave velocity (C), and wave length is 


L = CT (2) 


The approximate equation for the wave profile is 


y= cos (3) 


in which x and y are the horizontal and vertical coordinates of the water 
surface and H is the wave height. Values of x are positive in the direction of 
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wave motion and y is measured positively upward from the still water level 
(s.w.l.) as illustrated in Fig. 5. The expression 27x/L is the phase angle 
in radians which is designated as 6. In degrees, @ varies from 0 ata crest 
to 360° at the preceding crest, as shown in Fig. 5. The phase angle can also 
be expressed in terms of time as 27 t/T where t varies from 0 to T during 
a wave cycle. 

The equations for the horizontal and vertical components of the orbital 
velocity are 


7H cosh 27 (d + z)/L 


(4) 
and 
7H sinh 27 (d 
u = 7H Sin ( TA sin 0 (5) 


z T sinh 277 


In these equations, z is the vertical distance from the still water surface to 
the point where the velocity is being determined, upward being taken as posi- 
tive. The corresponding expression for horizontal and vertical accelerations 
can be determined by differentiating Eqs. (4) and (5) partially with respect 

to time. The resulting expressions are 


ou, 27 7H cosh 27 (d + z)/L 

and 

du, -277H 277 (d + 2)/L (7) 

Ot 7 T? sinh 27 d/L 


Oscillatory Wave Forces 


The forces produced on submerged structures by oscillatory waves are of 
two types: the drag force, resulting from the orbital velocity; and the in- 
ertial force, resulting from the orbital acceleration. The drag force can be 
expressed as follows, 


(8) 
d 


in which Cg is a drag coefficient, P is the density, A is the cross-sectional 
area of the structure in a plane normal to the force, and u is the component 
of the velocity in the direction of the force. The flow pattern referred to is 
that which would occur if flow were not disturbed by the presence of the 
model. Because of the nature of the velocity variation, the horizontal drag 
force varies from a maximum downstream value under the crest of a wave 
(x/L = 0) toa maximum upstream value under a trough (x/L = 0.50) as 
illustrated in Fig. 5. 

For steady flow, the drag coefficient has been examined for many different 
shapes and has been found to be a function of the form and of the Reynolds 
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number.(15,16) However, it was not known whether such values could be 
used for the periodic unsteady flow produced by oscillatory waves. An ad- 
ditional complication arises for the case of large structures because of the 
variation of the motion within the dimensions of the body. 

The inertial force is expressed by the equation, 

F, = ot (9) 
in which C,, is the coefficient of inertial resistance, V is the displaced 
volume and 0u/@t is the component of the acceleration in the direction of the 
force for undisturbed flow. This type of force occurs for two reasons. The 
acceleration of the fluid is caused by a pressure gradient so that the pressure 
on the upstream face of the submerged structure is not the same as that on 
the downstream face. The resulting force can be shown to be equal to the 
product of the acceleration and the mass of the displaced volume. Also, 
because of the presence of the body, the accelerations in the flow around the 
body are increased, thus requiring an additional force. This force divided 
by the acceleration of the ambient flow is called the virtual mass, and is 
conveniently expressed by its ratio to the mass of the fluid displaced by the 
body. Thus, C,, is the sum of a term, approximately unity, for the pressure 
gradient force, and the coefficient of virtual mass. 

The variation of the horizontal acceleration, Eq. (6), is such that the 
horizontal inertial force varies from a maximum downstream value at 
x/L = 1/4 to a maximum upstream value at x/L = 3/4 as illustrated in Fig. 5. 

Values of the virtual mass coefficient have been determined theoretically 
and experimentally for bodies of various shapes. For irrotational flow with- 
out separation, values are available for elliptical cylinders and ellipsoids(17) 
and for infinitely long bodies of rectangular cross-section(18) as shown in 
Fig. 6. The corresponding values obtained experimentally, (19,20) by 
oscillations at high frequency and small amplitude, correspond fairly well 
with the theoretical values. Values of C,, obtained in this manner may not 
be applicable to flows in which separation plays a predominant role. In fact, 
the development of a wake has been found to cause marked changes in both 
Cq and Cm.(9) Variations in flow over the length of large bodies cause 
further uncertainties. Consequently, the applicability of theoretical results 
to flows produced by oscillatory waves must be tested experimentally. 

The total force is the sum of the two forces given by Eqs. (8) and (9). 


Because the two forces are 90° out of phase, the maximum combined down- 
stream force may occur at any value of x from x/L = 0 to x/L = 1/4, depend- a 
ing on the relative magnitude of the components. The magnitude and location . 
of the maximum force for the case of cylindrical piles has been discussed 
by Morison and others. (2,3) 

For large structures, it is convenient to express the horizontal inertial a 
force in terms of the pressures exerted on the body. If the rectangular out- 
line shown in Fig. 5 is assumed to contain only the undisturbed fluid, and if 
convective accelerations are neglected, then an application of the equations 
of motion to this volume of fluid in the x-direction yields the following 
expression 


(Py - Pg) Ay (11) 
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Cm FOR TWO DIMENSIONAL FLOW 
PAST RECTANGULAR BODIES 
Fig. 6 


in which p; and pg are the pressures on the two sides of the volume of fluid, 
Ax is the cross-sectional area of the volume taken normal to the force, and 
O ux/@t is the average horizontal acceleration of the fluid. Substitution in 
Eq. (9) results in the following expression for the inertial force 


Fi = Cy (P) - Pg) Ax (12) 
From the equation of motion written for the vertical direction, tie follow- 
ing relationship can be obtained, 
p aa 3 ug 
in which w is the unit weight of the fluid. 


The substitution of the value of Ou,/ @t from Eq. (7) and iriegration 
yields 


p = -wz + emi cosh 27 (d + z)/L 


cos 8 +C (14) 
sinh 277 d/L 


The constant C can be evaluated because of the fact that p = 0 for points on 


the water surface, i.e., for z = y. The expression for the pressure then 
becomes 
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m7 LH cosh 27(d+z)/L - cosh 27(d+y)/L 
sinh 27° d/L 


p = w(-z+y)+0O cos (15) 


The value of (p; - Pg) to be used in Eq. (12) is found by inserting in Eq. (15) 
appropriate values of x and y for the two sides of a submerged body. 


For the particular case of a rectangular barge centered at the quarter 
point of a wave, substitution in Eq. (15) yields 


27x, 2 77H cos 277 x,/L 
= wH cos —;—|K - 2 sinh (16) 


in which K is defined as follows, 


cosh 27(d + z)/L 


cos a/L 


(17) 


and x, is the value of x at the upstream side of the model as shown in Fig. 5. 
Eq. (16) indicates that the pressure difference on the two sides of the body 
results from the difference in the water surface elevation, 4 H, 
(aH = H cos 2 7x /L) modified by the term in the brackets which corrects 
for vertical accelerations. 

The difference in the pressures occurring under a crest and a trough can 
also be obtained from Eq. (16) by letting x; = 0 and z=-d. Then 


(Pp. - P,) 

- K - 2 sinn? TH (18) 
in which p, is the pressure under a crest and Pt the pressure under an 
adjacent trough. Folsom(21,22) measured values of the ratio shown on the 
left side of Eq. (18) as part of a study of a pressure wave gage and compared 
them with K. He found that the measured values of the ratio were between 10 
and 20 per cent less than K. The correction provided by the additional term 
in the right side of Eq. (18) improves somewhat the correspondence between 
theory and the results of Folsom’s experiments. 


Wave Forces ona Flat Plate 


The results of the tests for horizontal forces on the flat plate are pre- 
sented in Table 1. The wave heights (column 4), measured forces (columns 
5 and 7), and location of model at the time of maximum force (columns 6 and 
8) were determined by averaging values from three or four successive waves. 
Heights of individual waves were taken as the average height of a crest above 
adjacent troughs. Separation of the upstream force from the downstream 
value was based on the assumption that position of the oscillograph pen cor- 
responding to zero force did not shift during a test. Although this change 
was small, it may have caused minor inconsistencies in individual values. 

The wave lengths shown in column 3 were computed by means of Eq. (1). 
The measurement of wave lengths was not part of the routine testing program. 
Nevertheless, a number of wave lengths were measured and the observed 
values corresponded with the computed values within 5 per cent. 
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In the analysis of the wave forces on a flat plate, both the drag force and 
the inertial force were found to be significant. The total horizontal force is 
determined from Eq. (10) by using the horizontal components of the velocities 
and accelerations. Values of F were computed by inserting into Eq. (10) the 
quantities which apply to this case. Values of Cg and C,, were selected to 
give the best correspondence with the measured phase angles and magnitudes 
of the maximum forces. A value of 3.5 was found for Cg and the value of 
CmV was found to be 1.75 times the volume of the circumscribing cylinder, 
which is 1.75 times the theoretical value(20) for an infinitely long thin flat 
plate. Because these values are large, an effect of the unsteadiness of the 
flow on the shape of the wake is indicated. Values of ux and Ou,/Ot were 
computed from Eqs. (4) and (6), respectively. 

For any particular value of wave height and wave period and for a par- 
ticular model location, F is a function only of the phase angle and Eq. (10) 
can be expressed as follows: 


=<,’ cos7e + sin @ (19) 
in which 
2 2 
cosh 277 (d+z)/L 
Cc," «GA (20) 
sinh 277 d/L | 
and 
277 24 cosh 27 (d+z)/ 
sinh 27° d/L (21) 


The phase angle corresponding to the maximum force is obtained by setting 
dF/d@ = 0, with the result shown in Eq. (22). 


sin@ (22) 
Computed values of 6,4, expressed in terms of x/L are shown in column 9 
of Table 1 and the corresponding values of Fg, Fj, and F are shown in 
columns 10, 11, and 12, respectively. Computed phase angles and forces are 
presented together with measured values in Table 1. These values are also 
compared graphically in Fig. 7, in which computed values are plotted against 
both upstream and downstream measured values. Although there is consider- 
able scatter in the data for the phase angle (Fig. 7a), the trends indicate good 
correspondence between measured and computed values. The correspondence 
between measured and computed forces, shown in Fig. Tb, is good for 
z/d = -0.25 and -0.50. For z/d = -0.75 the measured values exceed the 
computed values by about 30 per cent. 


Horizontal Forces on Barge Models 


Results of tests for barge models are presented in Tables 2 and 3. In 
Table 2 are shown results of a number of tests for which the data were 
analyzed in considerable detail. The values of wave height, maximum force, 
and phase angle at the time of the maximum force were determined by 
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x/L Computed 


x/L Meosured Measured Force in Pounds 


(9) (b) 


FLAT PLATE DATA 


Fig. 7 


averaging the results for three or four successive waves. The downstream 
forces were separated from the upstream values and are shown separately 
in Columns 8 and11l. The upstream forces were, on the average, 6.6 per 
cent smaller than the downstream forces. These values are subject to the 
same instrumental limitations as are those for the flat plate. 

The difference in water surface elevation on the two sides of the barge, 4 
oH, at the time of the occurrence of the maximum force were scaled from 
the oscillograph charts. The values are shown for the maximum downstream 
forces in column 9 and for the maximum upstream forces in column 12. 

In Table 3 are shown the results of a number of tests which illustrate the 
effect on the force of the position of the barge in the vertical. These tests 
were conducted on the rectangular model for three different wave heights 
with the wave period held constant at approximately 1.1 seconds. The location 
of the center of the barge (-z) was varied from 0.08 feet to 0.92 feet below 
the water surface. In the latter position, the barge was very close to, but 
not touching, the bottom. Three different wave heights were used as shown 
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TABLE 3 


HORIZONTAL FORCES ON RECTANGULAR BARGE 


For T = 1.1 Seconds oe 
Average Maximum Force in Pounds 


. H = .27 ft. H = .21 ft. H = .071 ft. 


-.08 2.99 2.85 1.07 
3.21 1.22 
-.25 2.87 94 
3. 2.54 .89 
-.42 2. 2.25 . 86 
-.50 2. 2.17 .83 
-.58 2. 2.05 75 
- .67 2. 2.03 .73 
«28 2. 2.09 75 
2. 2.10 .79 
> 2.25 .73 


in Table 3. Individual wave heights varied by less than five per cent from 
the average values for each group. The forces shown in Table 3 are the 
average of the upstream and downstream forces. From the results shown 
in Table 2, it can be estimated that the downstream forces were about 3 per 
cent larger than the average values shown. 

The maximum horizontal forces were created when the quarter points of 
a wave (x/L = 1/4 or 3/4) were approximately at the center of the barge as 
shown in Table 2, columns 7 and 10. Hence, the resistance caused by the 
velocity was negligible, and the maximum horizontal force was of the inertial 
type expressed by Eq. (9) or (12). Computed values of oH are shown in 
Table 2, column 14, together with the measured values for the rectangular 
and slotted barges (columns 9 and 12). The average of the computed values 
is 2 per cent higher than the average of the measured values for the down- 
stream forces and 11 per cent higher than the measured upstream values. 

Values of C,, were computed for each test from Eqs. (12), (16), and (17) 
and are shown in column 13. The value of x, is (L/4 - L'/2) in which L' is 
the width of the barge. For the case of the rounded barge an average width 
of 0.800 feet was used instead of the maximum width of 0.833 feet. The 
method of determining the values of Cy, is illustrated by the following com- 
putation for the test recorded at the top of Table 2. The value of K, from 
Eq. (17), is 0.784, and 4H = 0.101 (from column 14). The value of (p; - pg), 
from Eq. (16) is 4.94 psf. Finally, Cy, is found from Eq. (12). 
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Fy 3.07 
(P] - Pg) Ax 4.94 x 0.39 ~ 1-59 


For the slotted barge the displaced mass is used as a reference. In the 
computation of Cy, the forces on the inner faces of the slots (py" and pg‘) 
must be considered. The expression for the force is as follows: 


F = Ch - Pg) Ax - (py' - pa") Ax’| (23) 


in which A," is the area of the inner faces of the slots. 

The test data are summarized graphically in Fig. 8. The coefficient of 
inertial resistance for the rectangular barge is plotted against the ratio of 
wave height to wave length for three values of z/d with T constant at 1.1 
seconds in Fig. 8a. The values of C,, decrease significantly as the waves 
become steeper. The curves indicate the approximate trends. 

In Fig. 8b, values of C,, for the rectangular barge are plotted against the 
relative wave height for three wave periods, the barge position (z/d) being 
-0.50. In Fig. 8c, Cy, is plotted against relative wave height for three differ- 
ent barge shapes, for constant values of both z/d and T. 

The wave force is plotted against z for three different wave heights in 
Fig. 9. These are the data from Table 3 with all force values increased by 
three per cent to transform them from average to downstream forces in ac- 
cordance with the results presented in Table 2. The theoretical curves for 
the three wave heights were determined by means of Eqs. (12) and (16). The 
several constant values of C,, used in determining these curves were selected 
to give a good average fit. 

The tabulated data and the graphs indicate a general conformity between 
theory and test results. However, forces produced by the largest waves were 
less than would be expected from the results for the smaller ones. This is 
illustrated in Fig. 8 by the decrease in C,, with increasing wave steepness 
and also in Fig. 9 in which the value of C,,, which produces the best fit for 
the largest waves, is much less than that for smaller waves. Some deviation 
from the theoretical curves is also indicated in Fig. 9 for tests with the 
model very near the bottom or very near the water surface. 

The forces on the rounded barge were approximately 14 per cent less than 
those on the rectangular barge. The reasons for the lower force on the 
rounded barge can be explained with reference to Eq. (9), which shows the 
force to be proportional to the product of C,, and the displaced mass. Four 
per cent of the reduction is the result of the smaller displaced mass of the 
rounded barge. The remaining reduction in force was caused by the in- 
creased streamlining with the resulting lower values of Cy, as illustrated 
in Fig. 8c. 

Although the displaced mass of the slotted barge is 26 per cent less than 
that of the rectangular barge, the horizontal forces were only approximately 
19 per cent less. Thus, the reduction in mass is partially offset by the ad- 
ditional disturbance attributable to the odd configuration of the slotted barge 
and the resulting increase in C,, illustrated in Fig. 8c. 

Values of C,, for the rectangular model varied from 1.31 to 1.76. The 
value obtained from theoretical considerations of two-dimensional flow for 
this shape is 1.33 as shown in Fig. 6. Differences are expected because of 
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+ 2t/d*-25 

z/d*-50 + T#Q.87 Sec. 

1.8 e Tetl.tl Sec. 
x Ts 1.41 Sec. 


2-.75 


1.6 1.6 x 
Cm Cm 
1.4 1.4 
Rect. Barge ° Rect.Barge e 
T=1.1Sec. t/d= 


(b) 


+ Slotted Barge 
e Rectangular Borge 
* Rounded Borge 


COEFFICIENTS OF 


1.6 INERTIAL RESISTANCE 
For Horizontal Forces 
on Barge Models 
1.4 Fig. 8 
Tet.iSec. 
.75 
1.2 
02 04 06 


effects of periodic wake formation and because the effective dimension of 
the barge is an appreciable fraction of the wave length. 


Vertical Forces on Barge Models 


Test results for the determination of vertical forces are shown in Table 4. 
In this case, the dynamometers were not sensitive enough to warrant the 
separation of upward forces from downward forces. However, the records 
do indicate that upward and downward forces were nearly equal for inter- 
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Hs 0.07 Ft. 
Hs: 0.21 Ft. 
H= 0.27 Ft. 


Cm=* 1.47 


Rect.Barge 
T21.1 Sec. 


2 
Force in Pounds 


HORIZONTAL FORCES 


on Barge Models 
Fig. 9 


mediate values of z, and that the upward force was greater than the downward 
force if the barge was located near the water surface or near the bottom. 

The values of vertical forces, shown in Table 4, are the averages of the up- 
ward and downward forces. The location of the barge with respect to a wave 
crest at the time of occurrence of the maximum force is shown for a number 
of tests on each of the barges in columns 6 and 7. Wave lengths, shown in 
column 4, were computed from Eq. (1). 

The maximum downward vertical forces on the barge occurred if the crest 
was directly above the center of the barge (x/L = 0) as shown in column 6 of 
Table 4. The maximum upward force occurred in most cases if ihe trough 
was nearly over the center of the barge (x/L = 1/2), as shown in column 7. 
Thus, the vertical forces are also primarily inertial. Vertical forces can be 
analyzed by means of Eqs. (7) and (9). 

Values of C,,, obtained from Eqs. (7) and (9) and enumerated ia Table 4, 
can be compared with the classical values obtained from theory shown in 
Fig. 6.(18) Although theoretical results for rectangular bodies are available 
only if the length is infinite (two dimensional flow), experiments which have 
been conducted (19,20) indicate that the end effects for the lengths used are 
reductions in the virtual mass coefficient of five to ten per cent. For the 
proportions of the barge used in these tests, the theoretical value for two 
dimensional flow is 6.1, which is nearly the same as the average of the ob- 
served values. The average observed value for the ,ounded barge (referred 
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TABLE 4 


VERTICAL FORCES ON BARGE MODELS © 


1 


x/L Av .Max 
Type of T L H Downward | Upward Force 
Model Sec | Feet] Feet Force Force 


Pounds 
Rectangular] -0O. .09 -12]0.215] -0.10 0.26 


-0. .09 .12}0.198 .42 
.09 .12]}0.198 50 
.08 .06}0.196 .49 
-0. .13 .306 2? 
-0. .252 .39 
-0. . 066 .46 
-0. .210 41 
-0. 110 .45 
-0. .207 
-0. .187 
-0. .193 
-0. .199 
-0. .194 
-0. .205 
-0. .229 
-0. .230 
-0. 228 
-0. .147 
-0. .230 
-0. .135 
-0. 
-0. .148 
-0. .147 
-0. .152 
-0. 
-0. .322 
-0. .145 
-0. .151 
.289 
.316 
.317 
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to the somewhat smaller displaced mass of this shape) was 5.9, which is not 
significantly different from that of the rectangular barge. In comparison, 
for an elliptical section with the same proportions, an extreme in rounding, 
a theoretical value of 6.33 is obtained for two dimensional flow. Finally, an 
average value of Cy, of 3.8 was obtained from the tests on the slotted barge 
using its displaced mass as a reference. Although no theoretical value can 
be obtained for this shape, the weighted average of the theoretical values for 
two dimensional flow, for the unslotted portion (6.1), and for the legs (2.5), 
is 4.9. In view of geometrical differences and possible effects of separation, 
the theoretical and observed values of C,, correspond surprisingly well. 

The data are also shown graphically in Fig. 10. The forces on the rec- 
tangular and rounded models were plotted for various values of z in Fig. 10a. 
The group for the rectangular barge for a wave period of 1.1 seconds com- 
prises the values shown in Table 4 for which the wave heights were approxi- 
mately 0.20 feet. The group for the rounded barge is made up of the values 
from Table 4 having wave heights of approximately 0.23 feet. For this group, 
the forces were reduced by the ratio 0.20/0.23 before plotting, because, for 
a small range in wave height, the forces vary linearly with the wave height 
in accordance with the theory. The curves shown in Fig. 10a were deter- 
mined from Eq. (7) and (9), using a Cy, of 5.5. This value was selected to 
give the best fit for this group of data. In Fig. 10b are shown groups of 
points for the rounded and slotted models. These values are the ones shown 
in Table 4 for wave heights of approximately 0.15 feet. The curves in Fig. 
10b were determined from Eqs. (7) and (9) using the values of C,, shown. 

Good agreement was found between analytical and test values of the verti- 
cal forces on the three models used in this study. Values of C,, correspond 
quite well with the theoretical values for two dimensional flow except for the 
cases in which the models were very near the surface or the bottom. In 
these cases, there was some indication that the upward forces were larger 
than the downward ones so that the values reported here (which are the 
average of the upward and downward forces) should be increased somewhat 
for design purposes. Further studies are required to determine the exact 
nature of the forces in these locations. 


SUMMARY 


Forces exerted on submerged bodies by waves can be reliably predicted 
as a combination of the inertial and drag forces if the appropriate coeffi- 
cients can be determined. The tests on the models of barge-like structures 
showed that the inertial forces were predominant so that only the coeffi- 
cient of inertial resistance was needed to predict the maximum force. For 
the horizontal forces, the measured coefficient was approximately 1.5 
compared to 1.33 as estimated from theory for flow without separation. For 
the vertical forces, the measured and estimated values of the coefficient 
were approximately six for the unslotted barge. Some systematic variation 
in the coefficient was found with relative wave height and with the vertical 
position of the model. Measured values for the slotted and rounded barges 


were primarily inertial, the rounding of the edges of the barge did not cause 
large changes in the forces. 


The maximum horizontal forces on a vertical plate placed normal to the 


were also in good correspondence with theoretical values. Because the forces 
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—~+— Rectangular Model 
Rounded Mode! 
Model 


4 
/ 
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“7 
Te -8 T= 1.1 Sec. 
He 0.2 Ft. H= 0.15 Ft. 
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VERTICAL FORCES 


on Barge Models 


Fig. 10 


direction of wave propagation were predictable if a drag coefficient of 3.5 
and a virtual mass equal to 1.75 times that of the circumscribing cylinder 
were used. Both of these values are large because the growth of the wake in 
unsteady flow increases the drag coefficient and because the presence of any fs 
wake increases the virtual mass. The phases of the wave which produced the ae 
maximum forces corresponded satisfactorily with the predicted values. 
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SNOWMELT RUNOFF 


J. Harold Zoller, ! M. ASCE and Arno T. Lenz,2 M. ASCE 
(Proc. Paper 1834) 


SYNOPSIS 


The primary factors relating to the melting of snow were evaluated for the 
snowmelt periods of 1938 through 1952 for the Big Eau Pleine River basin in 
Wisconsin. Convection and condensation melt potentials were computed from 
equations developed by previous investigators. A nomograph was devised for 
computing radiation snowmelt potential from insolation, as related to sun po- 
sition determined by date, modified by cloud cover and albedo. The as- 
sumption was made that the runoff for any day would be a function of the sum- 
mation of melt potentials for that day. A snowmelt unit hydrograph was de- 
veloped and relative daily hydrographs were computed for each day. These 
relative hydrographs were then combined by a trial and error procedure to 
reproduce the recorded snowmelt hydrograph for each year. From this 
process the daily runoff values which contributed to the total hydrograph were 
computed. These daily runoff values were compared with the total melt po- 
tentials for the same days, and the correlation noted was expressed by an 
equation. 


INTRODUCTION 


Spring runoff caused by the melting of the winter’s accumulation of snow 
often provides an opportunity for the annual storage of a major part of the 
regional water resources. Reservoirs may be drawn down to produce water 
power during the high demand winter months with certain knowledge that melt- 
ing of the existing snowpack will replenish the reservoir supply in the spring. 


Note: Discussion open until April 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1834 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 6, November, 1958. 

1. Prof. of Civ. Eng., Univ. of New Hampshire, Durham, N. H. 

2. Prof. of Civ. Eng., Univ. of Wisconsin, Madison, Wisc. 
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Research has revealed techniques by means of which the quantity of water 
which will become available in the snowmelt runoff of any basin may be esti- 
mated with considerable accuracy if certain minimum data are acquired. 

The pattern of key research in the past has been to concentrate in small 
control areas the instrumentation necessary to measure every conceivable 
quantity suspected of being an agent in the melt phenomenon. This approach 
has yielded, and is continuing to yield, information of a most vital nature. 
However, to be of any practical value, this information must also be applicable 
to large basins where average instrumentation exists for the recording of 
pertinent data. If such applicability can be universally demonstrated, then the 
ultimate goal of the researchers will be attained. 

Most of the concentrated research in the field of snowmelt has been con- 
fined to the mountainous areas of western United States, where much of the 
annual precipitation occurs in the form of snow. Here, snow may accumulate 
to great depths in the winter months and its melting season may extend from 
March to July or even later at the higher elevations. In these areas melting 
is almost always confined to the daylight hours and occurs chiefly as a result 
of insolation, air temperature and air turbulence produced by wind and to- 
pography. 

In contrast, the snowpack in central Wisconsin is seldom over two feet in 
thickness at the time of initial melting and is commonly nearer one foot in 
thickness. It lies largely in open fields on gently rolling topography. Its melt- 
ing season which also begins in March, never extends beyond April and has a 
continuous duration of from four to twenty days. During periods of high melt 
potential, melting may occur throughout the night as well as during the day- 
light hours, and the melting achieved by condensation of moisture on the cold 
snow surface may also be a major factor. Thus the problem of forecasting 
the rate of runoff from melting snow may be recognized as more difficult in 
many ways under the conditions imposed by Wisconsin’s climate than would be 
the case for many regions farther west. 


Melt Potential Theory 


The primary factors relating to the melting of snow are: Temperature of 
the air, temperature of the soil, wind velocity and air turbulence, humidity, 
rainfall and cloud cover. These factors are measurable, within certain limits, 
and from their recorded values the potential melting capacities attributable to 
convection, condensation, radiation, conduction and rainfall may be deduced. 

In addition, certain basic characteristics such as elevation, aspect, topography, 
cover, latitude and soil conditions are influential in the final melting capacity 
of a given melt potential. 

A study was made of the melt potential—runoff characteristics of the Big 
Eau Pleine River Basin of Central Wisconsin. In this study the assumption 
was made that runoff in the snowmelt season was a direct function of the heat 
exchange relations between the snowpack and its environment. Accordingly, 
the total melt potential was determined for every six-hour interval in the 
snowmelt periods of the years 1938 to 1952. This total melt potential for any 
interval was considered to be the algebraic sum of the radiation melt po- 
tential, the convection melt potential and the condensation melt potential for 
that interval. It was further assumed that the runoff which occurred for any 
day would occur in a distribution indicated by the variation in melt potential 
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for that day. For example, if the total melt potential on a given day for the 
six-hour interval beginning at 9 a.m. were 1.5 inches and the total melt po- 
tential for the interval beginning at 3 p.m. were 0.5 inches, all other melt po- 
tentials for the day being negative, then it was assumed that the runoff hydro- 
graph for that day’s melt was a combination of the 9 a.m. and the 3 p.m. hydro- 
graphs in the ratio of 3 to 1 respectively. This procedure involved the de- 
termination of a snowmelt unit hydrograph and the determination of average 
daily discharges due to periods of melting beginning at different hours in the 
day. Negative melt potential summations were ignored because their effect is 
minimized by the low conductivity of the snowpack and because the time du- 
ration of negative melt potentials was short. Finally, it was reasoned that, if 
the theory and assumptions listed above were correct, a correlation would be 


found between the daily melt potential and the daily runoff for the various 
years being investigated. 


Convection Melt Potential 


The convection melt potential is the quantity of water in inches which theo- 
retically would be melted from snow of 100 per cent quality in a given time 
interval by the convective transfer of heat. The factors which are of im- 
portance in determining convective heat transfer are air temperature and 
turbulence, as indicated by wind velocity. Because preliminary studies with 
intervals of time equal to one day proved ineffective, the unit of time finally 
selected for investigation was six hours. Air temperatures and wind velocities 
were recorded at the Wausau CAA airport at 12:30 a.m., 6:30 a.m., 12:30 
p.m. and 6:30 p.m. for the period of study. The 9 a.m. to 3 p.m. six-hour 
period is generally the most important snowmelt period of the day and the 
12:30 p.m. temperature and wind velocity data may be considered fairly repre- 
sentative of that period. Thus, the four quarter-day periods were selected to 
begin at 3 a.m., 9a.m., 3 p.m., and 9 p.m. It was assumed in all cases that 
the wind velocity and temperature readings were average values for the six- 
hour periods in which they occurred. 

Although the Wausau CAA Airport is located completely outside the upper 
Big Eau Pleine Basin, approximately 30 miles east of the center of the basin, 
the Wausau data are considered applicable. Temperatures at Wausau and 
Medford, 20 miles northwest of the center of the basin, exhibit a very good 
correlation for the snowmelt period. 


The convection melt potential was determined from Wilson’s (1941) 
equation: 


M = KV(T-32°) 


“where M is the depth of water in inches melted from snow in six hours, V is 
the wind velocity in mph, T is the dry-bulb temperature in degrees F and K 
is a constant involving the latent heat of ice, exposure of instruments, air 
density, conversion of units, and certain considerations involved in the theory 
of turbulence.” A value of K = 0.002 was used. This value of K corresponds 
closely to Light’s (1941) value of K = 0.00176 as computed for the physical 
conditions imposed by the Big Eau Pleine Basin. The convection melt po- 
tential data for the snowmelt period of the year 1943 are tabulated in Table I. 


Condensation Melt Potential 


The condensation melt potential may be defined as the amount of water in 
inches liberated from snow, plus condensate, as condensation occurs on the 
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TABLE I 
CONVECTION AND CONDENSATION MELT POTENTIALS - 193 


Date Period Dry Bulb Dew Wind Convection Melt Condensation 
Temperature Point Velocity Potential in Melt Potential 


in °F in °F in mph Inches in Inches ‘ 
March 
23 1 18 15 2 206 -.02 
2 39 26 20 02 
3 39 29 16 ece 7 
33 29 8 202 
2h 1 33 30 10 202 -.92 
2 46 32 6 016 ) q 
43 32 8 218 18) 
3h 32 2 
25 1 30 30 8 
2 ko 36 295 
3 43 37 230 
39 38 10 ol 205 
26 1 34 2h 20 208 012 
2 36 19 20 018 -18 
3 3h 17 18 208 -38 
27 16 u 6 -.07 4 
2 26 13 12 012 
28 13 8 
23 14 6 
28 19 16 12 
2 29 25 12 
3 30 23 1h 
29 25 15 10 
29 1 28 2k 
2 38 29 18 e22 -.05 
3 37 7 202 
39 16 036 010 
30 1 L6 6 016 207 
2 61 47 12 «70 
3 61 10 058 
53 1h 258 
31 1 us 2 203 
2 47 37 22 256 299 
3 us 35 12 230 203 
k 38 26 12 e -.06 
April 
1 1 33 30 6 202 -.91 
33 31 7 292 
3 36 22 10 208 207 
28 19 12 
2 1 23 19 
2 31 16 18 -.18 
3 30 uy 10 
22 10 2 =.02 
3 1 18 pin 4 
2 35 18 16 10 
3 37 25 18 018 
h 36 23 12 010 -.08 
1 32 26 8 
2 36 28 298 209 
3 48 26 28 290 
39 2h 20 228 


* Periods 1 to i at 3 asme, 9 aeme, 3 Pome, and 9 pom respectively. 
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surface of the snow. The factors which determine the condensation melt po- 


tential are the dew point and the wind velocity, as related by Wilson’s (1941) 
equation: 


M = a V (e-6,11) 


where M is the depth of water melted plus condensate per six hours in inches, 
a is a constant involving the latent heat of ice, exposure of instruments and 
conversion units, and is equal to 0.0032 for basins with elevations less than 
3000 feet, V is wind velocity in miles per hour, e is the vapor pressure of the 
air in millibars, and 6.11 is the saturation vapor pressure in millibars at 

32° F. Records of dew point and wind velocity at Wausau were again con- 
sidered applicable to the Big Eau Pleine Basin. The condensation melt po- 
tential data for the snowmelt period of 1943 are tabulated in Table I. 


Radiation Melt Potential 


Adequate data for the direct determination of incident solar radiation for 
the Big Eau Pleine Basin were not available. Thus it became necessary to 
reconstruct the insolation characteristics from the six-hourly cloud cover 
data, recorded at the Wausau CAA Airport. Since records of both cloud cover 
and incident solar radiation were available for Madison, Wisconsin, a corre- 
lation between these factors was established for use with the Wausau data. 

Cloud cover is an estimated quantity, reported in tenths by a ground ob- 
server. A cloud cover of zero represents a cloudless sky, 10 represents a 
complete overcast and intermediate values represent intermediate conditions. 
It must be emphasized however, that although two days may be reported 
correctly as having 10-tenths cloud cover, the insolation reaching the earth 
for those two days may be vastly different since cloud canopies differ in their 
capacity to filter out solar radiation. A further indication of the inaccuracy 
involved in the use of cloud cover data in determining incident radiant solar 
energy is illustrated by several periods in which the noon cloud cover report- 
ed was 10-tenths and the duration of sunshine reported for the hours preceding 
and following the noon cloud cover observation were each 60 minutes. In other 
words, although the sky was completely overcast, the clouds were of sucha 
nature that the sunshine was visibly able to penetrate them continuously. 
Obviously, the insolation corresponding to 10-tenths cloud cover on such a day 
will not correlate properly with the insolation corresponding to 10-tenths 
cloud cover on a day when heavy clouds are extant. 

The Madison insolation records for several 15-day periods were checked 
for the 19-year period 1931-1949 inclusive, and the highest insolation value 
for each daylight hour within each interval was determined. From these 
values were determined the maximum composite insolations which are equal 
to the daily insolation in langleys which would have occurred at Madison if the 
highest insolation values for each hour within a given period had all occurred 
on the same day. These values were plotted against time to give the Maximum 
Composite curve shown on Fig. 1. This curve represents the daily insolation 
summation which is approximately the highest insolation that could possibly 
occur at Madison for each day represented. The average insolation experi- 
enced in Madison is seen from Fig. 1 to be a function of this maximum com- 
posite insolation. 

The latitude correction to adapt the Madison data to the Big Eau Pleine 
Basin is 0.97. This correction was considered insignificant in comparison 
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with uncertainties involved in other assumptions and was not used. Its neglect 
results in a constant error for which compensation may be made in the se- 
lection of the proper albedo value. 

Insolation values from Madison were converted to percentages of the maxi- 
mum composite insolation by dividing the insolation value recorded by the 
maximum composite insolation for that day as determined from Fig. 1. These 
data were then plotted against various combinations of cloud cover data. It 
was found that a significant correlation existed between the per cent of maxi- 
mum composite insolation per half day and the sum of either the 6:30 a.m. 
and 12:30 p.m. cloud cover values or the 12:30 p.m. and the 6:30 p.m. cloud 
cover values, expressed in tenths. Table Il summarizes the results of this 
correlation study. It is significant to note that errors resulting from approxi- 
mating the insolation from cloud cover data are later minimized materially 
by corrections for albedo and back radiation. It is also emphasized that the 
greatest accuracy of insolation determination occurs during periods of mini- 
mum cloud cover, or during periods when maximum melting from insolation 
would occur. 

Thus, from the cloud cover data at Wausau for any day the per cent of the 
maximum composite insolation for that day may be determined from Table IL. 
For example, on March 30, 1943, the 6:30 a.m., 12:30 p.m. and 6:30 p.m. 
cloud covers were 0.4, 0.5 and 0.8, respectively. The sum of the 6:30 a.m. 
and the 12:30 p.m. cloud cover values is 0.9 and the sum of the 12:30 p.m. and 


TABLE II 


INSOLATION-CLOUD COVER RELATIONS 


Summation of Percent of Summation of Percent of 
6:30a & 12:30p, Maximum 6:30a & 12:30p Maximum 
or 12:30p & 6:30p Composite or 12:30p & 6:30p Composite 
Cloud Cover-Tenths Insolation Cloud Cover-Tenths Insolation 
per half day per half day 

AR 10 3h 

1 43 1 33 

2 43 12 31 

3 k2 13 29 

27 

5 kl 15 25 

6 ko 16 23 

7 38 17 21 

8 37 18 19 

9 36 19 17 
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the 6:30 p.m. cloud cover values is 1.3. From Table II the morning insolation 
was 36 per cent of maximum for the half day and the afternoon insolation was 
29 per cent of maximum, so the daily total was 65 per cent of the maximum 
composite insolation for that day. From Fig. 1 the maximum composite inso- 
lation (March 30) was 658 langleys. Thus, the computed insolation for March 
30, 1943 was 428 langleys. 

The radiation components effective in producing snowmelt potential may be 
considered to follow the pattern shown in Fig. 2. The shortwave incident radi- 
ation (428 langleys on March 30) consists of the direct and diffuse solar radi- 
ation. The net longwave emissive radiation consists of back-radiation from 
snow, partially balanced by incident diffuse longwave sky radiation. This net 
longwave radiation has been shown to be emissive (Fig. 2b) and essentially 
constant at 0.10 langleys per minute. During hours of darkness when there is 
no insolation the energy loss has been termed non-concurrent radiation, Dur- 
ing daylight hours this emissive (concurrent) radiation must be deducted from 
the incident radiation to obtain the net radiation. It is evident that the duration 
and amount of the concurrent and non-concurrent radiations are functions of 
the total daily insolation. Thus, the insolation effective in producing melt for 
any day is the computed insolation for that day minus the concurrent back 
radiation. 

A variable amount of the incident solar radiation is reflected from the 
snow surface and is thus unavailable for melting. The percentage reflected, 
called albedo, to be used in reducing the effective insolation is a function of 
the density, structure, moisture content, depth and purity of the snowpack. 
Because sufficient records were not available to permit a proper evaluation 
of albedo variation, a constant albedo value of 60 per cent was assumed for 
the entire snowmelt period. For an albedo of 60 per cent, the computed inso- 
lation for March 30, 1943, of 428 langleys would become 171 langleys. 

The net concurrent radiation is a function of the envelope of effective radi- 
ation and has been determined from the authors’ emperical equation: 


r = 38.5 + 0.0633 p = 49 langleys 


in which: r is the net concurrent radiation in langleys per day and p is the 
net insolation as reduced for albedo in langleys per day. Subtracting the net 
concurrent radiation of 49 langleys would leave 122 langleys, which would be 
sufficient to produce 122—~(80 x 2.54) = 0.60 inches of melt. During the same 
day, the net non-concurrent radiation was (60 x 24 x .1) - 49 = 95 langleys, 
which is the amount of heat liberated in freezing 0.47 inches of water at 32° F. 
The insolation nomograph of Fig. 3 provides an easy means of solving for 
the radiation melt potential. To solve for the melt potential using Fig. 3 it is 
necessary to know the date, the per cent of total insolation and the albedo. A 
straight edge is placed on the nomograph from the date to the per cent of total 
insolation and the intersection with the net insolation line is noted. A hori- 
zontal line through this net insolation will determine the melt for various 
albedo values and the amount of water frozen by non-concurrent radiation. 
Utilizing the data of March 30, 1943, a straight edge placed from March 30 to 
65 per cent indicates a net insolation of 428 langleys per day. A horizontal 
line through this net insolation value intersects the 60 per cent albedo line at 
a melt of 0.60 inches per day and indicates that 0.47 inches of water would be 
frozen by non-concurrent emissive longwave radiation. The radiation melt 


potential calculations for the snowmelt period of the year 1943 are tabulated 
in Table Il. 
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The radiation melt potential values thus determined could not be utilized 
directly in the snowmelt calculations however, because they had to be assimi- 
lated with other melt potential factors which were determined for one-quarter 
day intervals. From a study of daily insolation graphs it was decided arbi- 
trarily to consider that two-thirds of the daily positive radiation melt potential 
occurred in the interval from 9 a.m. to 3 p.m. and one-sixth occurred in each 
of the intervals 3 a.m. to 9 a.m. and 3 p.m. to 9 p.m. It was further concluded 
that the negative melt potential should be equivalent to 36 langleys or 0.18 
inches of melt for the period from 9 p.m. to 3 a.m. and that the remainder of 
the negative melt potential should be equally divided between the 3 a.m. to 
9 a.m. period and the 3 p.m. to 9 p.m. period. Using this scheme of division, 
the melt potential for March 30, 1943 would be subdivided as indicated in 
Table IV, according to which two-thirds of the .60 inches of melt potential, or 
0.40 inches, is assigned to period 2 and the remaining 0.20 inches is equally 
divided between periods 1 and 3. A total negative melt potential of 0.47 inches 
is divided so that 0.18 inches is applied to period 4 and the remainder of 0.29 
inches is equally divided between period 1 and period 3. 


Total Melt Potential 


The total melt potential for any time interval is merely the algebraic sum 
of the individual melt potential factors. Table V gives the total melt potentials 
as computed for the snowmelt season of the year 1943. 


Snowmelt Runoff Theory 


The streamflow hydrograph resulting from melting snow must reflect the 
variation in total melt potential. The hydrograph is also a function of a snow- 
melt unit hydrograph because, like runoff from rainfall, once an increment of 
flow is started on its path it proceeds until it has passed the measuring 
station, unless extremely severe freezing conditions follow melting. 

By preparing a “relative hydrograph” for each day having a positive melt 
potential and making these relative hydrographs reflect the daily variation in 
melt potential as well as the snowmelt unit hydrograph pattern it is possible 
to divide the daily runoff into the increments which have been contributed 
from snowmelt potential factors existing on previous days. 


To illustrate this procedure in detail numerous examples have been used 
in the following discussion. 


TABLE IV 
RADIATION MELT POTENTIAL FOR MARCH 30, 1913 


Period Positive Melt Negative Melt Net Melt 

Potential Potential Potential 

(1) 3am to 9am 010 ine in. ins 
(2) 9am to 3pm 40 
(3) 3pm to 9pm 210 
(4) 9pm to 3am -.18 
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Date Period 
March 
23 1 
2 
3 
2h | 
2 
3 
25 1 
2 
26 i 
2 
3 
27 1 
2 
3 
28 1 
2 
3 
29 1 
2 
3 
4 
30 
2 
3 
4 
31 1 
2 
3 
April 
1 1 
2 
3 
4 
2 1 
2 
3 
4 
3 1 
2 
2 
4 1 
2 
3 
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TABLE V 


TOTAL MELT POTENTIAL = 19)3 


Melt Potential in Hundredths of an Inch - 
Convection Condensation Radiation Positive 


Total 


51 
15 


26 


Rain in 
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Inches 


293 


298 


bis 
+4 -2 
28 -10 33 
2 -2 -18 
16 0 Lh 60 x 
18 0 -3 15 | 
-1 0 
22 5 53 80 
30 6 0 36 | 
8 -12 | 
18 -18 53 53 
=22 -16 -18 
=32 -12 = 
-10 -7 -18 
“12 -9 -2 g 
18 3 5 
36 10 -18 28 
58 18 72 
58 20 -18 60 
a 66 9 19 94 : 
30 3 22 
2 -1 43 bb 
-26 0 = 
-18 53 31 
-4 -18 : 
10 -15 39 3 
18 -10 -4 4 
10 8 -18 
: 98 3 52 159 
28 -12 -18 2 
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Snowmelt Unit Hydrograph 
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A trial snowmelt unit hydrograph was determined by isolating the runoff 
from one day’s melt in a manner similar to that used in determining a rainfall 
unit hydrograph. Base flow was removed from consideration by utilizing the 
fall recession curve of Fuhriman (1952). In addition, the effect of the duration ‘ 
of the melt potential producing runoff was considered. The resulting snowmelt 
unit hydrograph, shown on Fig. 4a, is in close agreement with the rainfall unit 
hydrograph (Beck, 1946), leading to the conclusion that, at least in certain in- 
stances, the same unit hydrograph may be used for rainfall and snowmelt run- . 
off. The average daily discharge in cfs for melting beginning at four different 
times during the day are given in Table VI. Fig. 4b compares the recorded 
hydrograph with the reconstructed hydrograph for the snowmelt period of 1939 
by means of the procedures outlined herein. 


Ice Index Method of Altering Published Data 


According to the United States Geological Survey Water-Supply Paper No. 
1085 (1947) the runoff records for the Big Eau Pleine River near Stratford, 
Wisconsin are “. . . good except those for periods of ice effect and those below 
10 second-feet, which are fair.” The notation “good” “. . . indicates that, in 
general, the error in the daily records is believed to be . . . less than 10 per 
cent; ‘fair’ less than 15 percent;. . .” 

Several times during each year actual discharge measurements are made 
at Stratford and compared with the apparent discharge as determined from the 
concurrent gage height and the rating table for the station. Frequently a head 
correction of a few hundredths of a foot is determined from these actual 


TABLE VI 
AVERAGE DAILY DISCHARGE FOR 
SNOWMELT UNIT HYDROGRAPH 


Discharge in cfs for Melt Beginning at: 


3 am 9 am 3 pm 9 pm 
(Period 1) (Period 2) (Period 3) (Period k) 


5 


1 2753 1603 326 5 
2 2218 3068 3937 3667 
3 667 868 1137 1547 
4 218 282 377 507 
5 1,0 168 
6 46 59 68 80 | 
7 18 22 28 38 
8 1 4 7 12 


Total 6018 cfs days 6025 6020 602); 


Ratio of 2753/326 1603/326 326/326 5/326 
"First Day" or or or or 
Discharges 1.00 0.0153 
(Based on 

Period 
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discharge measurements which correction is then applied to the computed dis- 
charges for a period of time as determined by earlier or subsequent discharge 
measurements. It is quite natural that these corrections should be required 
as a stream bed tends to erode and redeposit its sediment in a slightly differ- 
ent fashion following each passage of high water. The variations however are 
slight and no alteration of the rating curve has been necessary since 1948, at 
which time only minor changes in the discharges were made for gage heights 
less than 4.5 feet. 

During the winter season ice forms on the surface of the stream and the 
discharge characteristics change from those of open channel flow to those of 
enclosed flow. During periods of prolonged cold weather the apparent dis- 
charge as computed from the recorded gage heights and the rating table will 
tend to increase markedly, while true discharge as determined by field 
measurements is actually decreasing. This is explained by the fact that the 
resistance to flow is markedly increased by the surface contact with the ice, 
while at the same time, the weight of the ice supported by the water causes an 
increase in the water level measured in the float well. As an example, an 
apparent discharge of 110 cfs may be determined from the gage reading at 
Stratford, while an actual measurement by current meter may indicate that 
the discharge is 10 cfs. 

Usually two or three discharge measurements are made each year during 
the period when ice effects are noted. The critical perjod during early snow- 
melt runoff, however, is not well represented by actual measurements because 
rising water makes the ice unsafe for measurement operations through the 
ice while complete cover exists, and floating ice renders the taking of velocity 
measurements impractical after the break-up begins. Correct discharge 
values for the early snowmelt runoff season are the most difficult to de- 
termine and are the ones most likely to be in error in the published data. The 
USGS engineers determine the discharge values reported in the Water-Supply 
Papers from the gage height readings, concurrent air temperatures, reservoir 
storage and discharge records, and hydrographs of previous runoff. The 
method of accounting for ice effects presently being used on the Big Eau 
Pleine River by USGS engineers consists of plotting the apparent discharge, 
as determined by recording-gage readings, on semi-log paper, along with the 
actual discharges determined by the winter discharge measurements which 
are made about two times each winter. A line is then drawn through the 
measured discharges approximately parallel to the apparent discharge and the 
discharge values to be recorded are read from this line. Fig. 5 shows the 
manner in which winter discharge values were determined by this process for 
the winter 1950-1951. Fig. 6 shows the same discharges plotted on an 
arithmetic plot for the break-up period of 1951. It will be noted in Fig. 6 that 
at the higher discharges, the use of the USGS method results in extremely 
large reductions in discharge for ice effects. It has been deemed advisable, 
therefore, to alter certain critical discharge values from those given in the 
Water-Supply Papers for the purpose of this investigation in accordance with 
the ice index method outlined below. 

All winter discharge measurements on the Big Eau Pleine River were ana- 
lyzed. The gage reading which would have been required to produce the actual 
discharge as determined by measurement (without ice effect) was obtained 
from the rating table. This reading was subtracted from the measured gage 
reading to give the “ice index”. It is reasoned that the ice index is equal to 
the equivalent water depth of the ice cover plus an increment of thickness to 
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compensate for the reduction of flow caused by the drag of the water on the 
ice cover. Values of ice index were computed for all days when winter dis- 
charge was measured, and are plotted against the date of the measurement in 
Fig. 7. The maximum ice index value determined by actual measurement is 
1.65 feet. 

It is reasoned that the ice indices of Fig. 7 may be applied to the recorded 
gage heights to determine the probable discharge for those days when only the 
gage readings are available. For example, on March 17, 1945, the average 
gage reading of 12.95 feet was recorded. This is equivalent to a discharge of 
5890 cfs without ice correction as determined from the rating table. Ifa 
maximum ice index of 1.60 feet were considered to apply, the corrected gage 
reading of 12.95 minus 1.6 or 11.35 feet would provide a flow of 4315 cfs for 
this date. The discharge reported for this date in USGS Water-Supply Paper 
No. 1035 is 1500 cfs, which would occur normally with a gage reading of 7.64 
feet as determined from the rating table. Thus the computed ice index for 
this date would be 5.31 feet. Such an index could be correct if an ice jam oc- 
curred just below the gaging station and backed up water at the station; or if 
an ice canopy about 5.31 feet in thickness were formed over the stream at the 
gaging station for the full channel width. The absence of any obstruction to 
flow, at the gaging station or downstream, which might cause large ice dams 
to form leads to the conclusion that such reductions in flow as a result of ice 
effect are undesirable for the purposes of this investigation. 

When published discharge data were examined from all the years of record 
for which good gage readings are available in the initial snowmelt period, it 
was found that certain inconsistencies exist, relative to a study by the ice 
index method. During the days immediately preceding the occurrence of snow- 
melt runoff, the ice index had a reasonable value between 0.7 and 1.7 feet for 
all years. During the early snowmelt period the maximum ice effect based on 


7 
= 
a 
4 
. 
7 


SNOWMELT RUNOFF 


ASCE 


AMUYNNYS 


QON3937 


INIZD Td of 


3Wil — X3Z0NI 3D! 


X3GNI 


1334 Ni 


1834-19 
Q : 
4 
a \ 
4 \ \ 
7 
tb a 
\ 
> 
a 
\ 
N \ \ 4 
N 
0 
\ 
» 
a 
| \ 
\ 
. 
> 
\ 


1834-20 


HY 6 November, 1958 


published discharge data, varied from year to year, at times reading maxi- 
mum values as high as 7.10 feet as in 1938, and as low as 0.74 feet in 1946. 
Furthermore, in some years the ice index was considered to become zero on 
the rising side of the hydrograph, at a discharge of 109 cfs as in 1946, or 
8950 cfs as in 1945; and in other years it became zero on the falling side of 
the hydrograph, at a discharge of 2240 cfs following a peak apparent discharge 
of 8950 cfs as in 1938. In addition, the ice index for some years was erratic, 
alternating between low and high values on successive days. 

To eliminate these variations, the following rules for reconstituting the 
average daily discharge for days in the early snowmelt period were adopted. . 


1. In the hydrograph summation computations the ice index correction at 
the beginning of snowmelt was assumed to be about equal to that de- 
termined from a previous measurement and was never more than 1.65 
feet. 

2. The ice index was assumed to approach zero on the rising limb of the 
hydrograph at a gage discharge of less than 5000 cfs. 

3. The exact discharge for the last day of ice effect was determined from 


the hydrograph summation analysis of the total seasonal flow, as out- 
lined later. 


Determination of Daily Melt 


The analytical determination of the amount of snowmelt required each day 
to reproduce the recorded snowmelt hydrograph is a trial and error pro- 
cedure. The complete analysis for the year 1943 is presented to illustrate the 
technique developed. 

Basically, the method consists of applying a runoff for any given day which 
will make the total runoff for that day exactly equal to the average daily dis- 
charge as determined from the USGS Water-Supply Papers, altered, as neces- 
sary, for assumed ice conditions. The runoff so applied for any day is then 
followed by a sequence of discharges for succeeding days as determined from 
what has been termed the “relative hydrograph”. Thus, the magnitude of first 
day’s runoff, resulting from a given period of melting, is of primary im- 


portance because it automatically determines the magnitude of the discharges 
for succeeding days. 


Relative Hydrograph 


The relative hydrograph for any day may be defined as the succession of 
average daily discharges resulting from melt occurring in the proportions 
given by the variation in melt potential for that day, and so related that the 
runoff corresponding to period 2 is equal to one inch. Since the melt potential 
variation for each day is usually different from that of every other day, it 
follows that each day will have a different relative hydrograph depending upon 
its melt potential variation. Table VII lists the melt potential data for 1943 
which is necessary for the computation of relative hydrographs, and a summa- 
ry of the final runoff data computed by the method of hydrograph summation, 
discussed later. The relative hydrograph computations are tabulated in Table 
VII and are determined in the following manner: 

From Table VII, beginning with March 25, 1943, it is noted that the melt 
potentials for periods 2, 3 and 4 are 0.80, 0.36 and 0.01 inches respectively. 
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TABLE VII 
MELT POTENTIAL AND RUNOFF - 19,3 
Melt Potential in Inches 
Date Period Total Cuwmilative Tota 
(1) (2) (3) (4) Total 
March 
23) 263 063 
) Ripening Period 
1.38 
25 280 036 1.17 2055 202 
26 053 053 3-08 208 
27 ell ell 3219 e22 
28 026 026 3045 220 
29 055 095 228 098 1.33 
30 019 1.28 ote 260 2079 
31 222 1.16 8.38 258 
April 
1 of oth 8.82 202 
2 e31 e31 9.13 015 
3 038 9-51 207 
4 1.59 076 2035 11.86 226 
TABLE VIII 


RELATIVE HYDROGRAPHS - 19143 


Relative Average Dail pescrartge - cfs 
Date Period First Second Third Fourth Fift n ies nth Eighth 
Day Day Day Day Day .~ Day Day 
March 
25 2 1603 3068 868 282 119 59 22 h 
512 170 63 31 13 3 
19 6 2 1 
Total 1750 4886 1399 58 18) 91 35 7 
29 2 1603 3068 868 282 119 59 22 4 
3 25 303 87 29 n 2 1 
4 2 1579 666 218 72 3, 16 5 
Total 1630 950 1621 529 202 98 Me) 10 
30 1 08 329 99 32 uw 7 a 
2 1603 3068 868 282 119 59 22 4 
3 183 = 2215 640 212 79 38 16 4 
A 2 1719 725 238 79 38 18 6 
Total 2196 7331 2332 764 291 142 59 wy, 
31 2 1603 3068 868 282 119 59 22 A 
3 76 921 266 88 33 16 7 2 
Total 1679 3989 1134 370 152 75 29 6 
April 
3 2 1603 3068 868 282 119 59 22 4 
3 38 463 134 kh 16 8 3 1 
Total 1641 3531 1002 326 135 67 25 5 
4 2 1603 3068 868 282 119 59 22 hk 
3 156 1882 543-180 67 33 13 3 
Total 1759 950 1412 62 186 92 35 7 
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Referring back to Table VI it will be noted also that one inch of runoff for each 
of periods 2, 3 and 4 would produce average “first day” discharges of 1603 
cfs, 326 cfs and 5 cfs respectively. To combine these discharges in the ratio 
of the melt potentials listed above, the “ratio of first day discharges” line of 
Table VI is utilized. The “first day” discharge of period 2 is 4.92 times as 
large as that of period 3, and the “first day” discharge of period 4 is 0.0153 
times that of period 3 for equal runoff, so the “first day” discharge for the 


relative hydrograph will be a combination of discharges from periods 2, 3 and 
4 in the ratio: 


4.92 x .80: 1x .36: 0.0153 x .01, or 
3.94 : 0.36 : 0.000153. 


For convenience in determining the relative hydrographs the runoff for 
period 2 was assumed equal to one inch in all cases. Thus the relative “first 
day” runoff for period 2 is 1603 cfs. Proportionally: 


1603 _ period 3 discharge 
3.94 0.36 


or Period 3 discharge = 04 (0.36) = 147 cfs. 


and for period 4 the “first day” discharge is: 


1603 
3.94 (0.000153) = 0.062 cfs 


Period 4 computations are carried to at least two significant figures because 
of the extremely large increase in second day discharges for this period as 
compared to the first day discharge. Thus, the relative hydrograph for March 
25, 1943 is composed of first day average discharges for periods 2, 3 and 4 
of 1603, 147 and 0.06 cfs respectively, as shown in Table VIII. For period 2 
the average discharges for succeeding days are: 3068, 868, 282 cfs for the 
second, third and fourth days respectively as determined directly from Table 
VI, since it was assumed that runoff for period 2 was equal to one inch in all 
cases, For period 3 the discharges for succeeding days are “f” times the 
tabulated discharges for period 3 given in Table VI, where: 
147 

{= 326 0.450. 
Therefore: 0.450 x 3937 = 1772 cfs for the second day, and 0.450 x 1137 = 512 
cfs for the third day, etc. as shown in Table VII. For period 4 the same 
process is repeated, with f = 0.0126, and the discharges for succeeding days 
are 46, 19, 6, 2 and 1, as shown in Table VII. When the above discharges are 
summed up by days for periods 2, 3, and 4, the relative hydrograph for March 
25, 1943 is obtained as indicated by the “Total” line. Thus the relationship 
between daily discharges as a result of melting occurring on March 25, will 
be in the proportions given by the relative hydrograph. 

The “f” values determined above are equal to the runoff contribution of 
each period to the relative hydrograph. The total runoff represented by the 
relative hydrograph, therefore, is: 1.000 + 0.450 + 0.013 = 1.463 inches. It 
will be noted that the ratio of runoff per period, 1:0.450:0.013 is exactly 
equivalent to the ratio of melt potential per period, 0.80 : 0.36 : 0.01 which 


J 
yt 
a 
i 
‘ 
bd 


ASCE SNOWMELT RUNOFF 1834-23 


was determined originally from the melt potential variation for March 25. 
Tables VII and VIII list the complete computations for relative hydrographs 
for the year 1943. 


Hydrograph Summation 


The process of hydrograph summation, by means of which the daily contri- 
butions to the total seasonal snowmelt hydrograph are determined, can also 
be explained best by the presentation of an example. Therefore, the trial and 
error procedure is presented below as it was computed for the year 1943. 

The average daily discharge for each day of the snowmelt season was tabu- 
lated in sequence as shown in Table IX, Trial No. 1. These discharges are: 
11, 24, 58 and 204 for March 23 through 26, respectively. Base flow dis- 
charges were assumed in accordance with principles not discussed herein. 
Since a base flow of 8 cfs was assumed for March 23 and the total flow for 
that date was 11 cfs, there remained 3 cfs which was presumed to have been 
supplied by melting snows. Since this melting occurred in Period 2 (Table 
Vil), discharges for subsequent days were proportional to the discharges for 
period 2 listed in Table VI and were 6, 2 and 1 cfs for the second, third and 
fourth days, respectively. These values placed in Table IX left a deficiency 
in the discharge for the following day, March 24, of 12 cfs, which was pre- 
sumed to have been contributed from snowmelt. Subsequent discharges, pro- 
portional to period 2 discharges in Table VI of 28, 8, 3 and 1 cfs were comput- 
ed and listed in Table IX. The deficiency in discharge for March 25 was then 
23 cfs, which presumably resulted from melting. The relative hydrograph for 
melt occurring on the 25th is given in Table VIII, as discussed previously. 


Since the first day’s flow for the 25th of March was 23 cfs, the second day’s 
discharge was: 


23 
1750 (4886) = 64 cfs, 
and the third day’s discharge was: 
23 
1750 (1399) = 18 cfs, 


the values 4886 and 1399 being the second and third day’s relative hydrograph 
values as given in Table VIII. 

The process outlined above was continued until April 1, at which time the 
summation of discharges was found to exceed the recorded discharge by about 
558 cfs. This excess was presumed to be the result of an error in the as- 
sumed base flow distribution, or an error in the average daily discharge data. 
The assumption was made that the discharge value most likely to be in error 
was the discharge corresponding to the last day of ice effect, which was 
March 28, 1943. The discharge for March 28, used in the first trial, was 1066 
cfs, which was the discharge computed from the gage reading for an ice index 
of 1.37 feet. Experience with the method indicated that an increase in the dis- 
charge for March 28, of about 34 cfs might result in a balance and a second 
trial was begun using a discharge of 1100 cfs for this date. 

Results of the second trial are also shown in Table IX. It will be noticed 
that all discharges from melt which occurred before March 28 were unaffected 
by the change under consideration and could, therefore, be lumped together as 
a summation of previous discharge. The second trial resulted in a discharge 
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summation of 2394 cfs for April 2, which was 584 cfs more than the recorded 
discharge for this date. It was evident that an over-correction had been ap- 

plied to the discharge of March 28, and a discharge of 1085 cfs was selected 

for this date for the third trial. 

In the third trial the discharges were reasonable until April 6, at which 
time the summation of discharges was 31 cfs too large and a fourth trial was 
undertaken, based on a discharge for March 28 of 1084 cfs. 

The change of one cfs in the assumed March 28 discharge was an over- 
correction and the fourth trial resulted in excessive discharge for April 5 of 
75 cfs. Since it was impractical to apply a further correction to the discharge > 
of March 28, the discharge for the second day resulting from melt which oc- 
curred on March 30 was changed from 7535 to 7545 and the fifth trial was 
computed which balanced out satisfactorily. It was assumed that all discharge 
required for balancing after April 7 was the result of base flow and rainfall 
runoff. 

The flow of 1084 cfs for March 28 is comparable to an ice index of 1.11 
feet. Discharges prior to March 28 were computed using an ice index of 1.37, 
which was determined by actual measurement on March 24. The ice index 
was zero after March 28. 

To compute the runoff attributed to any day, the ratio of the “first day” dis- 
charge from Table IX to the “first day” relative discharge from Table VIII 
was multiplied by the total runoff represented by the relative hydrograph. 

Thus, on March 25 the runoff was: 


23 
1750 


and on March 30 the runoff was 


2257 

2196 (2.180) = 2.241 inches 

After the daily runoff values were computed they were added cumulatively 
from the beginning of snowmelt as recorded in Table VII for the year 1943. 
The positive melt potentials were then added cumulatively for a period begin- 
ning two days before the occurrence of snowmelt runoff, as listed in Table 
Vil. 


(1.463) = 0.02 inches, 


Precipitation 


The occurrence of precipitation during the snowmelt season increases the 
difficulty of the melt potential-runoff studies many-fold. This is the results 
of several factors. First, there is no recording precipitation gage within the 
Big Eau Pleine River Basin. Therefore, the precipitation data, used from 
adjacent stations, may be grossly erroneous. Secondly, records of snowfall 
and snow on the ground are available only after October, 1949. Thus it is 
often difficult or impossible to determine whether the recorded precipitation 
was rainfall available for immediate runoff, or snowfall for which a signifi- 
cant melt potential would be required in order to produce runoff. 

Significant amounts of rainfall occurring in the snowmelt period were 
handled in two ways. Whenever possible the rainfall was added directly to the 
melt potential factors for the same period in the analysis. Hourly precipi- 
tation values were used for the periods for which they were available, and 
total precipitation was divided equally for the day for those periods where 
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hourly data were not available. For some years this procedure would not 
permit a solution and it was necessary in these cases to consider that all of 
the rain falling in a given period ran off. Accordingly, the hydrograph of daily 
values was calculated for each increment of precipitation and their sums sub- 
tracted from the total runoff to leave the runoff attributed to melting snow. 
The years 1942 and 1944 were omitted from further consideration because the 


precipitation during the snowmelt season was greater than the amount of snow- 
melt runoff. 


Correlation of Results 


As mentioned earlier, if the techniques developed herein are to be con- 
sidered valid, a reasonable correlation must be evident between the runoff and 
the melt potential. When the cumulative melt potential for the snowmelt 
season is plotted against the cumulative runoff, as tabulated in Table VII for 
1943, a marked similarity among the curves for the various years is noted. 
Fig. 8 shows the curves plotted with the occurrences of significant initial run- 
off at zero on the abscissa scale. It will be noted that the approximate maxi- 
mum initial rate of runoff is 2 inches over the watershed area for the first 
1.5 inches of melt potential, while the minimum initial rate is approximately 
1.35 inches of runoff for the first 1.5 inches of melt potential. In all cases the 
steepest runoff gradient coincides with the earliest portion of runoff, as would 
be expected, because the snow cover is most extensive in the early melt 
season. For most years the melting is complete, or nearly so, at a melt po- 
tential of 7 inches beyond the occurrence of significant initial runoff. If the 
curves of Fig. 8 are reconstructed with a common ordiante at a melt potential 
of 7 inches the curves of Figs. 9 and 10 result. The seven years having the 
most reliable data are plotted in Fig. 9 and a very good correlation is evident. 
In Fig. 10 the years of less reliable data are plotted. The snowmelt-runoff 
equation, represented in Figs. 9 and 10 by the small circles, is: 


Y = 27.4% - 1.88x2 


where Y is the per cent of the runoff at a melt potential of seven inches and 
X is the cumulative melt potential in inches after the beginning of the runoff. 
It was noted from a study of the melt potential runoff data that initial run- 


off did not occur at a fixed accumulation of melt potential, but varied within 
rather wide limits. 


CONCLUSIONS 


1. The relationship between the cumulative melt potential and the total 
snowmelt runoff, taken as the runoff at a melt potential of seven inches, 
may be expressed by the equation: 


Y = 27.4X - 1.88x2 


for the Big Eau Pleine Basin, where Y is the per cent of the runoff at a 
melt potential of seven inches and X is the cumulative melt potential 
after beginning of runoff. 

2. Apparently, factors causing snowmelt may be measured with less than 
ideal precision at locations somewhat removed from the point of appli- 
cation without losing their effectiveness. This may be due largely to 
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the compensating effect of the various errors which are encountered in 
Hydrologic data. 

3. The general problem of ice effect on discharge rates in the spring 
breakup requires further study. The assumption seems valid that actu- 
al discharge rates for the Big Eau Pleine River at Stratford were higher 
than the rates published by the U. S. G. S. in its Water-Supply Papers; 
however, the ice index method presented must be subjected to further 
study on other rivers before any sweeping claims may be made for its 
validity in other areas. 

4. Rain falling on the snowpack early in the snowmelt season has not pro- 
duced rates of runoff as high as those produced by snowmelt alone, for 
the period of record examined. It must be emphasized, however, that 
rainfall has produced higher rates of runoff than snowmelt in this basin, 
particularly in the fall of the year. 

5. The occurrence of rainfall in the snowmelt season disrupts the normal 
snowmelt-runoff relationships. This phenomenon requires further study 
and perhaps more intensive instrumentation on small plots. 

6. It appears that the coefficient K in Light’s (1944) equation: 


D = KV(T-32°) 


has a reasonable value of K = 0.002 for the melt conditions of the 
Wisconsin river valleys. 

7. The nomograph presented in Fig. 3 is a useful tool in determining the 
snowmelt caused by insolation, and is as accurate as the cloud cover 
data which is used in it. It is believed applicable to latitudes from 40° 
to 50°. 

8. The snowmelt unit hydrograph developed from snowmelt runoff may be 


considered equally applicable to the study of snowmelt runoff and rain- 
fall runoff. 
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CHARACTERISTICS OF A LARGE THROATED SIPHON@ 


Closure by J. C. Stevens 


J. C. STEVENS,! M. ASCE.—The writer is grateful for Mr. Rice’s very 
illuminating discussion. It is particularly valuable since the laboratory model 
studies were made by him, at Oregon State College. 

Mr. Rice has resurrected the ghosts of the definitions of discharge coef- 
ficients and efficiencies of a siphon. Webster’s International Dictionary de- 
fines the latter thus: “1. Efficient quality, power or action. 2. Mech. The 
ratio of the energy or work that is got out of a machine, a storage battery, an 
electrical distributing system, or the like to the energy put in.. .” 

In the writer’s earlier writings (footnotes 4 and 5) he thought he could 
settle these moot questions once for all but perhaps that was expecting too 
much. Mr. Rice concedes that the writer’s definition of efficiency there given 
is the “least objectionable, even though it varies with varying heads on the 
structure.” But how many devices have a constant efficiency? The writer 
can think of none. 

The discharge coefficient will vary with each change of cross-section. For 
this reason it seemed desirable to fix on one, and the outlet was chosen. It at 
least offers means of comparing siphons under comparable conditions—for 
instance, under maximum and minimum heads. 

Testing a model under full head is an excellent method of forecasting the 
characteristics of a projected prototype. Had the laboratory model been 
strong enough for full-head operation much more valuable data could have 
been secured in advance of construction. Table III indicates some of the po- 
tentialities in this method even crude as the facilities were for testing and 
coupled with the weakness of the model itself. Such tests render unnecessary 
the use of controlled atmospheric cabinets for sub-atmospheric pressure 
models—such cabinets have been built. 

The writer cannot agree that the method of correcting the pseudo-energy 
head of Fig. 5 is necessarily erroneous. Mr. Rice indicates that the maximum 
value of the alpha coefficient for velocity head is 2.0 and that any higher values 
would imply a reversal of flow. The writer firmly believes that that was ex- 
actly what happened. The highest value in Table VI is 3.5 at the inlet. Here 
the velocities were moderate but the sharp right-angle turn the water had to 
make to enter the siphon (Fig. 1) produced marked eddies with pronounced 
reversals of flow. A value of 3 is indicated at piezometer sets 3 and 4. No. 3 
is just below the sharp curve at the end of the crest. Here again marked 


a. Proc. Paper 1198, April, 1957, by J. C. Stevens. 
1. Cons. Engr., Stevens & Thompson Engineers and Leupold and Stevens 
Instruments Inc., Portland, Ore. 
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eddies with flow reversals surely obtained and doubtless these continued on to 
Set No. 4. 

The writer will agree that this method of constructing the energy line for 
water passing through a tortuous siphon has elements of arbitrariness in it. 
However there are two facts that cannot be ignored. (1) The energy line must 
slope continuously downstream since energy is being lost continuously, and 
(2) there is hardly any point within the siphon where the true mean velocity 
equals Q/A, hence the coefficient alpha is everywhere greater than unity, even 
though Table VI shows it as 1.00 at piezometer set No. 2. This is in the 
middle of the long crest where it would naturally be nearer unity than any- 
where else in the siphon. 

, Mr. Rice’s velocity traverses of flow in the laboratory model cannot possi- 
bly represent what actually occurred in the prototype, with its violent eddies 
and direction reversals of velocities. It is possible however that had it been 
practicable to have made such traverses in the high head model a greater 
conformity to the prototype might have resulted. 


ERRATA 


Paper 1198 as printed in HY 2 April 1957 contains errors as follows: 


Page 1 under synopsis 3rd line change primary to priming 
2, line 16 change floods to boards 


line 10 from bottom « change cro cross-section to profile 


5, 3rd line after the table change Figs. 1 and Figs. 1 and 2 to Fig. 2 
11, 1st line of the table change 56 in date to 55 


Eq. (6) change final zero to letter oO 


12, Table IV heading col. (1) ) change Mt to m 
15, Table V Ist col. change 6 to 60.0 — 
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TRANSITION FROM LAMINAR TO TURBULENT FLOW IN A PIPE 


Discussion by J. M. Robertson 


J. M. ROBERTSON,! M. ASCE.—The subject of laminar-turbulent tran- 
sition has received little attention in ASCE publications, although it has been 
discussed in the general fluid mechanics literature since the time of Osborne 
Reynolds and Lord Rayleigh. This lack of concern of the civil engineer stems 
from two main causes. Firstly, most theoretical approaches have been in re- 
gard to problems (such as boundary layers and laminar-flow airfoils) of 
primary concern to the aeronautical engineer. Secondly, the experimental 
information on transition of primary interest to the civil engineer interested 
in pipe flow was obtained over half a century ago. In fact, Reynolds’ classic 
experiments were conducted seventy-five years ago, while over a hundred 
years ago G. H. L. Hagen noted a temperature effect on the occurrence of the 
transition between, what may now be termed, laminar-viscous and turbulent 
flow.2 Although Reynolds did not use our present day terminology, his studies 
did indicate that for “engineering pipes” the critical Reynolds number has the 
(now) commonly accepted value of 2000 to 2500. This agrees with Hagen’s 
test results (for pipes with square entrances) which have been shown to indi- 
cate a critical value of 2200 to 3000. Since the subject appears to be rather 
hoary, this writer will attempt to justify current interest in the transition 
problem. 

It is usually stated that there is an upper critical Reynolds number above 
which the flow is always turbulent. Actually there is some uncertainty in this 
regard. Experimentally, it is observed that as we decrease the disturbance 
entering with the flow, the critical number increases rapidly. Reynolds him- 
self, in what we now call a “Reynolds apparatus” (transparent pipe with bell- 
mouth entrance from a rather large tank), had little difficulty achieving lami- 
nar flow at Reynolds numbers R = VD/V of as high as 12,000 to 13,600. Other 
experimentors have been able to reach even higher values with no great 
trouble. Thus in 1909 using Reynolds’ original apparatus, but taking particu- 
lar care with the bellmouth and allowing the water in the tank to settle for up 
to half a day, V. W. Eckman(3) raised the upper critical number to a value of 
30,000 to 50,000. According to Table III the author has found laminar flow at 
Reynolds numbers of about 15,000 to 51,000 depending upon the spot source 
and other flow factors. These high values are a tribute to his experimental 
skill, since there is some evidence that a timewise developing flow may be 


a. Proc. Paper 1450, December, 1957, by M. R. Carstens. 

1. Prof. of Theoretical and Applied Mechanics, Univ. of Dlinois, Urbana, Ill. 

2. The studies of Hagen, Reynolds and others have been carefully analyzed 
and presented by L. Schiller(1) and L. Prandtl.(2) 
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less stable than a steady flow. One wonders how high a Reynolds number he 
might have reached without turbulence if the disturbance strengths of his spot 
sources had been further reduced. 

The care required to achieve laminar flow at such high Reynolds numbers 
may appear impractical to the hydraulic designer—however the water ina 
reservoir often has ample time to become quiescent. Furthermore, laminar- 
viscous flow may occur even in pumped systems at Reynolds numbers con- 
siderably above the accepted critical values. Thus, we have the results re- 
ported by the University of Iowa 4) with a system in which oil is pumped con- 
tinuously through a pipe; it is impossible to obtain turbulent flow even at the 
upper limit of the pump capacity (R about 7800). The turbulence must be arti- 
ficially triggered to make it appear in this case. It would appear that ap- 
preciable amounts of laminar-viscous flow may well occur in some of our hy- 
draulic structures of clean design. 

The author’s study of transition is therefore of potential value to designers 
as well as being of intrinsic interest to fluid mechancians. One has only to 
glance at the discussions concerning laminar-flow shapes in aeronautical 
publications to visualize the potential gaines. Fifty to one hundred-fifty per 
cent increases in aircraft ranges from the use of suction-stabilized laminar 
boundary layers are suggested. If comparable achievements are possible in 
some of our hydraulic structures, it might be possible to justify the special 
construction and precautions needed to achieve the results. Only through 
fundamental studies such as reported by the author can one obtain the 
knowledge needed to assess what is possible. Also, we need to be able to tell 
if our pipe flows may occur with appreciable laminar flow regions. In many 
cases one may not design for laminar flow but one should consider its likeli- 
hood of occurrence. Differences in terms of pressure drop, heat transfer, 
etc., between a design based on the conventional assumption of turbulent flow 
(at Reynolds numbers above 2500) and a flow with a considerable laminar- 
viscous region can be significant. 

Study of the laminar-turbulent transition in pipes is of considerable in- 
trinsic interest, since there is some question whether it should occur at all if 
there is no disturbance. As noted above, the experimental evidence shows no 
upper limit to the critical Reynolds number as the external disturbances are 
reduced, One therefore turns to theory to see what says about this question. 
The author has discussed his observed transition results in terms of the pre- 
diction of the small-disturbance theory for the stability of Tollmien- 
Schlichting waves, but for the case of Blasius flow, i.e., for the laminar- 
viscous boundary layer on a plane with zero pressure gradient. The result is 
certainly suggestive, even if not applicable to the present problem. This 
small-disturbance theory has been quite useful in predicting the instability of 
laminar-viscous flows. The instabilities and resultant wave length—Reynolds 
number relations (cf Fig. 9 of the paper) do not immediately indicate tran- 
sition to turbulence, but they have been found to exist. Although such pre- 
dictions were long considered of doubtful physical reality, they were first veri- 
fied for flat plate flows in 1940. The phenomenon of selective amplification of 
disturbances of the proper wave length and frequency in laminar-viscous flows 
is now accepted as a prerequisite to transition, except where relatively large 
disturbances are present to override this occurrence. 

Prediction of a critical Reynolds number for pipe flows on the basis of the 
small-disturbance theory has been a rather elusive thing, considering the 
practical interest in transition and the available experimental information. 
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Amongst the early analysis was W. M. F. Orr (1906) who utilized the ap- 
proaches of Lord Rayleigh (1879 et seq) and O. Reynolds (1895)—in fact the 
basic relation of the small-disturbance theory is often called the Orr- 
Sommerfeld equation. Orr found a critical (for instability of the wave dis- 
turbances) Reynolds number of 117 for pipe flow. However, more refined 
analyses have indicated this result to be inaccurate. Thus, C. C. Lin(5) state: 
“no theoretical indication of instability has ever been found for the fully de- 
veloped parabolic profile in the case of the circular pipe.” 

Only recently has the apparent dichotomy, between our expectation of a 
finite critical Reynolds number for pipe flow and the theoretical prediction, 
been explained. The clarification stems from consideration of the fact that 
fully developed flow (parabolic velocity profile) in the pipe does not appear 
full blown at the entrance. Instead it develops in an entrance region whose 
length depends upon the mean flow Reynolds number. About 0.06 R diameters 
is required for the entrance length; thus 120, 600 and 3000 diameters are re- 
quired for the flow to develop at Reynolds numbers of 2000, 10,000 and 50,000, 
respectively. In this entrance region the flow is of a boundary-layer type with 
a small favorable pressure gradient. Because of this boundary layer type 
flow the author has some justification for comparing his results with the sta- 
bility predictions for Blasius flow. For this latter case the most refined 
calculations indicate a critical Reynolds number (lowest R value on neutral 
stability curve of the author’s Fig. 9) of R5+ = 420. The boundary layer of the 
entrance-region flow differs from the Blasius boundary layer due to the 
pressure gradient. 

The possibility that the entrance-region boundary layer may be unstable 
even though the following Boiseuille flow is stable was investigated by T. 
Tatsumi(6) in 1952. In terms of his solution for the velocity profiles of the 
entrance-region boundary layer, his calculations yielded neutral stability 
curves, i.e., indicated the existence of flows unstable to certain disturbance 
waves, The minimum critical Reynolds number (R = VD/v) was found to be 
19,400 at a station rather near the entrance, namely at about x/D = 0.0004 R 
(0.67% of the length of the entrance region) or about 7 D from the start of the 
pipe. For flows at higher Reynolds numbers there is a portion of the entrance 
region which is unstable and in which small disturbances, if of the proper 
wave length or frequency, will be amplified. Thus for a Reynolds number of 
25,000 the flow is “unstable” (capable of amplifying certain small dis- 
turbances) for x/D distances between about 1.4 and 20 as may be seen from 
study of Fig. A of this discussion. If not enough amplification of the dis- 
turbances occurs in this distance and they do not become large enough to 
cause transition to turbulence, then Tatsumi’s analysis would say the flow is 
stable. It is apparent that the higher the flow Reynolds number (above 20,000) 
the more opportunity there is for amplification of small disturbances and the 
greater is the possibility that transition to turbulence will occur. A study of 
the author’s transition data, in terms of the parameters of Fig. A, indicates 
that the transition Reynolds numbers were generally above the minimum criti- 
cal value for stability and tended to increase with x/D. However much larger 
critical x/D values occur. 

The author discusses the fact that the disturbances must be considerably 
amplified after they pass the critical Reynolds number (for neutral stability) 
before turbulent flow sets in. His discussion is in terms of the stability analy- 
sis of the Blasius flow rather than the actual case in hand. It would appear 
more appropriate to compare the experimental results with the predictions of 
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the stability analysis of Tatsumi even though this is only for the steady flow 
case. Actually the comparison should be with neutral stability curves for the 
profiles calculated from Eq. (3) for the temporally and spacially growing 
boundary layer. However, neutral stability curves are not easily determined. 
Since the stability condition is dependent upon the shape of the velocity pro- 
file, a suggestion as to which known set of solutions might best apply should 
be apparent through comparison of the velocity profiles. Lacking this infor- 
mation, it is still possible to glean some useful information from comparison 
with the Blasius case solution. The remainder of the writer’s remarks will 
follow the author in this regard. . 
The author’s presentation of his transition data in Fig. 10 is interesting and 
suggestive. However there are several points about the possible neutral sta- 
bility curves shown that the present writer cannot follow. These points have 
to do with the wave length A chosen as the basis of calculation for the curves. 
Thus 2.1 ry is chosen for spot 2 as “the longest disturbance wave length which 
originates at the pipe junction.” We know very little about the discontinuity 
presumed to exist at the pipe junction and if we did we still could not predict 
anything about the wave length of the disturbance which it would produce. 
Hence, the present writer is forced to conclude that the author chose a wave 
length which would pass the possible neutral stability curve through the ex- 
perimental point having the highest transition Reynolds number. The justifi- 
cation for this assumption is that transition is not expected once the dis- 
turbance has passed the upper branch of the curve. It is also stated that the 
appearance of some of the points (actually seven of the nine points for spot 2) 
below the lower branch of the curve causes no conflict. It is hard to accept 
this statement, since for points in this region the disturbances would have had 
no amplification whatsoever. There seems to be no reason why the wave 
length being amplified should be a unique function of the spot source. Actually 
it is a function of the boundary layer flow as shown by the neutral stability 
curve in Fig. 9. 
The amplification situation for the simpler Blasius case is clarified through 
reference to Fig. B, which is nothing but the author’s Fig. 9 with some ad- 
ditional information. Some of the amplification contour lines (as found by 
S. F. Shen(7)) mentioned by the author have been added, together with straight 
dashed lines indicating the progress of constant wave length disturbances 
across the amplification ridge. The encircled numbers on some of these 
dashed lines indicate approximately the total amounts of amplification that 
have been calculated for these certain wave lengths. The trends of the lines 
of constant A are quite suggestive of the run lines on Fig. 10. They are in- 
clined straight lines which would pass through the origin, if extended, since as 
the wave progresses along the flow 6* increases and the ordinate and abscissa 
values increase in constant ratio. Natural transition in the low disturbance 
flows (free stream turbulence intensity less than 0.3 per cent) is found at Rx« 
values above 2000. In the few cases in which natural Tollmien-Schlichting 
waves have been measured in terms of their frequency and wave length charac- 
teristics, they have been found to occur approximately along the upper branch 
of the neutral stability curve at Reynolds numbers above 2000. It is apparent 
that these have been subjected to amplifications of several hundred times at 
least and that they must cause transition rather quickly since damping sets in 
after they pass the upper branch of the neutral stability curve. It would ap- 
pear that amplifications of one hundred or less are inadequate to lead to 
transition—the amplification occurring at the critical Reynolds number of 420 
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is small indeed. It is not surprising that we find a considerable gap between 
this critical Reynolds number and that for transition to turbulence. 

Attempts to fit the above concepts and results for small disturbance ampli- 
fication and transition to the results shown in Fig. 10 are not encouraging. 
Two possibilities suggest themselves as explanation for the poor agreement. 
The first is that the disturbances causing the transition are not small. The 
second is that the neutral stability curve for Blasius flow is a poor repre- 
sentation for the actual flow situation. This latter is beyond doubt true, but it 
behooves one to also check on the former possibility. Information on the ef- 
fect of roughness elements in causing premature transition has been well 
summarized by H. L. Dryden. (8) Treating the joint discontinuity as a flat 
strip element, the important factor telling us whether it will trigger transition 
is the ratio of the height to the boundary layer disturbance thickness 6* evalu- 
ated at the location of the roughness element. For turbulence spots 2 and 3, 
calculations based on the data of Table III indicate that at transition this ratio 
was of the order of 0.03 to 0.07. Since the change in boundary layer thickness 
could not have been large in the distance from the pipe joint to the transition 
location and since Dryden shows that no appreciable effect occurs for ratios 
below 0.1, it is concluded that indeed the joints were only sources of “small 
disturbances.” In the case of the source for spot 1, the situation is quite 
different although exact calculation is not possible with the available infor- 
mation. The sand roughness initiating this spot was on the entrance bellmouth 
where the displacement thickness is extremely thin. The only question in this 
case is why the author chose to complicate his study with this avoidable source 
of disturbances and, perhaps, why the sand did not always lead to immediate 
transition. 

The author has been very successful in ferriting out of his rather indirect 
study considerable information on the occurrence of transition. It is unfortu- 
nate that it was not possible to measure the velocity directly. J. Rotta(9) was 
able to use a hot-wire anemomter to record the velocity fluctuations for tran- 
sition studies with air flow at Reynolds numbers up to 18,000. No evidence of 
Tollmien-Schliching waves appear in his results, probably because the hot- 
wire was located only at the end of the pipe. It would be most enlightening to 
see how these appear in the regions where their amplification might be ex- 
pected. Rotta’s measurements do give some interesting indications of an 
intermittent type of laminar-turbulent flow such as noted and discussed by 
Prandtl and Tietjens. (2) 

The writer finds himself at odds with the author in regard to Fig. 7 and the 
discussion thereof. The turbulence cannot be “diffused into the laminar 
boundary layer,” since it comes from a breakdown of the laminar-viscous 
flow. Such a breakdown, as opposed to an outside effect, is the whole point of 
the small-disturbance theory to which we have been trying to tie the transition 
occurrences. Fig. 7 furthermore shows the turbulent flow after transition as 
completely filling the pipe cross section where the laminar boundary layer 
flow did not. This is not necessarily so. Assuming that the laminar boundary 
layer profile differs only in detail from that for the Blasius case we may take 
6 = 2.9 6* and, except for run 20, 6 is therefore less than one-quarter of the 
pipe radius. Although the turbulent boundary-layer may well become four 
times the laminar one in thickness, this does not take place immediately. In 
flat plate flows the length of the transition region is often found to be as great 
as that of the preceding laminar flow region. Apparently the simplification of 
Fig. 7 is satisfactory for the purpose of interpreting the data. 
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In Table III, it appears that the boundary layer thickness and Reynolds 
number occurring for transition of spot 2 is less than that for spots 1 and 3 
(except for run 20). The author’s explanation for this would be appreciated. 
In terms of the indications that the disturbances which are amplified and then 
breakdown to turbulence starts from Rg* values of about 1000 and hence 
2 776*/ values of about 0.17 the author’s discussion of the number of wave 
lengths between the disturbance source and transition would seem to be in 
error. For spot 2 it would appear that the distance is only 3 to 15 wave 
lengths rather than 9 to 40. Also is not the x/D value listed for spot 1 or run 
20 in error? Finally, it seems impossible that for spot 1 of run 21, the 


boundary layer can have a finite thickness at x/D = 0, so this value must also 
be in error? 
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FLOOD FREQUENCIES DERIVED FROM RAINFALL DATA@4 


a Discussion by C. O. Clark 
| = Closure by J. L. H. Paulhus and J. F. Miller 


C. O. CLARK,! A. M. ASCE.—The criteria for designing structures, dams 
and bridges, subject to rain and floods can be derived from a look at a small 
budget from an easy chair, or from the arduous, costly, and sometimes un- 

# comfortable investigation of all data at and all around the scene in all kinds of 

4 weather, or from any intermediate combination of observation and non- 

. 3 observation. The consequences of the former, although prevalent, are cata- 

F strophic. The mere fact that the subject paper deals with flood frequencies 

a6 derived from something other than floods, and from data outside of the water- 

_ shed to which it refers, invites a critical examination as to whether it has 

; a really produced the flood frequencies. 

o The authors have used rainfall data near a site, but outside the subject 

, watersheds, to derive a series of calculated flood flows, which calculated 

ad flows are not greatly different in their mean value and variation from, and 

q hence indicate about the same frequency of flooding as, the actual peaks of that 

aa same short period. But to reconcile their estimate of normal flood frequency 
with one observed flood, the authors have to conclude, “data suggest that the 

1933 observed flood peak of 11,200 cfs is an unusual event with a return period 

probably much greater than 100 years”. 
4 The writer finds instead that the 11,000 cfs flood is not at all rare, and that 
4 it is only luck that there were not two or three as large in the sample period. 

7 This fundamental difference in views is attributed to faults of a common sta- 

= tistical process by which judgments about frequency are reached by exami- 

a nation of a small number of observations, and which in this case was probably 

used with both sets of data, the set belonging to the stream and the set derived 

from rainfall. 

The statistical process is understood to be applicable to some uses with a 
representative sample of a large population, and has been thought to yield 
reasonable estimates of size of floods of rare flood frequency when the sample 
_— . was fairly large and nothing indicated that it was a poor sample, or non- 

- representative, of the whole. 
In the record of the subject stream, Codorus Creek, there are several 
reasons to think that the record is a non-representative sample, some of 
‘ which reasons involve examination of paths of great storms along the Atlantic 
coast and noting how few hit Codorus Creek as compared with other places. 
Hence the data are neither in the rainfall nor the flood experience of a normal 


a. Proc. Paper 1451, December, 1957, by J. L. H. Paulhus and J. F. Miller. 
r 1. Cons. Hydr. Engr., Tulsa, Okla. 
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sample of storm events; both kinds of data, then, are “biased” or unrepre- 
sentative. Neither does the trend or scatter of the observed data suggest that 
the sample may be successfully fitted with a probability curve of a standard 
type derived from these data alone. This, too, may be correctly said of the 
record of the stream which the writer introduces to demonstrate that the 
normal flood experience of the region, and probably of Codorus Creek, is 
much larger. But neither has this record been used by the writer to select 
the form of frequency curve which he says is regional in application and ap- 
plies to Codorus Creek, too. Of the records, either of Codorus Creek or the 
one introduced, neither is of a sample size sufficient to disprove the form of 
the applied or indicated frequency. All that is shown is how much small 
samples can vary from a curve without proving the curve is wrong or right, 
and without being capable of reliably indicating what curve is right. 

Among the streams of this general area is another stream which from an 
identical quantity of runoff in six hours will produce the same size of peak as 
the authors assign to Codorus Creek. It is the Passage Creek near Buckton, 
Virginia, some 90 miles SW of the Codorus Creek, but well within the general 
area of the problems and basins selected by the authors. There is one signifi- 
cant difference. This stream had its floods; Codorus Creek did not. At least, 
not in the period 1933-54. 

A comparison of the flood experience of these streams with similar flood 
capacity, but different accidental exposure to opportunity for collecting a 
record during a specific period, is shown on Fig. 1. The circles and crosses 
are the experience, and the only facts. The lines are methods of portraying 
the normal experience of an infinite period of time in relation to these 
samples. The line attributed to the authors is plotted from Table 4. The lines 
attributed to the writer are from flood frequency studies made in the subject 
area about 1942, by James B. Saylor and the writer, for the purposes of the 
Potomac River Flood Control Survey Report, the essential results of which 
were published by the District Engineer, Washington District, Corps of Engi- 
neers, in his report on that subject to the Chief of Engineers in April 1944. 
This application is in accord with that study and flood producing capacity of a 
stream of 87 sq. mi. of drainage area and concentration time of 6 hours. 
(Data to this effect appear on p. 1484 of Transactions ASCE, 1945, Vol. 110.) 

The authors give 74.3 sq. mi. as the watershed of Codorus Creek, and Fig. 
5 indicates a time of concentration of about four hours. So the slightly larger 
drainage area of Passage Creek is to some extent offset by the apparent, but 
surprising, greater “flashiness” of Codorus Creek. In keeping therewith, the 
authors have found Codorus Creek capable of producing a flow of 6700 cfs 
from 1.0 inches of surface runoff uniformly distributed in six hours, and the 
writer had estimated that Passage Creek could produce 6900 cfs from such an 
event. On this basis, the streams are selected as essentially of the same flow 
producing capacity, or as having flood frequencies no more different than the 
two almost coincident lines on Fig. 1. 

Now let us return to the evidence that the larger flood of Codorus Creek is 
rare. Or to the evidence that because it lies well above a line fitted to the 
remaining or great majority of points that it should be so regarded. 

Except for the two greater values of the flows for Passage Creek, the ma- 
jority of those data are no greater than the majority of the values for Codorus 
Creek, and there, too, a line fitted to the majority of the data and extended 


would indicate a phenomenal rarity to the larger two of the actually observed 
items. 
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The writer has suggested that there is a pitfall which waits for all of the 
unwary in frequency calculations from short records, namely that there is no 
record sufficiently normal to be used, nor a sufficient number of observations 
to compensate for errors induced by the abnormalities of a few occurrences. 
What is there to indicate an abnormality in the record of Codorus Creek? 
First of all, a comparison like that here discussed. But what is there to indi- 
cate that it is Codorus Creek and not Passage Creek, or the writer’s idea of 
normal frequency, that is abnormal? Part of that is in the review of the 
record of this area as a whole. 

Such a review shows that few of the great storms of 1930 to 1950 passed 
very close to Codorus Creek. A further review shows that a storm prior to 
1900, probably 1889, not only passed nearby, but created rumors of great 
floods that persist to this day. Perhaps it was the same storm which sent the 
Potomac so high that backwater in the gorge below the Falls backed right over 
the Falls and submerged Delcarlia dam above the Falls and established a high- 
water mark, there engraved in rock (it was a mark that was discarded as 
falsely high until a lesser flood of September 1944 demonstrated that the Falls 
could, indeed, be drowned out, whether the dam could be or not). 

Continuing in such a review, the 11,200 cfs flood of Codorus Creek occurred 
in August 1933, in general a period of drought, but in a storm that established 
records on a number of nearby Maryland and Pennsylvania streams with short 
records, but not on any of the streams with moderate length of record at that 
time. The storm is listed as NA 1-24 among the great storms studied as 
“Storm Rainfall in the United States”, published by the War Department, and 
one of the bases for preparing Hydrometeorological Report No. 23, a report 
of the U. S. Weather Bureau on “Generalized Estimates of Maximum Possible 
Precipitation Over the United States.” 

The storm maps show 6.8 inches over 10 sq. mi. in six hours as the most 
intense part of this storm. Report No. 23 gives 24 inches as possible on 10 
sq. mi. both there and near Passage Creek. The writer does not consider the 
1933 rainfall to be a sufficient approximation of the merely possible as to be 
even unusual. How then, can the flood be unusual ? 

For further comparison, the two larger flood items for Passage Creek are 
for October 1942 and March 1936, Storms SA 1-28 and SA 1-27. In the 
October center, 6.0 inches is the value given for the maximum 10 sq. mi. in 
six hours, but it is shown as occurring after 14 inches of prior rain, as com- 
pared with 3 inches preceding the maximum rates in the NA 1-24 storm. 
March 1936 was such a commonly wet month and year that no rare signifi- 
cance can be assigned either to the 2.7 inches of six-hour rainfall or the 
magnitude of the flow which it produced and which is the second larger on 
Passage Creek. Compared with these not too rare storms, however, the 1889 
storm, SA 1-1, is presumed to have deposited 7.4 inches of water in the same 
period and area, a wetter season of the year, and somewhat to the northwest 
of Codorus Creek. 

What would Codorus Creek do with 7.4 inches of water in six hours, if one 
inch of runoff produces 6700 cfs? Fig. 3 indicates that from 2.5 to 5.5 inches 
would run off. Wouldn’t that make 17,000 to 36,000 cfs of flow? What frequen- 
cy do such flows have? Where would they fit on the suggested normal frequen- 
cy curves? 

Out of some 2500 flood items that lie behind the curves which the writer 
presents as normal are over a hundred as large as 2 inches of water in 6 
hours could produce, and a few are over three times that big. For the basins 
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the authors list in Table 1 as Nos, 14-19-20-15-18, values of 5.4-4.1-3.3-2.8- 
2.7 are found in some very short records. In the historical records of the 
same area, but without any guide as to the total number of items among which 
they would be found, are values as high as 8 inches of runoff could produce. 
All are in the area of study of the subject paper, and in the area for which 
Report No. 23 suggests 24 to 26 inches of rain as being possible on 10 sq. 
miles in six hours (Values for 100 sq. mi. are not given). 

It has been the purpose of this discussion to show that something more than 
the authors have presented is necessary to derive flood frequency either from 
rainfall data or from streamflow data, and that without it either or both kinds 
of data can lead to an erroneous conclusion. There is a wide amplitude of 
variation in the magnitude of floods which are possible in the subject region. 
The largest item cannot be discarded from the sample only on the basis of the 
remaining experience in the same small sample, nor for the convenience or 
judgment of the calculator. There is nothing needed to explain the wide vari- 
ations which the writer has presented except measures of accidental proximity 
to the storm centers and accidental variations of prior wetness and coinci- 
dence in the degree of storm intensity with that to which a particular stream 
is most sensitive. These alone are sufficient to explain wide variations in any 
small sample, or to explain the extremely long period of time that it would 
take any one area to experience all of the variations, or even a representative 
sample, of those which are possible. There are just more possible combi- 
nations for hitting a paying combination than there are on a common variety 
of conventional amusement machines, and as many statistical errors to be 
avoided in drawing a conclusion on the frequency of outflow as there would be 
from a similarly small sample of plays on such a machine, 

Since none would draw conclusions from a small number of personal obser- 
vations of action of such a simple machine, what recommends doing so about 
either rainfall or streamflow observations, either by the authors’ procedure, 
or by any other? Instead, let this profession be built on the processes by 
which no single observation alters the conclusion by a determinable amount. 


J. L. H. PAULHUS,! A. M. ASCE and J. F. MILLER.2—There certainly can 
be no argument against the contention that a good relation can be developed 
between any two sets of randomly selected unrelated items if the items in each 
set are arrayed in like or opposite order. This fact is so well known that it 
does not seem to deserve the lengthy treatment given it by Riggs and Benson. 
The effects of the ranking procedure on correlation had already been explained 
only a short time before by a Dutch climatologist, C. Levert, in a discussion 
of *Relation Between Point and Areal Rainfall”, published in the Bulletin of 
the American Meteorological Society, Vol. 38, pp. 618-620, Dec., 1957. How- 
ever, the subject has no bearing on the authors’ paper as correlation was not 
involved—the objectives of the study being the same as those of any other 
method of extrapolating streamflow records, namely, the approximation of 
flood-frequency data and not of individual’annual peak discharges. Other 
comments in the discussion indicated that several important points in the 
original paper had been misinterpreted. 


1. Staff Hydrologist, U. S. Weather Bureau, Washington, D. C. 
2. Meteorologist, U. S. Weather Bureau, Washington, D. C. 
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The second sentence of the discussion implies that a relation was developed 
between synthesized and observed annual peak discharges and then used to ex- 
tend the record. The writers did not develop such a relation at any time. The 
observed and synthesized annual peak discharges were arrayed in like order 
to facilitate comparison for determining whether frequency analyses of the 
two would yield similar results. When it was determined that frequency analy- 
ses of the synthesized peaks would yield results within 20% of those for the 
observed, the record of synthesized peaks was extended by means of rainfall 
data, a rainfall-runoff relation, and a unit hydrograph. At no time was the 
record of synthesized peaks extended by a relation between observed and 
synthesized peaks for a specific period. 

The objection to the ranking procedure and the suggestion that “the authors’ 
method might have been evaluated by applying it to the latter part of a long 
flood record and testing the extended results against the known earlier part of 
the period” are very easily answered. Using the only two basins included in 
the original study that had records of 30 or more years, the writers have per- 
formed the suggested test while showing that the ranking procedure has no 
bearing on the results. 

In the top half of Table 1 are listed in chronological order the observed and 
synthesized annual peak discharges during the 16-year period, 1938-53, for 
the Whippany River at Morristown, N. J., (29.4 sq. mi.). The synthesized peak 
discharges were obtained from runoff computed by applying a rainfall-runoff 
relation for the Anacostia River, about 175 miles SW of the project basin, to 
1938-53 storm-rainfall data from the 2-station network, Plainfield and 
Elizabeth, N. J., roughly 25 miles SE of the project basin. The 6-hour unit 
hydrograph used for converting the runoff to discharge was based on data 
subsequent to 1939, as was the testing and adjusting of the Anacostia River 
rainfall-runoff relation. Application of the Gumbel method of frequency analy - 
sis to the observed and synthesized peak discharges for the 1938-53 period 
yielded the results shown at the right in the top half of Table 1. The synthe- 
sized frequency data are well within the 20 per cent limitation, and the 
rainfall-runoff relation and unit hydrograph were therefore considered adapta- 
ble to the project basin. The arrangement or order of the observed and 
synthesized peak discharges has absolutely no effect on the flood-frequency 
data. The discharges could be ranked in like or opposite order or scrambled 
in any fashion whatever, and the frequency data would remain the same. 

The lower part of Table 1 shows the observed and synthesized peak dis- 
charges for the 1922-37 period and their frequencies. The same rainfall- 
runoff relation and unit hydrograph used in synthesizing the annual peak dis- 
charges for the 1938-53 period were used to synthesize those for the 1922-37 
period, the only difference being that they were applied to the Plainfield- 
Elizabeth storm rainfall data for 1922-37 instead of 1938-53. Again, the 
frequency data for the synthesized annual peak discharges are well within 20 
per cent of those for the observed. 

Table 2 shows the same kind of test applied to synthesized annual peak dis- 
charges for the North Branch Raritan River near Far Hills, N. J., (26.2 sq. 
mi). In this case, a rainfall-runoff relation for the North Branch Raritan River 
near Raritan, N. J., which includes the project basin, was applied to Plainfield- 
Elizabeth (about 30 miles ESE) storm-rainfall data for the 1938-53 period to 
synthesize annual peak discharges for the same period. Data for developing 
the unit hydrograph and testing and adjusting the rainfall-runoff relation were 
all from the period subsequent to 1939. The synthesized flood frequency data 
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TABLE 1 


FLOOD FREQUENCIES FOR THE WHIPPANY RIVER AT MORRISTOWN, N. J. 


Annual Peak Discharges Flood Frequency 
Return 


Period Observed Synthesized 
(cfs) (cfs) 


Year Observed Synthesized 
(cfs) (cfs) 


772 
1259 
1951 938 470 10 1477 1582 
1950 375 380 | 25 1807 1990 
1949 1350 610 
1948 829 900 | 
1947 591 420 
1946 1020 860 
1945 460 690 
1944 640 1150 
1943 540 470 
1942 1220 1460 
1941 679 510 
1940 1500 830 
1939 529 470 
1938 1170 1990 | 
Mean 875 836 H 
0 358 | 0. = Standard Deviation 
Mode 696 614 
1937 491 560 2 661 617 
1936 ae 1170 5 30 983 
1935 7 390 10 120 1225 
1934 620 600 | 25 1484 1531 
1933 1060 750 
1932 485 430 
1931 456 400 
1930 466 420 
1929 830 430 
1928 1100 1320 
1927 394 1140 
1926 585 1130 
1925 680 380 
1924 830 650 
1923 600 430 
1922 700 480 
Mean 704 665 
299 333 
Mode 554 499 
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being within 20 per cent of the observed, the same rainfall-runoff relation and 
unit hydrograph were then used to synthesize annual peak discharges for the 
1923-37 period from storm rainfall for that period. Again, the synthesized 
flood frequency data are within 20 per cent of the observed. 

Nowhere in the two tests just described would any rearrangement of the 
order of the observed or synthesized annual peak discharges within any one 
period have any effect on the computed flood frequencies. 

The results of the above tests were satisfactory, considering the short 
periods of record used for controls and testing—20-year periods should be the 
minimum and longer periods would be preferable. Generally speaking, this 
type of test is not a fair one for evaluating the writers’ procedure because the 
off-basin rainfall network does not measure the basin rainfall storm by storm. 
With the rainfall network fairly well removed from the basin, it is impossible 
to distinguish between the errors of the synthesizing procedure and sampling 
errors in the statistical sense. It would be very interesting to see a regional 
procedure subjected to this same type of test. 

The best test of the writers’ procedure would be to use a basin with a long 
streamflow record, say 40 years or longer, and having a rainfall network for 
5 to 10 of those years. The basin storm-rainfall data and concurrent stream- 
flow record would be used for selecting or developing a suitable rainfall- 
runoff relation and for constructing a unit hydrograph. The relation and unit 
hydrograph would then be applied solely to the storm-rainfall data of an off- 
basin network to synthesize the annual peak discharges for the entire period 
of the basin’s streamflow record. Comparison of the observed with the synthe- 
sized flood frequencies would provide the best possible test of the procedure. 

Incidentally, Messrs. Riggs and Benson misinterpreted the purpose of Fig. 
6, which was to avoid the necessity of having several unit hydrographs for the 
same basin. The computed peak discharges were obtained by applying the se- 
lected rainfall-runoff relation and unit hydrograph to the basin rainfall actu- 
ally producing the observed peak discharges. The observed and computed 
peak discharges were then plotted against each other on a concurrent basis 
and not after being ranked, as should be obvious from the plotted points of 
Fig. 6. 

In the discussion it was stated that an annual flood peak record of about 15 
years is required for applying the writers’ procedure. This is not so. The 
basins used in the study were restricted to those having at least 15 years of 
streamflow record so that the observed annual peak discharges would serve 
as a sort of “yardstick” for evaluating the synthesized peaks, as in the above 
tests. Actually, concurrent measurements of streamflow and of rainfall from 
an adequate network over the problem basin for a reasonable number of moder- 
ate and heavy storms would be the only requirements for application of the 
technique suggested by the authors. The rainfall and streamflow data would 
determine the choice and adjustment of the rainfall-runoff relation to be used 
and permit the development of a unit hydrograph. One or two years of record 
would suffice in many cases although longer records would, of course, be 
preferable. Furthermore, the observations need not be continuous but may 
be made only when major storms or high discharge rates are expected. 

There is absolutely no reason why the synthesized and observed peak dis- 
charges should correspond on an annual basis—no more so than observed peak 
discharges for basins being used in developing a regional procedure. If the 
rainfall network were in the problem basin, agreement between synthesized 
and observed peak discharges might then be expected. The agreement would 
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TABLE 2 


FLOOD FREQUENCIES FOR THE NORTH BRANCH RARITAN RIVER NEAR FAR HILLS, N. J. 


Annual Peak Discharges Flood Frequency 
Observed Synthesized Return Observed Synthesized 


(cfs) (cfs) Period (cfs) (cfs) 


S. 0. = Standard Deviation 
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pe 1952 2010 1710 | 2 1436 1485 - 
= 1951 1340 1020 5 2443 2593 
1950 443 560 | 10 3109 3326 
; 1949 1520 780 | 25 3951 4253 a 
1948 1420 2260 
ial 1947 753 860 | 
1946 1430 1360 
1945 1260 800 
1943 887 720 | 
1942 3280 3260 
1941 1240 1370 
1940 3410 2070 
4 1939 657 1000 
2 1938 2700 3390 
Mean 1566 1629 
8, 906 997 
7 1937 1270 480 2 1109 1160 ‘ 
4 1936 1860 1960 5 1747 1925 
1935 732 760 10 2170 2432 
iba 1934 2630 2170 25 2703 3073 
a 1933 1480 1480 
4 1932 1000 970 i 
1931 1320 770 
ig 1930 615 870 q 
ie 1929 1040 710 
1928 1850 1790 
1927 857 1820 
1926 1230 2740 
1925 692 440 
1924 556 1020 
1923 750 890 
mA Mean 1192 1258 q 
S. 0. 574 689 
Mode 903 912 
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not be perfect because of inaccuracies in rainfall-runoff relations and unit 
hydrographs and the probably greater inaccuracies involved in estimating the 
true volume and the true time and geographical distribution of the rainfall 
from the network data. In the study described by the writers the rainfall net- 
works were not generally located in the problem basins, because one of the 
primary objectives was to determine if the rainfall data from one network 
could be used to synthesize peak discharges for several basins so that the 
costs of applying the procedure on a per-basin basis would be thus reduced. 
The only limitation on the use of the rainfall data from a given network to 
synthesize the peak discharges for a particular basin is that both the network 
and the basin must have similar precipitation characteristics. Specifically, 
they must have approximately the same rainfall-intensity frequencies and 
their intense rainfalls must occur at about the same time of year. It certain- 
ly is not required that rainfalls of relative magnitudes over the network and 
over the basin occur on the same dates. 

The reasoning of the writers was that over a long period the basin would 
have a rainfall record similar to that of the network provided both were ina 
region of similar precipitation characteristics. This is certainly a better 
basis for obtaining flood-frequency data than the regional methods mentioned 
in the discussion. Users of regional methods have to assume that the hydro- 
logic and geologic characteristics, as well as the precipitation characteristics, 
of the problem basin are the same as those of the basins used in developing 
the regional method. Of these, the precipitation characteristics are the most 
conservative and may therefore be interpolated with more confidence than the 
hydrologic and geologic characteristics for an ungaged basin. 

The writers’ procedure does make use of causally related data for the 
problem kasin, namely, for adjusting or developing a rainfall-runoff relation 
and for constructing a unit hydrograph. If any method can be labeled as 
abandoning the use of hydrologic data causally related, it is the regional 
method of the usual type. Generally, the only independent parameter is the 
size of the drainage area. Even if the problem basin has a 5-year record of 
rainfall and streamflow measurements, the data are usually ignored in region- 
al analysis. 

Another advantage of the procedure described by the writers is that it can 
be adapted to the latest hydrologic characteristics of the basin and might, in 
some cases, yield more useful results than even the record of observed 
streamflow. For example, a small basin with 35 years of streamflow 
measurements has gradually undergone appreciable changes in land use during 
the first 30 years of record, If the runoff characteristics were seriously af- 
fected, the streamflow for the earlier years would not be comparable to that 
for the last 5 years. A frequency analysis of this record would have little 
significance in design. Neither would the application of the usual regional 
method. However, by adapting a rainfall-runoff relation and constructing a 
unit hydrograph on the basis of the last 5 years of streamflow data, the 
authors’ procedure could reasonably be expected to yield more reliable results 
than could be obtained from either the streamflow record or the regional 
method. 

The writers admit that there is a great deal of work involved in the appli- 
cation of the procedure, and it would be uneconomical to use it in the design of 
isolated low-cost structures. However, only about two man-months are re- 
quired for an experienced person to make the storm-rainfall survey (50-year 
record), establish the time distribution of the rainfall, adapt a rainfall-runoff 
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relation, construct a unit hydrograph, convert the rainfall into peak discharge, #1) 
and make the frequency analysis. The storm survey and time distribution of q 
the storm-rainfall data require about one man-month, so there is a 50 per 
cent reduction in time required if these two items have already been es- 
tablished for an applicable network. There are certainly many structures for ES 

° which the expenditure of one or two man-months in obtaining more reliable 
design data would be justified. 

In view of the comments in the discussion, it may appear strange that the 
writers look upon their procedure as an aid to, as well as a substitute for, 

8 regional or generalized methods. The development of regional methods is 
often hampered by the relatively short streamflow records available, especial- 
ly for small basins. The writers’ procedure would provide 50-year or longer 
records as a basis for such methods and thus extend their range of reliable 
application to include design values for longer return periods than is usually 
possible. 
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FLOW CHARACTERISTICS ON THE OGEE SPILLWAY# 


Discussion by Donald P. Thayer 
Closure by Robert B. Jansen 


DONALD P, THAYER, ! M. ASCE.—The author has done a service to the 
profession in bringing up and proposing a solution for one of the most signifi- 
cant problems confronting the hydraulic designer; i.e., the dissipation of 
energy by friction on an overflow spillway or steep chute. 

The writer is in general agreement with the author as he writes, “It is 
reasonable to assume that the energy loss on the spillway face is determined 
by a simple natural law.” It is felt that the most logical approach would con- 
sist of making a theoretical analysis, taking into account the major factors af- 
fecting the phenomenon. Such an analysis should be of immediate use for pre- 
liminary design and subsequently should be a convenient framework within 
which to correlate experimental measurements. From such correlation there 
should result emprical coefficients bringing the formula into close agreement 
with experience. By “experimental measurements”, the implication intended 
is measurement of prototype behavior; the writer’s observation has been that, 
for quantitative evaluation of friction losses, conventional models, based upon 


Froudian similitude, offer little help. More will be said of this later inthis 
discussion. 


Theoretical Analysis 


A theoretical analysis of the friction loss on an overflow spillway will be 


attempted. The well-known Darcy formula will be used for evaluation of the 
friction head, 


(1) 
where: f = the friction coefficient 


L = the length considered 
D = diameter of the pipe 
v2/2g = velocity head 


The choice of this formula is based upon two considerations: (a) it is di- 
mensionally homogeneous, and (b) it appears to be thoroughly documented by 
experimental data. It should be noted that this formula is identical, mathe- 
matically, with the much earlier Chezy formula; the respective constants ap- 
pearing in a different form and differing by a constant ratio. For the purpose 
a. Proc. Paper 1452, December, 1957, by Robert B. Jansen. 


1. Asst. Div. Engr., Div. of Design and Constr., California Dept. of Water 
Resources, Sacramento, Calif. 
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of this analysis, however, this formula has two shortcomings: (a) it has been 
derived for and its use has been limited, almost exclusively, to closed circular 
conduits, and (b) its friction coefficient f is a function, not of the roughness 
alone, but of €/D, the relative roughness, where € = the “absolute roughness”, 
ideally a linear measure of the height of the surface irregularities, but practi- 
cally an experimentally determined quantity. This formula is often applied to 
non-circular closed conduits by letting D = 4R, where R is the hydraulic 
radius, and may, with perhaps somewhat less validity be applied to open 
conduits, (1) 

The channels considered in this discussion will generally be very wide 
compared to their depth so that without significant error we may say that the 
depth of flow, perpendicular to the bottom, is equal to the hydraulic radius, R. 
With this in mind, the Darcy formula may be written 


he ve (2) 
J 2g 
If the hydraulic slope is defined as s = hg/L, then the formula becomes 
(3) 


In Fig. 1 consider an elemental prism of the flowing water of weight W, 
shown cross-hatched, moving down the spillway slope, A-B, from position M 
to position N. Let the line C-D be the trace of the centroid of the mass as it 
moves. Measured at this centroidal line, the hydrostatic pressure potential 
energy of the elemental prism of water is zero and its total energy may be 
written as the sum of its potential and kinetic energies, 


Wye 


E= De + Wy (4) 


where: y = distance of the centroid above some datum, as yet undefined, 
with positive direction upward. 


As the mass moves down the slope the change in energy may be obtained by 


differentiation with respect to L. be 
Wy gv + 


For a frictionless condition, the right-hand member of Eq. (5) would be 
equal to zero; with friction it is equal to minus the hydraulic friction slope 
times the weight of the mass, 


ge = -Ws 6 
aL Ws (6) 
From geometry of Fig. 1, i 

dL = edy/sina (7) 

where: o = angle of inclination of line C-D, at the point considered, from - 

the horizontal. 

From continuity of flow, t 
R = a/v (8) 


where: q = specific discharge passing down the spillway. 


one 
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FIG.1- OVERFLOW SPILLWAY 


Inserting relationships, Eq. (3), (6), (7), and (8) in Eq. (5) and multiplying out 
dL, gives 


Dividing by W and rearranging 


vdv fv3 
Let 
fv3 
B3 = sine (11) 
Then 
5 f 2/3 1712/3 
(12) 
Let 
ft 2/3 
sin (13) 
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v2 (14) 


Note also, that vo2/g = De, the depth at critical velocity, and substitute this 
relationship, together with the first differential of Eq. (14) for the left-hand 
member of Eq. (10) to obtain 


(15) 


Integrating 


The denominator under the integral sign has the factors 


[B - 2].[B+ (1 +4 v3) .[B+ (1-1 V3) 


These factors become the denominators of partial fractions formed in the 
usual manner. These partial fractions upon integration and the combination 
of conjugate complex quantities give as a value of the left-hand member of 


Eq. (16) 
- 2V3 tan-l #| = (3B) 


Bot Btl (17) 
Then, using this notation, Eq. (16) may be written 


F De 


Eq. (14) may be written 


where: hy = velocity head at any supercritical velocity 
he = velocity head at critical velocity 


nou 


It should be noted that the function (B), Eq. (17), is inherently negative, vary- 
ing from -2.4179 when B = 0, to negative infinity when B = 2. 

Eq. (18) and (19) then give dimensionless relationships between y, measured 
in terms of De, and hy, measured in terms of he. These equations have been 
evaluated for several different values of F and the results plotted in Fig. 3. 

In all cases the value of C in Eq. (18) has been so chosen that the curve passes 
through the point (hy/he = 1, y/De = +0.5) thus placing the datum as shown in 
Fig. 1, at the elevation of the crest where critical velocity is attained. It is 
recognized that in many cases critical velocity and depth are not attained at 
the fixed crest. In this case it is a simple matter to compute a critical depth 
at which the given specific discharge would pass over the given crest or, more 
precisely, to compute a fictitious crest elevation over which the given specific 
discharge from the given pool elevation would pass at critical velocity. It is 
believed that in many cases the former method is sufficiently accurate. In 


this connection the following relationship may be useful; the equivalent critical 
depth, 


Dt = H c2/3 
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FIG.2-DIAGRAM RELATING no, R, f, x, AND F 


where: H = head from the given crest to stagnation pool level 
C = coefficient of crest, as inq = C H3/2 


Fig. 3 has been drawn to give a rapid graphical method of obtaining the re- 


lationship between Manning n, the Colebrook function for f in the Darcy formu- 


la, and F as used in this discussion. The use of this diagram should be self- 
evident, except for the angle, a, to be used. Theoretically this is the angle 
the median line of the overflowing nappe makes with the horizontal; actually, 
at the lower part of the spillway, where the greatest part of the friction loss 
occurs, this line is nearly parallel to the spillway and the latter angle should 
be satisfactory. The depth, R, varies continuously throughout the spillway so 
that a value approximating the mean in the region where the greater part of 
the friction loss occurs would probably be the best to use. 


Comparisons and Application 


For purpose of comparison, three of the points representing velocity 
measurements at Madden Dam as reported by Mr. Randolph(2) have been plot- 
ted on Fig..3. The correlation is not particularly good, although the trend is 
in the right direction. It must be pointed out that the plotting may be inaccu- 
rate since the values were scaled from Fig. 17 in Mr. Randolph’s paper, the 
writer being unable to find numerical values. Two points from the model 
study of the Madden Dam, also scaled from this same Fig. 17, are also plot- 
ted. These points fall on a much smaller value of F than any of those for the 
prototype, although the roughness of the model and prototype, as reported by 
Mr. Randolph were in approximately the correct Froudian relationship 
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(assuming, of course, that the Manning formula is applicable to both). This 
emphasizes the statement made earlier regarding the inapplicability of model 
results where friction is concerned. 
Curves drawn according to the author’s Eq. (2) are also plotted on Fig. 3; 
one with K = 0.10 and the other with K = 0.125. These curves cut sharply 
° across the F curves indicating considerably discrepancy. It might be pointed 
out that the author’s curves are not asymptotic to a vertical line, indicating 
no terminal velocity. For practical purposes, with small relative values of y, 
this should not be considered a serious defect; the point is more theoretical 
° than practical. 
The analysis made in this discussion disregards two factors which might 
have a considerable affect on the spillway velocities: 


(a) The friction of the atmospheric air on the overflowing sheet of water. 
(b) The effect of entrained air on the friction coefficient. 


As the author points out, engineers are not agreed on this point and there are 
conflicting data. Viparelli(3) states that for calculation of velocity the value 
of n appropriate for the channel surface should be used. He cites earlier ex- 
periments by Mongiardini with flow of an air-water mixture through a closed 
pipe indicating that for conditions of complete turbulence (Reynolds’ number 
greater than 50,000) the friction coefficient is the same as for water alone. 

It should be noted that these experiments used a rather small pipe, 20.8 mm 


It is believed that more measurements of prototype behavior are needed to 
clarify this point. 


‘ diameter. On the other hand, Straub and Lamb(4) show that velocities in an 4 
i experimental chute computed from the “normal” value of n are considerably a 
less than those measured. 


Future Investigations 
si No other prototype measurements are available to the writer to compare 

Bi. with his analysis. It would be very helpful if anyone having such data would 

4 plot it on Fig. 3 or forward the data to the writer for such purpose. 


The lack of reported prototype measurements applicable to this subject is 
rather unfortunate, and it is hoped that more data comparable to that collected 
by Mr. Randolph will be collected and reported. It is encouraging to note that 
the north spillway bay of the Detroit Dam has been arranged to permit such 


measurements, and it is hoped that data from this source will be made availa- 
ble in the near future. 
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4. Straub, Lorenz G. and Lamb, Owen P., Experimental Studies of Air En- 


trainment in Open Channel Flow, Proc. Minnesota International Hydraulic 
Convention, Sept. 1-4, 1953, p. 435. 


ROBERT B. JANSEN,! A. M. ASCE.—The subject of head loss on the spill- 
way face is deserving of more attention. Further experimentation and analy- 
sis should be directed toward determination of the mechanics of this process 
of energy dissipation. The rate of energy variation surely must be determined 
by the relationship of a relatively small number of factors. In the case of an 
ogee spillway, or any straight spillway with constant slope and roughness, the 
pattern of energy change obviously is established primarily by conditions at 
the crest. 

For a given spillway slope and roughness, the head loss should be a function 
of the distance of fall and the rate of discharge. It was on this basis that Eq. 
(3) was suggested. Application of this formula, with K equal to 0.125, gives 
results fairly close to the prototype observations on the spillway of the Madden 
Dam. Since the face of this spillway was reported to be rougher than usual, 
it was estimated earlier that the value of K for a typical concrete spillway 
might be closer to 0.10. Recent developments have shown this assumption to 
be questionable. 

The newly published work of Bradley and Peterka(1) contributes to an 
understanding of this phenomenon. Fig. 15 of their series of valuable papers 
shows the results of prototype observations on the spillways of Shasta and 
Grand Coulee Dams. Also presented are curves, drawn from the considerable 
experience of the authors, which enable the determination of velocity and depth 
of flow “with a fair degree of accuracy, for chutes having slopes of 0.8:1 or 
steeper, where computation is a difficult and arduous procedure.” The proto- 
type tests cover a range of crest heads from 1.6 feet to 13.6 feet with height 
of fall as great as 220 feet. These important measurements provide oppor- 
tunity for comparison with the Madden Dam prototype data, and for further 
analysis of basic relationships. 

The Shasta and Grand Coulee spillways have a face slope of 0.8:1, as com- 
pared with the Madden Dam inclination of 0.75:1. All three have drum-gated 
crests; Shasta with three gates 110 feet by 28 feet, Grand Coulée with eleven 
gates 135 feet by 28 feet, and Madden with four gates 100 feet by 18 feet. 

Since the reported prototype measurements on these three spillways were 
limited to crest heads not greater than 13.6 feet, there is need for more infor- 
mation on flow at higher heads. No other prototype tests are known to the 
writer. However, flow at high crest heads is in a range where, at least for 
relatively low heights of fall, model and prototype results should be in fair 
agreement. That is, discrepancies attributable to roughness and air en- 
trainment should not enter seriously. It is considered acceptable, therefore, 
to use the results of model experimentation for comparison in this range. For 
this purpose, reference is made to two works of the Waterways Experiment 
Station, Corps of Engineers, Vicksburg, Mississippi. These are Bulletin No. 
37(2) and the model studies for Center Hill Dam.(3 Energy losses for large 


crest heads can be calculated by using the nappe profiles presented in these 
papers. 


1, Supervising Hydr. Engr., Dept. of Water Resources, State of California, 


Sacramento, Calif. 
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The data from the above various prototype and model tests have been ana- 
lyzed by the writer in an attempt to establish a generally reliable formula for 
computing the head loss on the face of an ogee spillway. Logarithmic graphs 
were used to establish exponents. It was discovered that the head loss varied 
directly as about the 4/3 power of the height of fall, and inversely as the 1/3 
power of the critical depth at the crest. (These are representative values, 
with emphasis placed on the apparently more consistent Shasta Dam test data. 
The exponents derived from the Maddem Dam prototype results were 3/2 and 
1/2, as reflected in Eq. (3)). These findings lead to the following relation: 

hy =J 4 4/3 


I (10) 


wherein J is a coefficient, Z is the height of fall from the reservoir level and 


is equal to vy sum of P and H, and D, is the critical depth at the crest and is 
C2/3H 
equal to "gi/3 ° 


A value of J = 0.15 yields results which agree well with the prototype and 
model records now available. As a matter of interest, a value of K between 
0.125 and 0.150 in Eq. (3) would give computed losses in good agreement with 
these same test data. It appears therefore that the suggested reduction of K 
from 0.125 to 0.100 because of the unusual roughness of the Madden Dam spill- 
way is not justified. A more proper conclusion from the Madden Dam analy- 
sis only, in view of present knowledge, would involve a coefficient of 0.125 in- 
stead of 0.100 in Eqs. (2) and (3). 

Although Eq. (10) is supported by only a meager experimental foundation, 
it has merits which should be noted. It is simple, easy to apply, and gives 
apparently dependable results. Use of the term D, includes the effects of both 
C and H, which were treated individually in Eq. (3), On the basis of the spill- 
way observations reported to date, Eq. (10) appears to be in a more acceptable 
form than Eq. (3). It is probably as consistent a criterion as can be offered 
now. 

Mention should be made of the influences of the roughness and the slope of 
the spillway face. The former should not be cause for great concern. In the 
construction of a modern concrete spillway, enough attention is usually given 
to concrete quality and formwork to ensure a uniform roughness. The effect 
of the slope of the spillway face cannot be determined on the basis of the three 
prototypes which have been studied. There is not enough difference in their 
face angles. However, most ogee spillways will fall into a range not much 
greater than the one in reference (0.75:1 to 0.80:1). To extrapolate the known 
data into the range of flatter slopes, it might be assumed that the energy loss 
varies directly with the length of travel along the spillway face, or inversely 
as the sine of the angle between the face and the horizontal. If the term J in 


Eq. (10) is assumed equivalent to <i, then the head loss is given by the 
following expression: 


ny = Jj 
sink 173 (11) 


in which j is a coefficient and a is the angle of inclination of the spillway face 
with the horizontal. 


Eq. (11) is only speculative. The subject of loss on slopes flatter than the 
usual gravity-dam slopes is beyond the scope of this paper. 
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It should be noted too that the prototype measurements on which this study 
is based are limited to values of Z generally less than 200 feet. However, 
this does not necessarily place a limitation on the applicability of the results. 
In using Eq. (10), the calculations of energy loss should be performed pro- 
gressively at intervals down the spillway face until the terminal velocity is 
reached. This is the point at which the increment in energy loss becomes 
equal to the vertical increment of fall. Since the velocity head is assumed 
constant below this level, the incremental head loss beyond the arrival at 
terminal velocity is taken as equivalent to the increment of fall. This means 
that Eq. (10) need be applied only in the range where the rate of energy loss 
is varying; i.e., only down to the point where the terminal velocity has been 
attained. 

Mr. Thayer has made a remarkable contribution to the literature in his 
theoretical analysis based on the Darcy formula. Essentially, it appears that 
he has independently accomplished the same result as achieved by Bresse in 
his method.(4) Bresse performed his integration assuming a very wide rec- 
tangular channel and adopting the Chezy formula. Mr. Thayer has used the 
same assumption, and has employed a formula which, as he states, is mathe- 
matically identical with the Chezy relation. 

Much has been written about the limitations of the Darcy formula as applied 
in the way that Mr. Thayer suggests. Admittedly arbitrary premises are in- 
volved. Open channels are usually rougher and larger than pipes, with high 
values of the Reynolds number, corresponding practically to complete turbu- 
lence. In this case, the Darcy coefficient “f” would be practically independent 
of the Reynolds criterion by which it is commonly determined. The belief has 
prevailed that some other index, such as the Froude number, should be em- 
ployed for flow with a free water surface. : 

It would be hard to imagine more dissimilar conduits than an ogee spillway 
and a small pipe such as that for which the Darcy formula was developed. Any 
attempt to use the Darcy relation for analysis of a steep spillway would be 
accompanied understandably by serious doubts. However, comparison of the 
Thayer results with data from other sources reveals a relatively close 
agreement for flow with low crest heads. Discrepancies appear at higher 
heads. 

Although Mr. Thayer has performed an analysis closely paralleling the 
work of Bresse and has suggested a criterion basically the same as the method 
of Bresse, the new approach has advantages over the earlier work. The 
Thayer graphs are presented in such a convenient form as to make compu- 
tation by this means probably more attractive than by the Bresse tables. Mr. 
Thayer, in his integration of incremental energy losses defined by the Darcy 
formula, makes unnecessary much of the labor that would be involved in the 
common step method using the same formula. This is especially true if the 
energy relation at only a single point on the spillway face is concerned. If the 
cumulative loss is to be determined at several points, there is less advantage. 
For practical purposes, the result is the same whether integration is ac- 
complished by the calculus or by small finite steps. 

The two graphs of Mr. Thayer provide a rapid means of estimating the 
energy loss on the spillway face. In using his curves, it is necessary first to 
estimate the dimension R and to enter the first graph with this trial value to 
find the factor F. Mr. Thayer states that, “the depth, R, varies continuously 
throughout the spillway so that a value approximating the mean in the region 
where the greater part of the friction loss occurs would probably be the best 
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to use”. It would first appear that this assumption would inject a considerable 
measure of approximation into the Thayer method. However, the writer has 
found that the suggestion has merit. Also, the extent of trial involved in se- 
lecting the value F is not an important handicap. The convenience of the first 
graph makes the trial a simple matter. 

When the sheet of flow is thin, the graphical selection of the proper value 
of the Thayer factor F becomes rather difficult. In this range, F varies ap- 
preciably with small changes in depth. The term F can be expected to vary 
as much as 100 per cent on a typical ogee spillway. This is well demonstrat- 
ed by the first graph of Mr. Thayer. 

For large heads on the spillway crest, in general the second Thayer graph 
is too small in scale to yield sufficient accuracy. This is not an unusual 
disadvantage of a graph prepared for publication. In practice, the plot could 
be made to a larger scale, or the basic data could be placed in tabular form. 
However, it is not certain that the method of Mr. Thayer will give acceptable 
results when the ratio y/D,. is small. More experimental verification is need- 
ed in this region. 

Table 1 is a compilation of data from the several sources which have been 
mentioned. The tabulation shows that the energy losses computed from the 
actual observations on the Shasta spillway prototype are generally illogical 
for higher magnitudes of the dimension Z. For example, with a crest head of 
1.6 feet the listed energy losses from actual measurements on the prototype 
are 40, 59, 77, 92, 103, and 108 feet corresponding to heights of 70, 90, 110, 
130, 150, and 170 feet, respectively. The increments of energy loss are 19, 
18, 15, 11, and 5 feet, in the same order. Logic dictates that, since the ve- 
locity will continuously increase until the terminal velocity is attained, there 
should be no decrease in the increments of energy loss in proceeding down 
the spillway face. The prototype records presented by Bradley and Peterka 
for Shasta and Grand Coulee spillways manifest a general defect in this re- 
spect. In order to eliminate this inconsistency, the writer had adjusted the 
computed losses for the Shasta spillway by accepting the value of hy, corre- 
sponding to the lowest recorded value of Z for each crest head, on the as- 
sumption that the accuracy of measurement of flows on the prototype would 
decrease as the height of fall increased. In making the adjustment, the rate 
of energy loss was not allowed to decrease. In the case of the Grand Coulee 
spillway observations, since the initial value of hy, for a crest head of five 
feet appears obviously too low, the first listed energy loss has been made to 
conform to the Shasta prototype data. The other adjusted losses shown for the 
5-foot crest head fit the loss differentials of the original Grand Coulee 
measurements, except that again no decrease in the loss rate has been permit- 
ted. For a crest head of 13.6 feet, the measured energy loss for the lowest 
value of Z has been accepted, although it may be slightly high, and an incre- 
mental energy loss of ten feet has been adopted from examination of the 
pattern of the Shasta columns. 

The tabulated energy losses on the Madden Dam spillway were taken direct- 
ly from Fig. 17 of the paper by Mr. Randolph. A general nappe profile (Fig. 
8) in Bulletin No. 37 of the Waterways Experiment Station has provided an 
estimate of losses at crest heads of 20 and 30 feet. Similar data from the 
model studies of Center Hill Dam, performed at the same laboratory, give 
approximate losses for flow with 44 feet of head on the crest. These measure- 
ments were taken from the original model, which ‘iad a face slope of 0.61:1. 
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Presented also in Table 1 are energy losses determined on the basis of 
curves drawn in Fig. 15 of the Bradley-Peterka series. This graph shows 
that the authors themselves have adjusted for the obvious defects in the proto- 
type observations on the Shasta and Grand Coulee spillways. 

The last two columns in Table 1 are based on the Thayer method with 
Darcy coefficients corresponding to Manning n = 0.014, and on the Jansen Eq. 
(10) with the factor J equal to 0.15. It is believed that Table 1 demonstrates 
that Eq. (10) offered by the writer gives a fairly reliable estimate of the ener- 
gy loss on the face of an ogee spillway, on the basis of comparison with the 
prototype and model observations which have been discussed. As has been 
stated, the Thayer method when applied to large crest heads is at variance 
with the other listed references. However, information presently on record 
does not permit a definite conclusion in this regard. More observations on 
spillway prototypes are sorely needed. 

One avenue of investigation which can be explored, at least as an approxi- 
mate means of correlation, is application of the law of hydraulic model si- 
militude. For example, disregarding lack of similarity in roughness and in 
air entrainment, flow on a spillway 200 feet high with a crest head of 20 feet 
should involve approximately twice the energy loss as on a spillway 100 feet 
high with a ten-foot head on the crest. 

The superposition of the data of the writer on the second Thayer graph 
would have resulted in a more logical curve if attainment of terminal velocity 
had been taken into account. As has been stated earlier, the equation defining 
the energy loss should be applied only down to the point where the terminal 
velocity has been reached. Eq. (10) is recommended for this purpose, until a 
more comprehensive formula can be written. It is recognized that this 
equation is the result of a mere scratching of the surface of this important 
subject. Certainly a more satisfactory criterion ultimately will be discovered. 
In the meantime, the procedures herein described may serve as a guide. 

The profession is fortunate to have received the outstanding works of 
Messrs. Randolph, Bradley and Peterka. Their important contributions to the 
literature have laid the foundation for a broader understanding of spillway 
hydraulics. In closing this discussion the writer wishes to record appreci- 
ation of the notable treatise by Mr. Thayer. The time which he has given to 
the clarification of the subject matter herein has greatly enriched the paper. 
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" TABLE 1. - ENERGY LOSS ON SPILLWAY FACE 
Loss in Energy Head, hy in Feet { 
a 
¥ 8 gy 
i] 
in: in: : in 2 
Ft.: Fte: =o: Ft. 
1.2 133 3.0 0.79 131 lilo 119 110 
1.6 70 3.0 1.05 ho k0 ho 53 13 
1.6 90 3.0 1.05 59 59 s8 472 60 | 
1.6 110 3.0 1.05 77 78 77 92 77 5 
1.6 130 3.0 1.05 92 97 97 112 97 " 
1.6 150 3.0 1.05 103 116 117 (132 117 S 
1.6 170 3.0 1.05 108 135 137 152 137 2 
2.0 13h 3.0 1.31 108 97 112 9b 
le2 136 3.2 2.87 70 78 85 73 
80 3.1 3.07 33 33 30 33 «35 
110 3.1 3.07 53 54 53 8 55 
be6 130 3.1 3.07 65 68 67 77 68 ey 
150 3-1 3.07 73 82 86 82 
4.6 170 3.1 3.07 78 % 105 115 97 
5.0 100 3.1 3.34 19 bb kl b6 
5.0 120 3.1 3.3h 56 6h 59 
140 3.1 3.34 48 73 72 82 73 q 
Se0 160 3-1 3.34 64 89 90 101 87 4 
180 3.1 3.34 77 105 110 120 102 
5.0 200 3.1 3.3h 89 121 130 139 117 
5-0 220 3.1 3.3h 98 137 150 159 133 
6.0 138 35 59 65 62 66 
766 100 3.3 5.30 35 35 32.25 
7.6 110 3.3 530 WM 37 4s 
705 130 3.3 5.30 2 53 
766 150 3.3 5.30 60 65 63 55 69 
7.56 170 3.3 5,30 65 77 79 70 81 
7.78 1h0 3.6 5.75 48 ss 
10.6 130 3.5 7.68 3 43 37 29 SO 
13.6 130 305 9.86 499 32 19 16 
13.6 150 3.5 9.86 Sh 59 39 27 56 
13.6 170 3.5 9.86 49 49 3h 66 
20 60 15.85 15 3 #1 
30 90 heO 23.78 23 19 #1 21 a9 
120 heO 34.87 32 7 
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DISCHARGE CHARACTERISTICS OF RECTANGULAR THIN-PLATE WEIRS@ 
Discussion by M. R. Carstens 


M. R. CARSTENS,! A. M. ASCE.—The authors are to be commended for 
the study of rectangular-notch weirs. 

The authors have attempted to eliminate the effects of viscosity and surface 
tension from experimentally determined values of the coefficient discharge of 
a thin-plate weir by the simple expedient of an effective head and width. The 
successful elimination of these two effects should result in the coefficient of 
discharge of potential (irrotational) flow without surface tension. Thus, the 
equations for the coefficient of discharge for the full-width weir, Eqs. (25), 
(26), (27), (28), and (29), should be identical for these five series of tests. The 


small differences in these equations could be attributed to the following 
causes: 


(a) The use of the effective head does not completely eliminate the effect of 
viscosity and surface tension and the differences in the equations repre- 
sent differences between the five series of experiments in this regard; 

(b) The geometry in each of the five experiments was sufficiently different 
to account for the small differences; and 

(c) The differences are merely an indication of the order of accuracy with 
which the coefficient of discharge can be experimentally determined. 


Since the maximum deviation of the data from each experiment (Figs. 3, 4, 
5, and 9) is greater than the deviations between the equations, the writer be- 
lieves that (c), above, is the explanation for the small differences in Eqs. (25) 
- (29), inclusive. Therefore, the writer would have combined all the data (re- 
duced to potential flow) shown in Figs. 3, 4, 5, and 9. 

Despite the tenuous derivation of Eq. (15) the solution for the coefficient of 
discharge is amazingly good. Eq. (15) for a full-width weir is 


% 
2 2 2 
q = Ce | (15a) 


This equation can be placed in the usual form of the discharge equation by 
factoring the quantity, h3/2, from the parentheses. 


3 


2gh 29 


. Proc. Paper 1453, December, 1957, by C. E. Kindsvater and R. W. Carter. 
Prof., Civil Eng., Georgia Institute of Technology, Atlanta, Ga. 
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The quantity, V12/2gh, can be transformed as follows: 


q = V, (h + P) 
and 


q = 2 Cy N29 7% 
3 


from which 


Substituting into Eq. (15b) 


% 2 3h 
9 h/P +t 9 h/P +1 (15c) 


If the theoretical values of C, for a two-dimensional orifice are utilized in 
Eq. (15c), the equation can be solved by successive approximation for Cg as 
a function of h/P. In the following table, the theoretical values of Cy are tabu- 
lated for comparison purposes with the best available experimentally de- 
termined values, that is, Eqs. (25), (26), (27), (28), and (29). In addition, the 
experimentally determined values of Kandaswamy and Rouse are tabulated in 
order to compare the theoretical solution with an experimental solution ob- 
tained at larger values of h/P. 

The theoretical solution for the coefficient of discharge is within a few per 
cent of the various experimental solutions throughout the entire range of the 


weir, that is, 0< h/P < 10. The writer cannot agree with the authors’ 
statement, 


“The result is a complex equation which must be modified to make it 
practicable and adjusted to make it agree with a particular set of ex- 
perimental data. The fact that a large number of different empirical 
formulas have been evolved by this procedure is evidence that the 
“theoretical” equation is inadequate. The obvious reason for the 
inadequacy is the fact that the basic assumptions are inaccurate.” 


Rather it appears to the writer, that too many people write empirical 
equations. 
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AIR BINDING IN LARGE PIPELINES FLOWING UNDER VACUUM@ 


Closure by R. T. Richards 


| Ae RICHARDS, ! A. M. ASCE.~—The author agrees with Messrs. 
Templeton and Stelson that air evacuation equipment may be required at verti- 
cal drops. This is particularly true of low velocity flow. The equations they 
have submitted give an interesting comparison of the relative head losses at 

a vertical and sloping drop. However, the author’s observations on large di- 
ameter pipelines indicate that at vertical drops preceeded by relatively short 
radius bends and at water velocities of 6 to 8 feet per second, turbulence ap- 
pears to reduce air collection to a point where it causes no serious head loss. 
Such bends may be found, for example, at the majority of steam power plant 
main condenser discharges. While a few of these condensers become air 
bound unless continuously evacuated, the majority are self priming. Hy- : 
draulic efficiency here is of great importance. A sloping pipe from the same 
condenser discharge would probably be air bound. ‘ 

Messrs. Templeton and Stelson close their discussion by appropriately 
commenting that air binding can be a serious problem in conduits flowing 
under pressure as well as under the vacuum conditions discussed in the paper. 
The author’s firm has recently investigated head loss problems in two large 
diameter sloping power tunnels flowing several miles under pressure. The 
collection of air along the tunnel roof appears to have bee’: a major cd «ribut- 
ing factor in the observed excess head losses. Air relief taps similar to those 
described in this paper have been provided on several large diameter sloping 
pressure pipes designed by the author; however, to date the test data from 
these pipes is not sufficient to clearly establish the effect of air binding on 
these particular installations. If continuous air removal is required it can 
very easily be accomplished by automatic air relief valves. Removal of air 
from a pipe flowing under vacuum is far more of a problem, especially at 
pipeline locations where electric power is not readily available for air 
evacuation equipment. 


. Proc. Paper 1454, December, 1957, bv R. T. Richards. 
Civ. Engr., Ebasco Services Inc., New York, N. Y. 
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TURBULENCE CHARACTERISTICS OF THE HYDRAULIC JUMP2@ 


Discussions by A. J. Peterka and J. N. Bradley, James M. Robertson, 
Donald R. F. Harleman, Philip G. Hubbard 


A. J. PETERKA,! M. ASCE and J. N. BRADLEY,2 M. ASCE.—The authors 
are to be complimented for accomplishing an investigation which is difficult 
both in concept and in performance. Their study represents the ultimate ap- 
proach; to explain the workings of the hydraulic jump by means of its internal 
mechanics. Although the study is academic in nature, there are certain as- 
pects of the results which help the practicing hydraulic engineer to rationalize 
the use of stilling basins shorter than the hydraulic jump iength. 

It is common practice to make the stilling basin floor (usually concrete) 
shorter than the true length of the hydraulic jump. Hydraulic model tests and 
prototype verification tests have shown that, within limits, this is a desirable 
procedure since little in the way of performance is lost and the cost of the 
stilling basin is reduced appreciably. In analyzing the action in the jump, the 
authors have shown why it is possible to reduce the apron length to less than 
the jump length. The drawings in Fig. 13 indicate that most of the losses oc- 
cur upstream from the end of the roller and that the downstream end of the 
jump is relatively ineffective in accomplishing hydraulic losses. However, 
the velocity distribution curves in Fig. 6 show that a redistribution of velocity 
occurs after the major part of the losses have occurred. For this reason, it 
is not possible to reduce the apron length as much as appears possible in Fig. 
13. For example, in Fig. 13, for Fy = 2, it might be concluded that the apron 
need extend only about 3dg downstream from the start of the jump. Fig. 6 for 
F, = 2, however, shows that at jump length 3d9 there is a major redistribution 
of velocity yet to be accomplished. If the apron were ended at this point, the 
high bottom velocity might scour an unprotected riverbed. An end sill on the 
apron would help to redistribute the bottom velocity. Even at the end of the 
jump, 4.5dg, as indicated on Fig. 4, a slight redistribution of velocity is still 
to be expected according to Fig. 6. 

In tests reported by the writers3 the jump was found to be imperfectly 
formed for a Froude number of 2, (dg = 2.1d,). The total energy loss was 
only about 7 per cent in a jump measured to be 4.5d9 long. Thus, the length 
of the jump was useful more to redistribute the velocities than to accomplish 


a. Proc. Paper 1528, February, 1958, by Hunter Rouse, T. T. Siao, and 
S. Nagaratnam. 

1. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 

2. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 

3. Bradley, J. N. and Peterka, A. J., “Hydraulic Design of Stilling Basins”, 
Journal of the Hydraulics Division, Vol. 83, No. HY5, October 1957. 
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losses. Nevertheless, the apron length could not be reduced appreciably from 
the 4.5d9 given as the length of jump. 

For a Froude number of 4, the writers found the hydraulic jump to have 
better form since dg = 4.5dy. The energy loss was about 40 per cent. The 
incoming jet showed some instability, tending to flutter from top to bottom of . 
the jump. Surface waves downstream from the jump were more of a problem 
than scour tendencies. It was found that the apron length could be reduced 
considerably, however, without changing the jump action. With the addition of 
chute blocks and an end sill, the apron length was reduced to about 3.6d9. For 
a Froude number of 4, the authors show in Figs. 13 and 6 that most of the 
losses have already occurred at 3.6d9, and only a slight redistribution of 
bottom velocity will occur downstream from 3.6d9. 

For a Froude number of 6, the writers found the hydraulic jump to have 
ideal form, dg = 8d 1, and the energy loss in the jump was about 55 per cent. 
The apron length was reduced to 4d9 when cute blocks and an end sill were 
added to the apron; Figs. 6 and 13 help to justify this length reduction. 

It is unfortunate that the authors’ tests for F; = 8 could not be completed 
since these data would probably have shown that the roller occupied more 
than half the jump length as stated in the conclusions. It was found in the 
writers’ tests that as the Froude number increased, the roller size increased 
until the roller length became the controlling factor in determining the jump 
length. The roller length approached the jump length for Fy = 8, dg = 10.7d}, 
and the energy loss in the jump was about 65 per cent. A discharge of 28 
second-feet in a flume 4 feet wide produced a jump which was very turbulent 
and rough and the jump itself surged because the bottom jet did not carry 
through at all times for the full length of the jump. Bottom jet action was ex- 
tremely erratic, sometimes penetrating the jump for the full length and at 
other times, for only a short distance. At times, the jet appeared to jerk 
horizontally and flutter vertically. 

Also, the large difference in conjugate depths resulted in the jump having 
a steep profile. Slugs of water intermittently rolled down the steep face and 
fell into the high velocity jet. A surface surge was set up each time a slug 
made contact with the jet. This was strictly a gravitational effect which an 
air model may not duplicate. Actually, flow down the face is common to all 
jumps, but it is usually steady. These findings probably help to explain the 
difficulty found by the authors in adjusting their air model to obtain good data 
for Fj = 8. The writers, therefore, suggest that perhaps the agreement be- 
tween the air model and water model at Fj = 8 was better than suspected by 
the authors. In the writers’ opinion, it would be impossible with presently 
known equipment to make accurate measurements in a water jump of F, = 8 
having appreciable physical size. In addition, if the air model boundary was 
set to the jump length, curve of Fig. 4, it is possible that the dog test length 
was too short. In the writers’.tests, the jump length was found to be consider- 7 
ably longer for the higher Froude numbers than shown in Fig. 4. It is not 
clear in the description of the water model tests whether the length ratio 
Lj/dg was actually measured in the water model or whether the Bakhmeteff- 
Matzke curve was accepted as the true length of jump. In the air model, the 7 
shape of the boundary (the jump profile in a water model) could greatly influ- 
ence the action within the jump. The jump length given in Fig. 4 for F; = 8 
is 4.2d9; the writers’ tests showed 6.ld9. The reason for the difference is 
believed to be due primarily to the greater thickness of incoming jet and width 
of test flume used in the writers’ tests. Excessive friction in too small 
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models can lead to erroneous results. The authors found it necessary to 
widen their air model to decrease the wall effects noted in the 1-foot duct. 
This should serve as a warning to other experimenters who are tempted to 
use very narrow test facilities when wider sections could be used. 


JAMES M. ROBERTSON,! M. ASCE.—The study reported in this paper is 
not only of great interest because of the light which it sheds on that most in- 
triguing hydraulic phenomenon—the hydraulic jump—but also because it repre- 
sents a thorough study of a separated flow which reattached. Many flow situ- 
ations, besides the jump, occur with a separated zone of flow and a roller- 
like action of the separated flow on the main stream. 

The suitability of air flow in the carefully formed diverging channel as a 
model of the hydraulic jump is justified by the authors on the basis of expedi- 
ency and verification of the mean flow pattern. With no intention of deprecat- 
ing the value and significance of the study, the writer will point out certain 
aspects of the flow model wherein it cannot have been similar to a jump. Two 
features of the flow model differed from those of the hydraulic jump; one of 
these was inherent in the type of model chosen. The presence of a curved 
rigid boundary rather than the upper free surface of the jump certainly caused 
some modification in the nature of the flow while the initial mean and turbu- 
lence velocity conditions for the model jump differed from those normally ex- 
pected in open channel flow. 

The conditions of flow at the initial station of most hydraulic jumps must 
be those of fully developed turbulent open channel flow. The probable vari- 
ations in mean velocity, turbulent shear stress, and turbulence intensity for 
this situation are sketched in Fig. A. In contrast, the initial conditions tested 
in the air flow analogy involved relatively thin boundary layers at the upper 
and lower surfaces with an intervening core of essentially shearless, turbu- 
lence—free flow. In the F, = 2 case for instance, the shear stress curves of 
Fig. 7 indicate that the boundary layer thickness along the bottom was at most 
one-quarter of the initial flow “depth” even at relative distances x/do of at 
least 2. Since turbulent boundary-layer behavior in adverse pressure gradi- 
ents is sensitive to initial flow conditions, the actual flow following this con- 
dition must have differed considerably from that in a proper hydraulic jump. 
In this respect the condition studied would appear to conform closely to a 
jump following a sluice gate—a situation still of considerable interest. 

Since the pressure increases as the flow proceeds in the x-direction (cf. 
Fig. 5) the writer would expect to note evidence of separation along the “bed” 
of the channel. That none is in evidence might appear to be due to the ex- 
istence of an abnormally thin boundary iayer along this surface. Further con- 
sideration indicates that this explanation does not appear to hold water. The 
likelihood of turbulent boundary-layer separation occurring has been shown 
(1) to depend upon the boundary-layer Reynolds number at the start of the 
pressure gradient and upon the rate of pressure rise. The paper does not 
give sufficient information to properly determine the initial conditions but it 
would appear that at the initial station the momentum thickness Reynolds 
number probably lay between 103 and 104. If so, it is expected (1) that the 
maximum mean-flow velocity could not be changed by a factor less than 0.76 
to 0.69 for the Fy = 2 test and 0.78 to 0.67 for the F, = 6 test. Actually the 


1. Prof. of Theoretical and Applied Mechanics, Univ. of Dlinois, Urbana, Ill. 
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FIG. A VELOCITY AND TURBULENCE 
VARIATION IN UNIFORN\ FLOW 


apparent reduction in maximum velocity is nearer to 0.50 and 0.12 in the two 
cases. It is not clear why separation did not occur along the “bed” at an x/do9 
location of about 1.2. 

The writer’s second reservation concerning the “air modt1” of the hy- 
draulic jump has to do with the upper or “free” surface of the jump and the 
roller. Since the streamline separating the ror from the main flow and the 
bed pressure distributions are well verified, if@ppears that the model is ade- 


quate. This mean flow similarity is certainly a necessary condition for simi- at 
larity, but it is essentially geometric in character. From a dynamic: stand- we 
point two other factors enter. The first of these is that in the air model a pie 


finite wall shear must have occurred due to the air motion past the rigid upper 
boundary. In the actual jump no such shear occurs, at least if one neglects 
air friction. The other divergence of the upjier surface from the true hy- 
draulic jump is less simply presented. In the truly-hydraulic jump the surface 
and its roller appear to be far from steady. Instantaneous views of the upper 
surface of the jump tend to give the appearance of a breaking wave falling 
backwards on the oncoming flow. It is felt that the fluctuating nature of the . 
free surface may have a significant effect on the flow at some distance below 
the surface even though none is apparent in the mean dow characteristics. 

The upper boundary, or free surface, differences between the air model 
and water prototype may be likened to the differences between the outer 
(distant from the wall) portions of boundary layer and pipe flows. Early 
boundary layer analyses assumed that the velocity profiles (and by impli- 
cation, the turbulence pattern) in the two cases were the same if the boundary 
layer thickness was equated to the pipe radivs. More recently careful studies 
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ys have indicated significant divergences, both in mean velocity and turbulence 

ae characteristics, in the region away from the wall. G. B. Schubauer(2) has 
shown that this difference in turbulence is completely explained through the 
phenomenon of intermittency which occurs only at the outer or free stream 
edge of turbulent flows. As sketched in Fig. B, the division between the turbu- 
lently flowing fluid and that outside the boundary layer, does not always occur 
at the same location, as in the case laminar boundary layers. Instead the 
interface, between turbulent and non-turbulent flows, fluctuates with time and 
is conceived of as waves or wrinkles which travel along the outer edge cf the 
layer. 

The fluctuating nature of the outer edge of a turbulent boundary layer 
differs vastly from the situation that obtains at the center of a pipe where the 
fluid is always turbulent. A velocity indicating instrument with fast time re- 
sponse, when placed in the outer boundary-layer region will see patches of 
turbulent and non-turbulent fluid pass by. This condition is termed inter- 
mittency and the fraction of time that the flow is turbulent is defined as the 
intermittency factor Y. Schubauer(2) showed that if we consider only the time 
when the flow is turbulent in the boundary layer the turbulence in the outer 
region agrees well with pipe measurements. It is probable that the differ- 
ences in mean velocity profile are also tied in with intermittency. 

Intermittency cannot be immediately applied to our hydraulic jump. Even 
in uniform open-channel flow, the free edge of our flow is subject to a strong 
limitation and the air-water interface must inhibit if not eliminate the oc- 
currence of this phenomenon. In the case of the hydraulic jump it would ap- 
pear that the unsteady roller may result in an intermittency-like behavior. 
Certainly the upper surface of the roller is rather intermittent in shape and 
position. The occurrence of air bubbles in this region should also lead to 
intermittency. Asa further bit of suggestive evidence it may be noted that in 
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simple boundary-layer flows the wave or wrinkle length of the free boundary 
is of the order of two boundary-layer thicknesses. This is of the same order 
of magnitude as the relative length of the surface roller in the jump. 

Although the writer questions whether the authors have adequately 
“modeled” the hydraulic jump, he agrees that some very suggestive infor- 
mation has been obtained pertinent to the flow behavior in the jump. Actual 
turbulence studies of a wet hydraulic jump will be needed to ascertain the ex- 
tent to which the model studied was appropriate. 

Aside from its hydraulic jump connotation, the study reported is an ex- 
tremely valuable piece of basic research. Very few turbulence measurements 
have been inade of separating and reattaching flows. And yet, these are quite 
common occurrences (in valves for instance). The only comparable study of 
which the writer is aware was also conducted at the Iowa Institute of Hydraulic 
Research by A. A. Kalinske(3) and co-workers. The sudden expansion and 
large angle diffusers studied in that project involved separation and subse- 
quent reattachment. The experimental techniques used were much cruder 
than in the present instances, but somewhat similar conclusions were reached. 
Thus Kalinske stated, “the maximum total turbulence energy is a small part 
of the total energy change taking place” and “the principal loss of energy oc- 
curs at the outside of the high velocity stream as it passes into the expansion 
due to the ve'y high shear developed.” 

The relatively rapid disappearance of the extra turbulence energy generat- 
ed in their expanding flow is noted by the authors. Kalinske’s sudden ex- 
pansion results indicate almost as rapid a decay. Since to the observer of 
the hydraulic jump it still appears quite “turbulent” at its downstream end, 
where Fig. 10 shows a low level of turbulence, one may question whether to 
believe his eyes. Another possible explanation may account for this anomaly. 
In other flow situations (decay of isotropic turbulence) the turbulence energy 
is dissipated through viscous action in the smaller sized eddies while a rela- 
tively small portion of it associated with the large scale eddies, is little ef- 
fected. In the case of the jump at the downstream station the turbulence level 
may be low but large eddies could still be quite prominent even though con- 
taining little energy. This possibility could be verified through spectrum 
measurements of the turbulence. 

Study of Figs. 6 and 7 should give some idea as to the manner in which a 
separated and then reattached flow redevelops into a stable pattern again. 
Following the roller the turbulence levels are seen to be high, but the mean 
velocity profiles appear to ignore the upper surface. Analysis of mean- 
velocity data from other studies of reattaching flows have suggested to the 
writer that a new boundary layer forms along the wall more or less inde- 
pendently of what has gone on before. It may be that the measurements in the 
jump case did not proceed far enough downstream to indicate this. One of the 
flow cases analyzed, with the aid of H. R. Fraser, is depicted in Fig. C. This 
involved the flow of a jet into a partly throttled pipe. (4) Although the pitot 
tube used was not calibrated for the reverse flows measured, a significant 
separation zone is apparent. The low values of the momentum thickness 6 
and displacement to momentum thickness ratio (6 */g ) found at x/D of about 
8 suggest the start of a new boundary layer. Although not at all accurate the 
wall shear stress variation found in the separated region is of the proper sign 
and approximate magnitude. In contrast the apparent behavior of the wall 
shear stress along the upper curved surface in the authors’ investigation 
seems incomprehensible. In the separated zone the flow direction near the 
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surface is reversed and yet Fig. 7 indicates a negative wall shear stress 
within and following the roller. 

To calculate the turbulence dissipation the authors utilized the isotropic 
relation (constant K across and along the flow). They did adjust the value in 
a bulk fashion for the several Froude number models studied. Certainly K 
must vary throughout the flow field since the turbulence is far from isotropic 
and the kind or degree of isotropy can hardly be uniform. P. S. Klebanoff(5) 
in his very thorough study of the simple boundary layer, measured a majority 
of the terms appearing in the turbulence dissipation relation. In places some 
of these differed by a factor of ten from the isotropic values. A transverse 
integral of the dissipation, based on plausible assumptions for the unmeasured 
terms, indicated over twice the dissipation calculated with the isotropic re- 
lation. It is hard to believe that the more involved pressure gradient and 
separation flow case in the jump model would come out as close. Therefore 
the dissipation relations of Fig. 11 must be considered as of suggestive rather 
than absolute value. 

Even within the limitations which it is felt should be placed on the results 
presented, the paper reports on an extremely valuable and informative study. 
The writer could only wish for more detailed information on the quantities 
measured in order that they could be studied independently of the hydraulic 
jump framework in which the paper is cast. 


REFERENCES 


Robertson, J. M., “Prediction of Boundary-Layer Separation,” Journal of 
the Aeronautical Sciences, Vol. 24, 1957, p. 631. 


2. Schubauer, G. B., “Turbulent Processes as Observed in Boundary Layer 
and Pipe,” Journal of Applied Physics, Vol. 25, 1954, p. 188. 
3. Kalinske, A. A., “Conversion of Kinetic to Potential Energy in Flow Ex- 
pansions,” Trans. ASCE, Vol. 111, 1946, p. 355. - 
4. Carrier, W. H., and Misener, W. S., “Mixing and Equilibrium in Pipe + 
Flow,” ASME Paper 53-S-44, May, 1953. 
5. 


Klebanoff, P. S., “Characteristics of Turbulence in a Boundary Layer with 
Zero Pressure Gradient,” NACA Report 1247, 1955. 


DONALD R. F. HARLEMAN,! A. M. ASCE.—The authors have clearly 
demonstrated yet another aspect of the universal nature of fluid mechanics by 
their detailed study of the hydraulic jump by means of an air model. 

As a supplement to an otherwise comprehensive treatment the writer would 
like to set down the explicit equation for the exact depth ratio of the hydraulic 
jump. This will permit an investigation of the small but systematic deviation 


of the experimental depth ratio, dg/d,, from that predicted by the elementary 
momentum analysis given by, 


(a) 
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This deviation is shown graphically in Fig. 4 where the author’s hydraulic t, 
experiments, together with those of Bakhmeteff and Matzke, are compared 4 


with Eq. (a). At a Froude number Fy = 6, for example, the measured depth 
ratio is five per cent lower than that given by the elementary theory. 
The analytical expressions for the hydraulic jump depth ratio taking ac- 

c count of (1) boundary shear, (2) the excess of mean momentum flux due to non- 
uniform velocity distributions at the beginning and end of the jump and (3) the 
momentum flux due to turbulence have also been presented by Van Driest.(1) 
The author’s paper, however, permits for the first time a quantitative evalu- 

F ation of these terms from direct measurements. 

The exact equation for the jump depth ratio d9/d, can be obtained by evalu- 
ating the indefinite integrals of Eq. (13) at the end of the jump where d = do 
and x Lj; Eq. (13) then becomes, 


(Two typographical errors in the author’s Eq. (13) have been corrected, i.e. 
the dy appearing in the last term should be dx and the sign of the last term 
should be minus, since the bottom shear acts on the fluid in the same direction 
as the downstream hydrostatic force.) 

Upon introducing the dimensionless coefficients 


(b) 


where By and Bo are the momentum flux correction terms for non-uniform 
velocity distribution and Ig is the downstream turbulence flux correction term 
(Ij has been assumed equal to zero). After dividing through by Y Eq. (b) may 
be re-written as follows: 


By introducing the continuity equation, Uy dy = U2 dy and defining the Froude 
. and Reynolds numbers as 


Eq. (c) may be written as, ig 


[(T+iNT-1) + S] (a) 
A - (Ar.+1.) 
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and 
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Eq. (d) should reduce to the form of the elementary momentum solution of 
the hydraulic jump if 8; = 82 = 1 and Iz = 0 and Rj = ~ (frictionless flow, 
S = 0). Under these conditions, 


F*= (T+!) (e) 


which is in agreement with Eq. (a). 

The dimensionless quantities B;, 82, Ig and S can be evaluated from the 
author’s experimental data for each of the Froude numbers tested (F, = 2, 4, 
6) and the corresponding values of the depth ratio J can be calculated from 
Eq. (d). 

Table I gives the quantities determined from Figs. 7 and 8 of the author’s 
paper. Each of the downstream quantities B9, Ig and S were determined at 
the end of the jump x/dg = 5. 


Table I. 


Momentum, Turbulence and Shear Parameters 


By Bo I, J 
eq. (dj eq.(e) 


2 1.01 1.02 201 0.16 2.31 2.37 
1.01 1.05 206 1.05 5.02 5.18 
6 1.01 1.02 015 3.60 7270 8.00 


Also given in Table I are the values of the depth ratio J = dg/d; from the 
exact analysis, Eq. (d), and the elementary analysis, Eq. (e). The exact depth 
ratio is seen to be systematically less than the elementary. 

A final verification of the author’s measurements and the general validity 
of Eq. (d) may be demonstrated by comparing the three predicted (exact) 
values of the depth ratio with Fig. 4 of the author’s paper. This comparison 
is shown in Fig. A and indicates that the predicted values (shown by the dashed 
line) account for the departure of the experimental depth ratios from the 
elementary momentum equation for both the hydraulic and air model tests. 

The bottom shear S is by far the most important corrective term in the 
exact equation for the jump depth ratio. In fact, it can be shown that the mo- 
mentum flux and the turbulence flux corrections tend to cancel each other. 

In this respect, the effect of bottom shear is similar to the effect of baffle 
piers in reducing the downstream depth of a hydraulic jump. (2) In the 
presence of baffle piers the distributed boundary shear parameter S may be 


replaced by a term descriptive of the more concentrated force exerted by the 
piers on the fluid. 
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PHILIP G. HUBBARD, ! A. M. ASCE.—As stated by the authors, the analogy 
between the air-flow system and an actual hydraulic jump may seem to be 
rather artificial at first glance, so that some engineers may want further 
verification or at least an evaluation of the applicability of the results. Be- 
cause of the thoroughness with which the study was conducted and the great 
importance of the conclusions drawn, it would seem desirable to take ad- 
vantage of any technique which might supply additional evidence. The hot- 
film anemometer mentioned in that paper is one such technique, and some 
exploratory measurements with that instrument are presented herein. 
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Fortunately, the flume used by Dr. Siao for the measurements in water is 
a permanent item of equipment at the Institute, so that it was possible to 
reproduce the conditions of those earlier experiments. In order to obtain 
representative results in a brief period of time, measurements were made at 
each Froude number for which results were reported, but only for one repre- 
sentative zone in the flow. The zone chosen was near the downstream end of 
the roller, where the longitudinal component of turbulence was measured 
along a vertical line. Although a complete comparison cannot be made on the 
basis of this one component, it is probably the most representative single 
parameter. As mentioned by the authors, the total intensity of turbulence is 
very near equal to twice the intensity of the longitudinal component: 


When the measuring instrument was placed in the zone chosen, two perti- 
nent characteristics were immediately evident: The passage of innumerable 
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Fig. B. Record of Velocities Inside Roller, Showing Tranquil Intervals 


air bubbles past the probe had scarcely any discernible effect on the velocity 
indications, so that some earlier misgivings were found to be without basis; 
and the fluctuations were so erratic that only a rough indication could be ob- 
tained of either the mean velocity or the root-mean-square deviation of the 
instantaneous velocity from that mean. Although this problem had been en- 
countered to some extent in the air jump, it was greatly compounded in water. 
It had been expected that the frequency spectrum of the turbulence would be 
shifted downward because of the lower velocities used in water, and this shift 
was expected to cause difficulty because the time constants of the turbulence- 
analyzing instruments were marginally small. Even when the probable error 
due to inadequate frequency range in the root-mean-square analyzer was 
taken into account, however, the measured intensities of turbulence shown in 
Fig. A were discouragingly below those shown in Fig. 7 of the subject paper. 

In order to examine more carefully the variation of the fluid velocity, the 
output voltage of the hot-film control circuit was connected directly to the de- 
flection terminals of a cathode-ray oscillograph so that no error was intro- 
duced by coupling networks. Some typical records obtained in this way have 
been traced from photographs for reproduction in Figs. B and C, and it is be- 
lieved that they contain the key to the discrepancy already noted. It was ob- 
served that bursts of high-intensity turbulence were interspersed with periods 
of relatively tranquil flow. Strangely enough, this intermittency existed not 
only near the upper edge of the boundary layer, but was quite pronounced 


Fig. C. Record of Velocity Immediately Downstream from Roller 
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inside of the roller. If the data of Fig. A are modified by eliminating the 
periods of tranquil flow from the record, then the intensity of the remaining 
fluctuations is slightly higher than that indicated in Fig. 7. 

Any evaluation based upon the results reported herein can only be tentative 
pending a more complete investigation, but two conclusions seem warranted: 
(1) The diffusion characteristics of the hydraulic jump are probably repre- 
sented quite well by the results of the subject paper, but those flow charac- 
teristics which depend upon the smaller-scale components are probably 
slightly distorted by the higher indications of the intensity of fluctuation (2) 
Techniques and possibly definitions will have to be developed especially for 
flow situations where intermittency may be a factor. 
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THE TOTAL SEDIMENT LOAD OF STREAMS?# 


Discussions by R. J. Garde and M. L. Albertson, D. C. Bondurant, 
John L. Bogardi 


R. J. GARDE,! J. M. ASCE and M. L. ALBERTSON,2 M. ASCE.—This 
paper is a welcome addition to scarce literature on sediment transport in 
alluvial channels. It will be of use to those who are working in the field, irre- 
spective of the approach in which they are interested—empirical or theoreti- 
cal. 

Broadly speaking, the problem of evaluation of total sediment transport 


rate in alluvial channels can be handled in one or a combination of the follow- 
ing three ways: 


1. Based on the present-day knowledge of fluid mechanics and sediment 
transport, a theoretical relationship for the total load can be developed. 

2. By adding suspended load to the contact (or bed) load, the total sediment 
transport rate can be determined. Einstein’s procedure, and the Modi- 
fied Einstein’s procedure would fall under this category. 

3. Based primarily on dimensional analysis and intuition, some empirical 
relationships can be developed to express the total sediment transport 
rate as a function of parameters relating hydraulic characteristics of 
the flow and characteristics of bed material. 


Because of time and energy involved, and also because of the inadequate 
knowledge of sediment transport phenomena, very little work has been done 
in developing the first method. 

Einstein’s procedure of finding the total sediment load by adding suspended 
load to the bed load is quite often used. It is found that the results obtained 
by Einstein’s procedure do not agree with certain data obtained in the field. (8) 
Therefore, Hembree and Colby(8) have developed a “Modified Einstein’s Pro- 
cedure” to evaluate the total sediment load. 

In evaluating the total suspended load by integration of the sediment distri- 
bution curves, some important assumptions are made. Following are given 
the basic equations and the assumptions made in their development: 


Equation Assumptions 


(uv +é€m)d% (ev) That the density P of the fluid-sediment 
dy suspension is independent of y and that 
Cem AV the kinematic viscosity V is relatively 
oy insignificant compared with the eddy 
viscosity €m. 


a. Proc. Paper 1530, February, 1958, by E. M. Laursen. 
1, Graduate student, Colorado State Univ., Fort Collins, Colo. 


2. Director of Research Foundation and Prof. of Civ. Eng., Colorado State 
Univ., Fort Collins, Colo. 
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Equation Assumptions 
2. & . 2.02 That the logarithmic distribution of 
Ky 


velocity is applicable and that the 
Karman universal constant k has a 
constant value, 


That the sediment transfer coefficient 


€, is equal to the momentum transfer 
4, Em= €s = coefficient 


At present, work is being done at Colorado State University on the evalu- 
ation of total sediment load by adding suspended load to contact load. Thus 
far, this method has been found to have serious limitations because in a dune 
regime one does not know how to predict accurately either the concentration 
of suspended load at a given elevation or the exponent z in the actual sediment 


distribution curve. Furthermore, if one assumes a general law of velocity 
distribution as 


Y= Log 


Ve Ks (1) 


in which v_ is the velocity at a distance y above the bed 
k is the Karman universal constant 
Kg is some length parameter 
V* is the shear velocity 


it is found that both k and K, depend on the flow regime—such as plane bed or 
dune bed—and other flow parameters. 

Therefore, for immediate use, at least some empirical relationships are 
necessary for predicting the total load. With this view in mind, the writers 
carried out certain preliminary investigations in 1956. The results of these 
investigations are summarized in the following paragraphs. 

While studying the sediment transport through pipes, (1) the writers found 
by dimensional ‘analysis the following relationship to be applicable: 


ReVf (2) 
in which D is the diameter of the pipe 
d is the diameter of sediment particles 


Cy is the concentration of total sediment load expressed in per cent 
by weight 


is the Darcy-Weisbach resistance coefficient and 


Re is the Reynolds number = yD 


Vv 


For smooth pipes (as compared to standard:corrugated, or helically corru- 


gated pipes), for a given value of >. it was possible to obtain plots as shown in 
Fig. 1. 
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Although open-channel flow differs from pipe flow because of the existence 
of a free surface, it is reasonable to assume that a similar relationship 
should exist for the transport of sediment in alluvial channels. Therefore, 


Y*D (which is proportional to Ref ) was plotted against Cy. A typical vari- 


ation of “eB with Cy and » is shown in Fig. 2 using Gilbert’s data. 

For each set of data, a definite relationship existed between YD p and 
Ct. Plots were made for flume data from various sources such as Kalinske 
and Hsia,(3) Brooks,(5) Barton and Lin,(6) and Gilbert.(4) The range of size 
of bed material included in the study was from 0.011 mm to 1.71 mm. 

Using standard techniques, it was possible to express the total sediment 
transport rate C; by the following empirical relationship: 


(“eP. +, ) (3) 

in which ¥ is a function only of size of bed material. This functional relation- 

ship is shown in Fig. 3. 

Fig. 4 shows Eq. (3) plotted for the various data already mentioned. 
(Laursen’s data for 0.04 mm and 0.10 mm are also included.) There is an 
average scatter of 40 to 50 per cent. Possible causes for this scatter are 
that, in the development of Eq. (3) neither the influence of flow regime (such 
as dune bed or plane bed) nor the influence of the Froude number on the total 
sediment transport are taken into consideration. 

The investigations under progress at Colorado State University by the 
writers have shown that k, Kg, and z depend on the condition of the bed and, 
therefore, the rate of suspended load transport should depend on the flow re- 
gime. Similarly, Tsubaki and others(7) have shown that under conditions of 
equal tractive force the contact load is smaller if ripples exist on the bed than 
if the bed is plane. Therefore, flow regime at least should be taken into con- 


sideration in the development of formulae for transport of sediment in alluvial 
channels. 
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D. ©. BONDURANT, | A. M. ASCE.—There are several formulae available 
for the computation of sediment transport. Unfortunately, data available for 
the development of these formulae have seldom extended to the range of large 
streams, and it has been found that their use in respect to such streams re- 
quires a liberal application of experienced judgment. The writer’s primary 
objective in reviewing this presentation was, thus, to determine its values and 
limitations in respect to streams such as the Missouri River. It is concluded 
that the form of the author’s equation can be very useful but that some re- 
vision of the curve of is indicated. 

Data for the purpose of the review were available from measurements 
made in the Missouri River in connection with the sediment investigations of 
the Missouri River Division, Corps of Engineers. Of the more than 400 sets 
of measurements of vertical distribution of velocity and suspended sediment, 
76 were selected to provide the maximum range of data. Within these data, 
stream depths ranged from 7 feet to 22 feet, velocities from 3 fps to 9 fps, 
and discharges from 25,000 cfs to 90,000 cfs. Stream widths approximated 
1,000 feet and slopes were from 0.8 to 0.9 foot per mile. From 6 to 8 verti- 
cals were occupied at each range, and the ranges occupied were in 4 separate 
reaches located at (1) River Mile 640, above Omaha, Nebraska, (2) Mile 614.1, 
2-1/2 miles above the mouth of the Platte River, (3) Mile 608.2, 3-1/2 miles 
below the mouth of the Platte River and (4) Mile 365.6, below Kansas City, 

Mo. and 14 miles below the mouth of the Kansas River. Reaches were as 
nearly straight as could be found, and were unobstructed. Measurements were 
made from a boat held in place on a tag-line. 

Figs. 15, 16, and 17, respectively show typical plots of velocity, suspended- 
sediment, and bed-material. Velocities were measured with a Price-type 
current meter; suspended sediment samples were obtained with a U. S. P-46 
(point-integrating) sampler; and bed-materials samples were obtained with 
either a B.M.-48 or B.M-54 sampler or both. The B.M. instruments sample 
only the top one-inch (approximately) of the bed, and are so constructed that 
the sample is effectively protected against washing or contamination. 

Fig. 18 shows a comparative plot of the average suspended sediment 
concentrations at the three upstream reaches as determined from the samples 
and as computed by the author’s method. The measured concentrations are 
four times as large as the computed; however, the relationship is reasonably 
consistent. Furthermore, the measured and computed size distributions were 
in excellent agreement. Additionally, four points provided by Mr. J. C. Pyle 
of the U. S. Army Engineer District, Little Rock, Arkansas, showed almost 
exactly this same relationship of four. These latter points represented dis- 
charges of 25,000 cfs, 73,000 cfs, 135,000 cfs, and 190,000 cfs in the Arkansas 


1. Head, Sediment Section, U. S. Army Engr. Div., Missouri River. 
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River. The writer has no information as to the magnitude of other factors a 
involved. 
On the basis of the consistency indicated in Fig. 18, it appeared that ug 
further computations might be unnecessary. However, the bed-materials in 
the Kansas City reach, as may be noted in Fig. 17, are substantially larger ¥ : 
than those of the upstream reaches, thus they provided an opportunity to ex- . 3 
tend the range of the data used. With these data the consistency of the re- . 
lationship deteriorated rapidly, the ratio varying from 10 to 100 and indicating i 
the need for a more basic review. ” 
° Inspection of the computations revealed that the values ot ¥%/p/ associ- ef 
ated with the larger bed-materials of the Kansas City reach caused the corre- i 
sponding function values for the major size fractions to be much smaller than 
appeared reasonable. Accordingly, re-computations were made to derive 
values of the function corresponding to the measured data. a 

Figs. 19 and 20 show plots of f(U,/w) computed for the upstream reaches ” 
and Kansas City reach respectively. Since all the writer’s plotted data use a 
the symbol “U,” for the more formal V%/p , that symbol will be used for this 
discussion. It may be seen that the points plot with a reasonably limited 
scatter except for those points representing the two smaller size fractions 
and a few points representing the largest sizes in the respective samples. 
These deviations point to a systematic error in the data in respect to these ap 
fractions, and investigations made subsequent to the analyses of the samples Ms 7 
indicated that such error did exist. 

The samples were analyzed by sieves, for reasons outside the scope of the 
present discussion. Fall velocities at various temperatures were later de- 
termined by an extensive investigation involving both individual grain drops ' 
and analyses of individual sieve fractions by the visual accumulation tube . 
method. It was noted that the smaller size fractions periodically contained 7 . 
quantities of heavy magnetic minerals and, if the fraction represented only a . 
small percentage of the total sample, the average specific gravity of the 
fraction might be as much as 2.85 and in a few cases greater than 3. Had ; 
these points been corrected for the higher fall velocity which would have re- = 
sulted from the presence of these minerals, they would have been moved to 
the left and would have been essentially in line with the other points, bi: 

The less consistent variations for fractions of the larger sizes were found 7 
to be due to low density agglomerated particles. They were present only in 
sufficient quantity to affect a fraction which was less than about one per cent 
of the total sample. 7 

It might be well to note here also that the percentages pertaining to the 
size fractions of the bed-material samples were read from the plotted curves tte 
rather than from the original laboratory records. Inspection of the compu- 
x tations showed that the possible error thus introduced was negligible insofar 

as the concentration computations based on Missouri River data were con- 
cerned; i.e., those fractions for which the error of reading might be appreci- Bs 
able were of themselves such a small portion of the whole that a reading * 
" error of even several hundred per cent would make little difference in the a 
over-all result. Time did not permit a return to the original analyses for the a 
function computations. 

Fig. 21 shows individually the lines drawn as the approximate average of 
the plotted points for the upstream reaches and the Kansas City reach. The 
spread between these lines is not excessive, however, they are roughly paral- _ 
lel, indicating a systematic difference. A third line on this figure indicates an  e 
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average for all points representing fractions which were more than 5 per cent 
of the total respective samples, but this line is not presented as a recom- 
mended mean. The difference between the lines for the upper and lower 
reaches is believed to result primarily from the limits used to determine the 
average concentrations from the measured data. 

If one assumes that the logarithmic distribution of concentrations extends 
to y = 0, the concentration at that point becomes infinite; thus, one must as- 
sume a limit at which y equals some real value. For purposes of prior 
computations, this limit had been adopted as y = 2d with an arbitrary maxi- 
mum of 480 g/1 which was not to be exceeded at that level. The average 
concentrations thus determined were used for the present computations. In- 
sofar as the three upper reaches were concerned, these averages were 
reasonably consistent with averages shown by depth-integrated samples 
(which could not extend closer than 0.5 foot to the bed). In the Kansas City 
reach, where the material was larger and thus necessarily concentrated more 
heavily near the bed, the difference was substantial. 

Experiments of Einstein and Chien26 indicate that the logarithmic distri- 
bution of suspended sediment begins to break at the level where the total 
concentration approximates 100 g/1. Inspection of the computations for aver- 
age concentrations at Kansas City indicates that assumption of any limiting 
value below 480 g/1 would reduce most of the computed average concen- 
trations, and assumption of a value below 350 g/1 would reduce all of them. 
This would move the computed values of {(U,/w) for the Kansas City data to 
the right, more nearly approximating the position for the upstream data. 

Time has not permitted recomputation of the data to verify, quantitatively, 
the variations above discussed. Even with these variations, however, the 
Missouri River data provide a consistent plot within their range, without ex- 
cessive scatter, and verify the form of the author’s equation. A plot of the 
author’s data (roughly scaled from Fig. 7), included in Fig. 21, indicates that 
the average lines from the Missouri River computations fit the author’s data 
for lower values of U,/w as well as or better than the curve necessarily con- 
structed by the author to fit his data for the higher values of U,/w. For the 
higher values of U,/w the Missouri River results deviate sharply from those 
of the author. 

The Missouri River data are sharply restricted in their range in respect 
to both sediment size and slope of the energy gradient. The writer does not, 
thus, suggest that a function derived from these data be given arbitrary pre- 
cedence over a function derived from other data. It is to be hoped that other 
readers will further augment the available information to permit the author 
to provide an integrated result. 

In considering the fit of the Missouri River data, it is well to remember 
that all pertinent factors in these data represent measured values in which 
the correct values of other functions, such as roughness or turbulence, are 
implicit. It has been found that equally valid computations can be made with 
the theoretical type of equations (for example) if the measured values of k (or 
8k) are utilized. The author’s function has the advantage that velocities can 
be more simply measured than k values; thus, in cases where the velocity 
can be adequately established, this function should be the more valuable. In 


26. H. A. Einstein and Ning Chien, Effects of Heavy Sediment Concentration 
Near the Bed on the Velocity and Sediment Distribution. University of 
California Inst. of Engr. Research Series No. 33, Aug. 1955. 
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cases where velocity must be computed from assumed roughness factors, the 
necessary dependence upon judgment and the possible variations in results 
would be about the same as for other available formulae. 


JOHN L. BOGARDI.!—All credit is due the author for his successful in- 
vestigations into the subject of sediment transportation. In keeping with the 
latest concepts, bed-load and suspended load are not separated in the 
treatment, and the possibility for relationships of a more general nature is 
opened thereby. Numerous aspects of this complex problem are elucidated 

by clear and logical reasoning. The writer fully agrees with every statement 
made by the author and only wishes to contribute a few additional remarks, 

In the writer’s view, the relationship shown in Fig. 7 applies to the 
sediment-transporting capacity of the stream rather than to the actual rate of 
sediment transport. Quite obviously, the actual rate of transport will depend 
on a number of factors besides the transporting capacity, e.g. also on the 
availability of sediment in the stream. It follows that the actual rate of 
transport in natural streams is governed also by the hydrological conditions 
of both the watercourse and the catchment area drained thereby. The writer 
interprets the excellent agreement between predicted and actually observed 
rates of transport on the Niobrara River as an indication of the fact that 
hydrological conditions of the watercourse, as well as sediment supply, are 
such as to enable the stream to carry sediment in an amount equal to its 
transporting capacity. A rate of transport in excess of the transporting ca- 
pacity and observed on Mountain Creek is, as mentioned also by the author, 
due probably to artificial flash floods. Consequently, the author is fully justi- 
fied in arriving at the conclusion that the prediction of the so-called wash load 
is not feasible. 

Relationships (8) and (9) published by the author gave the writer the oppor- 
tunity to extend his investigations, conducted on the basis of results published 
by Liu(1) and by Albertson, Simons, and Richardson,(2) to total sediment 
transport. In his paper Liu established a definite relationship between U,/w 
and U,d/v for the state of incipient ripple formation. Albertson and his 
collaborators have further shown the existence of a similar relationship be- 
tween the two above dimensionless invariants also for the beginning of 
movement, the development of dunes, transition, and antidunes, respectively. 


The writer investigated the significance of the channel stability factor 
ns U 
* 


used in Hungary since 1942. This factor, being formally a Froude number, 
was named, following the example set by Albertson and in analogy to the shear- 
velocity Reynolds number, the shear-velocity Froude number. Using the 
channel stability factor b = 1/Fr, the writer succeeded in defining numerically 
various regimes of movement by a single parameter 8. Investigations of the 
writer, conducted using data and results published by Albertson and his col- 
laborators, showed a relationship between the parameter # and the shear- 
velocity Froude number that can be expressed as 

pa abet = 

hs Fr, 

1. Asst. Prof., Technical University, Budapest. Head, Experimental Section 
of the Scientific Research Inst. of Hydroeconomy. 
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Since for all regimes N was 0.882 = constant, the parameter 
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cd 
(3 = — = 


x Vx 
defines uniquely the bed formation. The relationship between the particle 
size d and 1/Fr, is shown in Fig. 15. It is to be seen from the figure that the 
parameter # (of the dimension cm-0.882) in effect defines uniquely the bed 
formation. 

A number of important conclusions can be arrived at from Fig. 15. 
Ripples, for instance, cannot develop unless the particle size is smaller than 
approximately 2 mm, and even for these particle sizes only if 66< 8 < 322. 

Fig. 15 relates to a water temperature of 20° Centigrade and to a sediment 


having a specific gravity v, = 2.65 gr/cc. For other temperatures, e.g. for 
a temperature of t°® Centigrade, 


2 2 
t 20 v 2 Fr, 


wherein V9 and vt denote the kinematic viscosity of water at 20° and t° 
Centigrade, respectively. 
For sediments having any arbitrary specific gravity V1 
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The interrelation between U, and d can also be determined from that es- 
tablished for B and Fr,. Thus 


a 


Since, however, 1 = 0.0582 = constant for every bed configuration, 


| 


of the dimension cm9-9418 sec™l, results in the relationship 


* 
This relationship, derived also on the basis of data published by Albertson, is 
shown in Fig. 16. As can be seen from the figure, the configuration of the 
bed is uniquely defined for any arbitrary particle size by the parameter 


introducing the parameter 


It should be noted that all sorts of bed configurations are uniquely defined 
similarly to parameters f and € by the parameter 


aA = 


5.505 
ne 
v 


The parameter € is also related to a water temperature of 20° Centigrade 
and to a sediment having a specific gravity vy = 2.65 gr/cc. For water 
temperatures and sediments having a specific gravity other than mentioned 
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above, the parameter can be determined from the following two equations: 


1856-77 


Details pertinent to the above investigation are to be found in the writer’s 
discussion of Liu’s paper. 


In the writer’s opinion, if Fig. 7 presented by the author actually expresses 


the relationship between the sediment-load parameter corresponding to the 


transporting capacity of the stream and U//w, the variation in the sediment- 
load parameter 
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can also be entered into Figs. 15 and 16. 

Values of U,/w can be established for all five bed configurations in the 
case of several particle sizes from Fig. “A” published by Albertson and his 
collaborators;(2) moreover, values of U,d/v can also be determined. The 
latter lends itself readily to the computation of values of U, and subsequently 
of those of the parameter gd/U,2. For each value thus determined, those of 


pertaining to various U,/w values can be correlated using Fig. 7 published 
by the author. 


The relationship between 


can thus be represented in Fig. 17 for particle sizes in Albertson’s Fig. “A”. 


Thereafter, relying upon Fig. 17, the isometric chart of the sediment load 
parameter 


can be traced in Fig. 15. Finally using Fig. 15, the isometric chart of the 
sediment load parameter can be entered into Fig. 16 as well. 

It should be remembered that Figs. 15, 16 and 17 were prepared on the 
basis of values obtained from Fig. “A” and Fig. 7 published by Albertson and 
the author, respectively, as no numerical values were available to the writer. 
Some discrepancy may easily be introduced thereby. 

It should be noted further that the writer completed Albertson’s Fig. “A” 
by extrapolating Shield’s curve towards smaller particle sizes, and, on the 
basis of Fig. 15, slightly modified at the diameter d = 0.0011 cm the curve 
characterizing the five bed configurations by transferring observation points 
of Kalinske and Hsia to higher U, /w values. 

The writer considers Figs. 15 and 16 very useful for design purposes, 
since for any given watercourse, having established the hydraulic radius 
(water depth), the slope, and the particle size of the bed-material, parameters 
B and € may be computed directly, and the value of one of the two is alone 
sufficient to define uniquely the bed configuration, and also the value of the 
sediment-load parameter, which in turn permits the computation of the total 
sediment-transporting capacity of the stream. 

It should be noted further that no accurate limiting values have been de- 
termined towards large and small particle sizes in deriving the relationships 
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a and in constructing the figures. These are felt to require further study and 1 
a experiments. For example, the line representing the formation of antidunes ‘" 
in Fig. 15 is likely to join that representing the limiting condition of a smooth 
bed through a curved transition instead of a definite point of intersection. 
Similarly, the determination of relationships for values smaller than 7 and 
. higher than about 1000 to 3000 of the sediment-transport parameter is likely t 
to require further investigation. " 
The number of particular and general conclusions that can be arrived at ‘ 
by means of the relationships described above and by Figs. 15, 16, and 17 con- ‘a 
< structed therefrom is great. It is, however, the intention of the writer to 
devote a special paper to these problems. 
REFERENCES 
- 1. Liu, H. K., “Mechanics of Sediment Ripple Formation”, Proceedings of if 
. : ASCE, Hydraulics Division, Vol. 83, No. HY 2, April, 1957. af 
‘ 2. Albertson, M. L., Simons, D. B., Richardson, E. V., “Discussion of Sedi- a. 
ment Ripple Formation” by H. K. Liu, Proceedings of ASCE, Hydraulics a 
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EISENHOWER AND GRASS RIVER LOCK MODELS#@ 


Discussion by Marvin J. Webster 


MARVIN J. WEBSTER,! M. ASCE.--The authors state that the degree of 
surging and consequently the magnitude of force imposed upon a craft in the 
lock chamber reaches a minimum value when the manifold system of ports or 
laterals occupies the approximate middle third of the lock chamber walls. 
Tests performed by the Bonneville Hydraulic Laboratory, Bonneville, Oregon, 
in a 1:25-scale model of The Dalles navigation lock, Columbia River, and in 
the prototype furnish additional data to substantiate this statement. 

The Dalles lock has a culvert system in which seven lateral culverts from 
each wall culvert are located in the middle third of the lock chamber (see 
Fig. 9). Model and prototype tests showed very low hawser forces on a com- 
mercial tow of 2780 tons.(1 

The tests in the prototype were made in December 1957.2 The tow consist- 
ed of two barges and a tug having a total displacement of 2780 tons. Instru- 
mentation for the tests was developed and constructed at the Bonneville Hy- 
draulic Laboratory. The bridles for holding the tow in position were con- 
structed of 5/8-in. cable and connected to the floating mooring bits, as shown 
in Fig. 10. The tow was centered over the lateral culverts. The longitudinal 
forces were measured with two 20,000-lb-capacity dynamometers, one at 
each end of the tow. No attempt was made to measure transverse forces in 
the prototype because of the difficulty of constructing a suitable bridle to 
measure transverse forces and because previous tests in the model on larger 
size tows showed that transverse forces were less than the longitudinal forces. 
Tests were made in the 1:25-scale lock model prior to the prototype tests 
and, although the displacement was simulated, the shape of the model tow 
differed from that of the prototype. 

Results of the tests during filling for a lift of 84.6 ft in the prototype and 
86.0 ft simulated in the model are shown in Fig. 11. The maximum lift of 
The Dalles lock is 90.5 ft. The test results show low hawser forces during 
filling operations on the tow used in the tests and good correlation between 
the model and prototype for filling. An emptying test in the prototype with the 
same tow showed a maximum longitudinal force of 1.1 ton. An emptying test 
was not made in the model. Hawser forces also were measured on the same 
tow in the McNary navigation lock, which has a culvert system similar to The 
Salles, for a lift of 87.9 ft (maximum lift is 92.0 ft). The maximum 


a. Proc. Paper 1582, April, 1958, by Martin E. Nelson and Harvey J. Johnson. 

1. Chf., Hydr. Design and Lab. Section, U. S. Army Engr. District, Portland, 
Ore. 

2. Ona tow furnished through the courtesy of the Tidewater-Shaver Barge 
Lines, Portland, Oregon. 
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longitudinal force measured during filling was 1.7 tons. The filling valves 
open in 7.6 min and the lock filling time is 15.0 min at McNary compared with 
4 min and 12.5 min at The Dalles, which probably accounts for the lower 
hawser forces at McNary. 

The writer wishes to acknowledge the use of material for this discussion 
from the U. S. Army Engineer District, Portland, Oregon. Instrumentation 
for the tests was developed and constructed under the supervision of C. 
Robert Grim, Head of the Electronics Unit under the general supervision of 
Harry P. Theus, Chief, Bonneville Hydraulic Laboratory and the writer. 
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1. Memorandum Report 1-1, Prototype Hawser Forces, The Dalles and 
McNary Locks, Bonneville Hydraulic Laboratory. 
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WATER DISTRIBUTION DESIGN AND THE 
McILROY NETWORK ANALYZER@ 


Discussions by Quintin B. Graves and Don Branscome, E. Shaw Cole, 
Claud C. Lomax 


QUINTIN B. GRAVES,! M. ASCE and DON BRANSCOME,?2 J. M. ASCE.— 
The authors have presented an interesting approach to the design of water 
distribution systems with particular emphasis on the advantages of the 
McIlroy Network Analyzer in providing a visualization of the hydraulics of a 
water distribution system. 

This paper and a paper by the writers on “Digital Computers for Pipeline 
Network Analysis” were presented at the October, 1957 ASCE Convention, 
but before different divisions and were published simultaneously in the Pro- 
ceedings. 

It is obvious, as the authors point out, that by visual inspection of the 
system as set up on a McIlroy Analyzer design changes will be quickly indicat- 
ed and may be easily adjusted at “Inflow rates well below the design rates”. 
This is an advantage for the McIlroy Analyzer over the electronic digital 
computer primarily in those cases where changes in the distribution of in- 
flows or outflows are indicated. In the case of the electronic digital comput- 
ers, computations on the basis of the largest rates of demand usually will 
prove to be advantageous. This may not be required in the case of the McIlroy 
Analyzer. 

No doubt, improvements in the techniques of utilizing the McIlroy Analyzer 
have been made since its first introduction, and more probably will be made 
with continued use. Also, improvements in programming for digital computer 
analysis of distribution systems will be made. It has been evident from ‘long 
hand” calculations with the Hardy Cross Method that convergence of the 
process does not always occur. Usually, by changes in the initial flow as- 
sumptions for the network it is possible to obtain rapid convergence. The 
writers are in agreement with the authors in regard to the use of “electric 
service power-system calculating boards with linear resistors”. 

One of the advantages of electronic digital computers over other types of 
analyzers is their wide distribution throughout the United States. One of the 
costs to occasional users of McIlroy Analyzers is that of travel and associated 
expense for those who do not have one available for use locally. For most 
users of digital computers, such costs will be substantially reduced due in 
part to the greater number of such computers available. As greater use is 


a. Proc. Paper 1588, April, 1958, by M. B. McPherson and J. V. Radziul. 
1. Prof. of Civ. Engr., Oklahoma State Univ., Stillwater, Okla. 
2. Sales Engr., Aluminum Co. of America, Pittsburgh, Pa. 
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made of digital computers for water distribution analysis, many improvements 
in programming will be made thus enhancing their value. 

The authors mention that 50 to 60 punched cards will be required for basic 
instructions, regardless of the size of the network being studied. It is possi- 
ble to reduce the number of basic instructions required by using an inter- 
pretive programming system. For most networks only 30 to 40 punched cards 
will be required for the basic instructions. 

The authors, in discussing Fig. 6, state “A revised storage rate requires 
a new storage curve”. The writers are of the opinion that this would be more 
intelligible if the statement read, “A revised uniform pumping rate requires 
a new storage curve.” 

Values of the Hazen-Williams coefficient “C” are not accurately known for 
many small water distribution systems. In such cases the variation in 
voltages or currents, which has been mentioned to the writers by some users 
of the McIlroy Analyzer are of little or no concern. In those few cases where 
basis information is more accurate, a variation of 5 per cent or greater may 
be significant. 

The authors’ statement that “(However, a 2.00 power rather than 1.85 
power is involved)” in reference to the electronic digital computer programs 
is no longer valid as programs are available utilizing the 1.85 power. 

The authors hold that “an error in programming or data instructions is not 
evident until the run is made and the analysis is completely summated and 
checked, ...” In general after a program is once checked out all errors 
should be eliminated and would not introduce errors when applied to other 
water distribution systems. Errors in data may occur for any given case, 
however automatic machine checking of data such as the summation of inflows 
and outflows at a junction may be included in a program prior to the analysis 
of the system. 


The authors are to be highly commended for this fine contribution to the 
literature on water distribution system design. 


E. SHAW COLE,! M. ASCE.—The authors have made an important contri- 
bution in water distribution design and have shown how complicated problems 
can be solved readily by use of the McIlroy Analyzer. The writer has been 
privileged to work with them on some of their problems and has used the 
analyzer on a number of design problems for other systems. He believes the 
authors have presented thoroughly the advantages of the analyzer, as well as 
the limitations in its use, but would like to mention several points that may 
not be apparent to the casual reader. 

The importance of accurate field data is mentioned. This should be obvious 
and was recognized by Dr. McIlroy at the outset of his work in developing the 
analyzer. Unless the C-Value is reasonably accurate, the answers given by 
the analyzer may be seriously in error, but accurate field data implies more 
than determination of C-Values. 

Frequently, water distribution systems do not function as intended, due to 
closed valves or obstructions that reduce carrying capacity far below the 
most conservative estimates. Some water works designers may not appreci- 
ate the difficulties of maintaining distribution systems at peak efficiency, but 


1. Pres., The Pitometer Associates, New York, N. Y. 
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the discovery of major difficulties is becoming more frequent every year. } 

The increasing complexity of our systems and the lack of responsibility, or 

complete records of the field operations, are major causes for this condition, Re 
The problem of determining the quantity of equalizing storage necessary a 

to permit uniform pumping is most interesting and illustrates the advantage - 

of the analyzer. In many water works systems, uniform pumping is not Be 

practical and pump units are added or subtracted periodically as the load 

changes. The pump changes may only occur a few times during a twenty-four 

hour period but this can be sufficient to make effective use of the existing 

storage. For medium size systems, the cost of pumps is relatively much less 

than storage, particularly if elevated, and it would be interesting to have a 

study made of the economics of variable pumping under various conditions 

vs. uniform pumping. 


CLAUD C. LOMAX,! A. M. ASCE.—This paper is a welcome addition to the 
literature on the analysis of fluid distribution networks. Use of the McIlroy 
analyzer to study the effect of demand ratio and pumping station outages, and - 
to balance pumpage and storage demonstrates its versatility. As engineers 
become more familiar with the potential of this instrument, increased use for 
other than the routine head loss and flow rate determinations can be anticipat- 
ed. The R. L. Albrook Hydraulic Laboratory of the Division of Industrial " 
Research, Washington State College, has had such an analyzer since 1951. Se 
This analyzer has been used for high and low pressure gas and for water o 
distribution networks. The usual analysis has been determination of head 
losses and/or flow rates for present and future demands. Unusual problems 
include analyses for Anchorage and Fairbanks, Alaska to insure all pipes had 
the minimum flow required to prevent freezing. The distribution of cooling 
water for Priest Rapids turbines was also analyzed. The City of Spokane in- 
cluded the pump characteristics as influenced by the system pressure. Bend, 
Oregon used this analyzer to assist in selecting a location and water surface 
elevation for a reservoir being added to their system. The City of Los 
Angeles used the analyzer to study drawdown characteristics of their reser- 


voirs, Pressure drops and distribution of bleed air for B-52 jets were a 
checked, 


Other possible uses include economic studies related to pump selection, a 
reservoir size and location, and main size. Stagnation potential of reservoirs 4 
can be evaluated. Leaks and blocking can be pinpointed when a few field + 
pressures have been taken. Hazen-Williams C values can be determined or a 
confirmed from field data for a badly tuberculated system. In time of dis- d 


aster this analyzer could be invaluable in re-routing gas and/or water for the 
most efficient use of the remaining system. A booster pump which will raise b 
pressure at a selected point in the network has been added to the R. L. Albrook J 
Laboratory’s installation. 

The authors’ comparison of the digital computer and McIlroy analyzer con- 
firms several of the conclusions of this Laboratory. The great value of the 
analyzer over a digital computer in presenting an immediate and visual analy- 
sis of a system should be emphasized. While such a picture is very useful to 
the engineer, it is also valuable in helping administrative personnel, fire 


1. Associate Hydr. Engr., R. L. Albrook Hydr. Lab., State College of 
Washington, Pullman, Wash. 
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chiefs and water superintendents understand the distribution systems. The 
analyzer is probably a better design tool than the digital computer since the 
effect of changes are easily and quickly evaluated. 

There are disadvantages to the McIlroy analyzer which are not easily 
solved. The number of pipe sections which can be included are limited by the 
available Fluistor sockets. The identical pipe sections which can be placed 
on the analyzer are limited by the Fluistor stock, and it is not practical to 
maintain a large stock of identical Fluistors for several pipe sizes. Solution 
of the Darcy equation takes more time, and a trial and error fitting of the 
Fluistors is necessary. 

A program for the IBM 650 digital computer to supplement the McIlroy 
analyzer has been developed by the R. L. Albrook Hydraulic Laboratory and 
the Computing Center at the State College of Washington. It is believed that 
the computer will be better for system analyses when there are few changes 
anticipated and when (1) the Darcy equation is used, (2) the piping system con- 
sists of many identical sections, (3) the number of pipes greatly exceeds the 
analyzer capacity and/or (4) great accuracy is required. The combination of 
the analyzer and the computer will provide complementary rather than com- 
petitive services. 

The possibilities of the digital computer and the analyzer have not been 
fully developed. As increased use of each continues it is hoped future com- 
parisons will be made. There seems to be an established need for each. 


Factors influencing the selection of one or the other are not clearly es- 
tablished at this time. 
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NORTHEASTERN FLOODS OF 1955: METEOROLOGY 
OF THE FLOODS4 
Discussion by H. Alden Foster 


H. ALDEN FOSTER, ! M. ASCE.—Hydraulic engineers and hydrologists are 
particularly interested in the effects of hurricanes or other large cyclonic 
disturbances on the production of floods; hence it may be in order to point out 
certain lessons that may be learned from the authors’ analysis of the storms 
of 1955. 

One interesting fact is the extreme magnitude and expanse of the forces 
involved in such cyclonic storms. Thus, the total volume of water carried 
over the continent in one 24-hour day of hurricane Connie is given as 562,000 
square-mile inches, or about nine cubic miles,—a lot of water even if spread 
over a large land area. Also, during a single week in the life of such a hurri- 
cane, a large portion of the atmosphere over the tropical Atlantic Ocean could 
be drawn into the storm pattern. Evidently, hurricanes mean not only ex- 
cessive winds but also immense volumes of water that may be released to the 
streams and rivers in a short period of time. 

There is an extreme variation in location of the tracks of these tropical 
storms. They originate over the open ocean far south of the northeastern 
states, but during their travels may end up almost anywhere in this area. 

Their intensity is maintained largely because atmospheric conditions near the 
North Atlantic coast during late summer or early fall may at times be simi- 
lar to those commonly occurring in tropical regions where these storms 
originate. 

The travel paths of hurricanes are quite unpredictable for more than a day 
or two in advance. These paths are controlled by the general tendency of all 
atmospheric masses in motion to curve to the right in the northern hemi- 
sphere, but they are also strongly affected by atmospheric conditions in their 
vicinity, such as occurrence of high or low pressure centers or troughs. 

The storm may pass over the open ocean throughout its journey to the 
northern latitudes, in which case it is most likely to cross the coastline of 
Long Island or southern New England. Or it may pass inland south of Cape 
Hatteras, and wear itself out while passing northward more or less parallel 
to the coastline. Where the storm track crosses the coastline, exceptionally 
high tides may be expected, not necessarily at the location of the path of the 
center or “eye” of the storm but more likely some distance ahead of this lo- 
cation,—that is, further north or east. 

The precipitation rates caused by a hurricane are not necessarily a function 
of the wind velocity. Even after the storm has grown old and is about to 


a. Proc. Paper 1661, June, 1958, by Charles S. Gilman and Kendall R. 
Peterson. 


1. Prin. Assoc., Parsons, Brinckerhoff, Hall & MacDonald, New York, N. Y. 
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expire, it may still carry an immense load of moisture in the circulating 
atmosphere and is quite capable of dumping this water in places where it can 
cause severe floods. 

Theoretically, the major portion of the precipitation should take place 
ahead of and somewhat to the right of the eye of the storm, although topo- 
graphic or orographic influences may tend to affect this tendency; or the air 
masses carrying excessive quantities of moisture may be deflected from their 
natural paths by atmospheric conditions in the general vicinity. As a result, 
the track of the storm does not necessarily indicate the location of the maxi- 
mum total precipitation during the storm. For instance, the total precipitation 
for Connie was generally to the right of the track, while for Diane most of the 
precipitation in southern New York and New England appears to the left of the 
track. This seems to be due to the fact that, on August 16 when Diane was 
moving northward over Virginia, the precipitation was mostly located far to 
the northeast, from Pennsylvania all across southern New England. By the 
next day, Diane had made an abrupt turn to the right and passed south of the 
area where the rain had mostly fallen the day before. If Diane had not changed 
her course so abruptly but had followed a track eastward near the southern 
boundary of Massachusetts, the precipitation in southern Connecticut might 
have been considerably increased by bringing in some of the water which actu- 
ally must have fallen over the ocean south of Long Island in the right-front 
quadrant of the storm as it moved along. 

This discussion brings up the important question as to whether the storms 
of 1955 produced the maximum precipitation that may ever be expected in the 
affected areas. Studies made in 1955 by Dr. A. A. Koch under the writer’s 
supervision show that no storm of record in these Northeastern states has 
yet closely approached the maximum possible value for precipitation. It is 
quite possible that greater precipitation and more severe floods may occur in 
these areas in the future, though the probability of occurrence of such floods 
may be very small. On the other hand, the maximum possible storm tides 
have been approached on several occasions at certain locations on the coast 
of southern New England, although other Atlantic coastal regions have not yet 
experienced the greatest possible storm tides. 

It is interesting to note in a news report of August 23, 1958, that the 
Weather Bureau has started a comprehensive study of storm tides along the 
3,000 miles of the Atlantic coast, from Texas to Maine. It is to be hoped that 
such a study will show whether the probability of occurrence of hurricanes in 
this region may be greater in the future than indicated by past records. 
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NORTHEASTERN FLOODS OF 1955: FLOOD HYDROLOGY@ 


Discussions by Gordon R. Williams and H. Alden Foster 


GORDON R. WILLIAMS,! M. ASCE.—Ample evidence is given in this paper 
that floods resulting from extraordinarily large concentrations of precipi- 
tation are hydrologically and hydraulically a different phenomenon from floods 
having recurrence intervals of 10 years or less. No longer should engineers 
and mathematicians compile frequency relations on the assumption that all 
floods comprise a homogeneous statistical population. The writer has previ- 
ously emphasized this point. (1) For river basins with areas of less than about 
4,000 square miles in rolling or mountainous country, the major floods often 
define frequency curves which have steeper slopes than the curves for ordi- 
nary floods. These relations are particularly pronounced in many 
California(2,3) streams and elsewhere.(4) Rowe‘3) has used the term bilineal 
frequency curve. Curves with this characteristic are also common in New 
England streams.(5) Unfortunately the number of points available to define 
the upper curves are too few so the problem of determining reliable recurrent 
intervals for rare floods is not solved. 

The most promising approach to frequency determinations is the regional 
or composite frequency curve method that has been adopted by the U. S. 
Geological Survey(6,7) and others. (3) Briefly this method uses the mean 
annual flood or the 5-year flood as an index to evaluate the magnitudes of 
floods of equal frequency. Within similar physiographic regions, floods of 
equal frequency have quite similar magnitudes when expressed as ratios to an 
index flood. Thus it is possible to use the sum total of flood experience in a 
region to develop a composite dimensionless frequency curve and to obtain 
more reliable estimates of the frequencies of major floods than can be ob- 
tained by a single record. Dalrymple, in another paper in this symposium, 
has used the ratio to the mean annual flood to indicate the relative magnitudes 
of the 1955 floods. 

Correlations between index floods and size of drainage area are improved 
by introducing factors such as the lagtime for the unit -hydrog raph(8) ora 
factor which is a measure of average basin slope.(7) 

At least two conclusions can be drawn from this paper. One is that the 
conservative meteorologic, hydrologic, and hydraulic design criteria(9,10) 
that have been used in the past by the Federal engineering organizations are 
now justified. Modifications to the criteria will be inevitable as operating 
experience is accumulated but the basic conservatism must be retained. The 
second conclusion is that experience and judgment in hydraulic engineering 
cannot be replaced by stereotyped design procedures. 


a. Proc. Paper 1663, June, 1958, by Elliot F. Childs. 
1. Prof. of Hydr. Eng., Massachusetts Inst. of Technology, Cambridge, Mass. 
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H. ALDEN FOSTER,! M. ASCE.—Mr. Childs’ paper brings out very clear- 
ly the fact that hydrological knowledge is never complete or final. The basic 
data for design of hydraulic structures need to be continually revised to con- 
form with new experience. Nature always demands the last word in such 
matters. 

The writer will confine his comments to the portions of the paper dealing 
with questions of probability or frequency of floods. 

The author refers to the “probable maximum precipitation” which is used 
as a basis for deriving the “maximum possible flood”. This seems to involve 
some confusion of ideas, although the latter expression has been commonly 
used in the past by the Corps of Engineers. It does not seem advisable to 
place any finite limit to future flood magnitudes. In any case, the “probable 


maximum precipitation” should be used to estimate the “probable maximum 
flood”. 


1. Prin. Assoc., Parsons, Brinckerhoff, Hall & Macdonald, New York, N. Y. 
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The difficulties of combining extraordinary floods with the “run-of-the- 
mill” variety when preparing a probability curve from a particular flood 
record were recognized by Allen Hazen nearly 30 years aga. (1) It is closely 
related to the question of “sampling errors” and ‘confidence limits”.(2,3) 

Any sample record of 25 years or more is likely to have one or more extreme 
values well outside the general trend of the rest of the sample. The writer 
has been unable to develop any method by which the floods caused by extra- 
ordinary storms, such as hurricanes, could be separated from the record, 
because some of the smaller floods may have been produced by hurricane 
storms at a considerable distance from the gaging station which would cause 
excessive floods in other locations. 

In preparing a flood-probability curve, one or two floods that obviously do 
not conform with the general trend should be ignored in selecting the curve- 
of-best-fit. This curve should only be used for such purposes as determining 
the mean annual flood damage where omission of the extreme floods would 
have relatively small effect on the results because of their obviously small 
probability. Of course the hurricane floods should not be ignored in spillway 
design, but here they should be supplemented with estimates of the “maximum 
probable flood”. The probability of the latter is necessarily small, but we 
are not concerned with whether it is the 1,000-year flood or the 100,000-year 
flood; as long as we are convinced that it may occur sometime, or even be 
exceeded, we cannot afford to ignore it. 

These remarks do not include the so-called “historical floods”. It should 
be generally possible to include them in the probability study by the methods 
described in the ASCE “Hydrology Handbook”, page 102ff. See also (4). 

The following comments are quoted from a report(5) on the Northeastern 
Floods of 1955 prepared under the supervision of the writer: 


“Climatic Changes. —The great increase in precipitation and stream 
discharge in the 1955 floods as compared with previously recorded 
storms indicates that the storms of August 1955 were of a very ex- 
ceptional nature. From a statistical standpoint, it appears that they 
would have a probability of occurrence much smaller than would be 
indicated by theoretical studies based on previously available records. 
Two alternative inferences may be drawn from this situation: 


(1) The 1955 storms were the result of atmospheric conditions of a 
most unusual character, following the simultaneous occurrence 
of a number of unfavorable controlling factors that normally 
would not be expected to occur together; or 

(2) There have been changes in general climatic or atmospheric con- 
ditions in recent years that tend to increase the frequency and 
intensity of storms of this type along the North Atlantic seaboard. 


“If the first inference is accepted, it would be possible to make a 
reasonable estimate of the probability of occurrence of such storms in 
the future, for use in determining the effect of the 1955 floods on the 
mean annual flood damage to be expected on rivers in the area under 
consideration. If the second inference holds, however, any estimates of 
flood probability based on previous records would have to be revised to 
include the effect of the assumed changes in meteorological conditions. 

. .. The recent apparent changes in climate affecting the paths of hurri- 
canes along the North Atlantic seaboard cover too short a time period to 
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justify the assumption of any long-term major climatic changes in this 
region. Nevertheless, the possibility of such cyclic variations in climate 
should not be neglected in estimating flood probabilities. 

“Probability of Occurrence of 1955 Floods.—If the question of cli- 
matic change is eliminated from consideration, the excessive magnitude 
of the 1955 storms must be treated as due to an unusual combination of 
most unfavorable controlling factors, but excluding any additional 
factors that have not been effective during the period of record. Asa 
result, if we have prepared a suitable probability curve representing the 
previously recorded floods, the 1955 floods should be located on the 
same curve. But since these floods are so much larger than any previ- 
ously observed, such a method of plotting will result in assigning a very 
low value of probability to them. If the present hurricane situation is 
assumed to be the result of climatic changes, the question of determin- 
ing the probability of occurrence of the 1955 floods becomes much more 
difficult. Such a change would have to be assumed as commencing not 
earlier than 1938, the year when the first severe hurricane in more than 
100 years struck the Long Island coastline. The intervening period of 
17 years (between 1938 and 1955) is much too short to permit any reason- 
able estimate of probability for the recent storms. All of this brings 
out the difficulties of estimating mean annual flood losses in any locality 
by statistical methods, even if the problem of errors of sampling could 
be satisfactorily taken care of.” 


Although the 1955 floods in many places were exceptionally large in com- 


parison with previous recorded floods, their effect on the determination of 
the Mean Annual Flood Loss for any particular property may be relatively 
small. Studies referred to above indicated that, for a number of properties 
in the City of Waterbury, Conn., the addition of the 1955 floods to the record 
available for 1930-1953 would have increased the estimate of mean annual 
flood losses for those properties by less than 3 per cent. 
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